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Tyrh"''^"' Field 

The present invention relates to screening and production of compound libraries for 
20 drug development. 

packpround 

Conventional methods to identify leads for drug development involve primary 
screening of compound libraries for activity "hits," foUowed by secondary screening to 
reduce the number of primary hits to a congeneric series of optimal leads for drug 
25 development. 
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The compound libraries, such as synthetic (e.g., combinatorial) and namral product 
(e.g., biological preparations and extracts) libraries, vary in size and complexity, ranging 
from hundreds, thousands, to millions or more of related or diverse compounds. The smaller 
libraries usually are well defined and each of the compounds frequently are contained in a 
separate storage or test vessel (e.g., dry or liquid form of the compound residing in a well of a 
multi-well storage or test plate with other members of the library). Larger libraries often are 
less well defined and typically contain mixtures or pools of compounds per vessel. For 
libraries containing pools of compounds, where an activity hit resides in one pool compared 
to the next, deconvolution and chemical analyses are typically performed in parallel to isolate 
and characterize the compound(s) responsible for the observed activity. Information gleaned 
from the initial screening and testing process also is used for subsequent rounds of analog 
synthesis (analog/focused libraries) and convergent screening and testing of particular 
analogs (i.e., iterative process). Computer-implemented theoretical or virtual compound 
libraries also provide a repository fixjm which activity hits are selected for known or 
predicted structure-activity relationships. 

Primary activity screening of compound libraries is based on selection of compounds 
that directiy or indirectly interact with a specific biological receptor(s) (i.e., receptor- 
dependent activity screening). Isolated receptors and cells expressing single or combinations 
of receptors chosen to mimic a particular biological system or disease state generally provide 
the context of an assay for receptor-dependent activity screening. For high-throughput 
screening of larger libraries, automated systems utilizing multi-well arrays representing 
isolated receptors or cells that express them are the standard. 

The driving force behind receptor-dependent activity screening as the primary 
approach for sifting through compound libraries is simple. Drugs 
(phamacological/toxicological agent) elicit a phannacological response through interaction 
with one or more biological receptors (drug/receptor-specific interaction). Thus, compounds 
that interact with a particular receptor or combination of receptors are presumed to be the 
most promising candidates for exhibiting some mutual activity in vivo and tiius targets for 
secondary screening. 
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Compounds identified from a primao' screen are then subjected to successively more 
focused and quantitative rounds of screening and validation to eliminate false positives and 
identify those exhibiting optimal biological activity against a target receptor(s) in an in vi^>o 
setting. This typically involves a combination of physiochemical and biological testing, 
including structural characterization and biological studies using cells, tissues and animals. 
Compounds with the most promising biological activity are selected as leads for drug 
development. 

Drug development involves scale up and detailed toxicity, pharmacodynamic and 
pharmacokinetic studies that are performed to characterize pharmacological efficacy. These 
) studies are conducted not only to gauge whether a test compound has activity in an in vivo 
setting, but also to examine bioavailability to assess possible route of administration, delivery 
formulations and the amount of a test compound necessary over time to produce a therapeutic 
effect with Uttle or acceptable side effects. A variety of cell, tissue and animal model assays 
typically are employed for such studies. A handful of compounds (e.g., 5-10) that pass these 
5 tests are then tested in scaled up animal studies for further characterization. A lead drug 
compound with the most promising results in animal studies is then tested in humans in 
clinical trials. 

Pharmacokinetic studies are conducted to characterize the time-dependent 
concentration of a test compound in the body, which collectively depends on absorption. 
0 distribution, metabolism and elimination (ADME) of the compound following administration. 
For instance, in order to reach the site of action, a lead drug compound that is administered to 
a subject must first be absorbed across epithehal barriers, usually by passive diffusion and/or 
active uptake, into the systemic circulation. In the case of intravascular administration, 
absorption is instantaneous and complete. However, all other routes of administration 
25 involve an absorption step with the potential that only a fraction of the administered 
compound may be absorbed into systemic circulation. 

Systemic blood then delivers the compound to cells and tissues in the body, where the 
likely receptor/site of action resides, but various parallel processes compete for the 
compound. The compound may reversibly bind with proteins (albumin, al-acid 
30 glycoprotein) in plasma, or in some instances with tissue proteins. This is important since an 
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unbound compound is typically the form taken up by cells and tissues. These processes 
determine distribution of the compound. 

In a process referred to as excretion or elimination, organs such as the kidney, lung 
and Uver are able to remove an unchanged lead drug compound from systemic circulation. 
Alternatively, the compound may be metabolized by enzymes frequently localized in all 
tissues but mainly in the liver. Such metabolism produces metabolites that are chemically 
different from the administered compound and generally are more readily excreted from the 
body (reduced lipid solubility). Often the phaimacologicaytoxicological activity of a 
metabolite is reduced compared to that of the parent compound. 

Thus while a lead or collection of lead drug compounds may continue to exhibit 
promising activity profiles early in the drug development process, most fail to make it as a 
1 drug product because of poor bioavailability discovered in animals, or worse poor 

Hi bioavailability not discovered until human clinical trials (e.g.. gancyclovir). This 

5 unacceptably high and expensive failure rate can be attributed in large part to the biased 

nature of activity-based screening to identify primary hits ultimately used as lead drug 
candidates. For instance, activity screening is pursued from the mindset that the greater and 
more specific the compound-receptor interaction/activity, the more potent a compound, and 
thus the smaller the dose required and consequent lower potential for toxic side-effects, as 
weU as cheaper product produced and sold. However, a potent compound exhibiting poor 
bioavailability might require a higher dose than a less potent compound exhibiting superior 
bioavailability; this less potent compound also may exhibit reduced dose related toxicity. 
Therefore, the majority of activity levels do not result in drug products. 

Receptor-dependent screening and testing also provides little to no information as to 
the probable route of administration for an activity hit. As an example, a test compound 
25 selected for activity may ultimately require intravenous administration, which is a less 
preferred route of administration. Here again a different less potent compound overlooked or 
discarded from an activity screen for lower potency may have been a good candidate for a 
preferred extravascular form of administration (e.g.. oral). An oral form would be cheaper to 
administer even if administered at a higher dose to compensate for lower potency. 
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The dogmatic process of screening compound libraries first by receptor activity 
likewise reduces the value of the libraries themselves. Newly obtained or previously 
screened compounds having true therapeutic potential due to superior bioavailabiUty 
properties are likely never to make it into the drug development pipeline if they fail to pass 
the primary activity screening process. Also valuable physical and chemical information 
from compounds otherwise possessing good bioavailability profiles that are discarded or 
overlooked for having less than some preferred activity level will be lost and unavailable for 
future development of structurally related activity leads or synthesis of new libraries. 

Accordingly, a need exists for identifying compounds that exhibit desired 
pharmacokinetic properties before the drug development process, as well as guidance for 
future synthesis. The present invention provides an unprecedented and counterintuitive 
approach to address these and other needs. 



Barlow et al. (WO 9716717) disclose a robotic system for automated in vitro 
15 measurement of cell permeability. Pidgeon et al. (J. Med. Chem. (1995) 38:590-594) disclose 
immobilized artificial membranes for permeability assays. Minth et al., (Eur. J. Cell. Biol. 
(1992) 57:132-137) disclose apparatus forperfiision cell cultures and in vitro assays. Various 
pharmacokinetic models of oral drug absorption are disclosed in Grass, G. (Advanced Drug 
Delivery Reviews (1997) 23:199-219); Amidon et al., {Pharm. Res. (1988) 5:651-654); 
20 Chiou. W.L., {Int. J. Clin. Pharmacol. Ther., (1994) 32:474^82); Chiou. W.L., {Biopharm. 
Drug Dispos., (1995) 16:71-75); Dressman et al., {J. Phami. Sci., (1985) 74:588-589); 
Lennemas et al., {J Pharm. Pharmacol., (1997) 49:682-686); Levet-Trafit et al.. {Life 
Sciences., (1996) 58:PL359-63); Sinko et al.. {Pharm. Res., (1991) 8:979-988); and Soria et 
al... {Biopharm. Drug Dispos., (1996) 17:817-818)). Grass et al., (Investigative 
25 Ophthamology & Vis. Sci. (1993) 34(7):225 1-2259) disclosure simulation model to predict 
aqueous humor and plasma phannacokinetics of oculatory appUed drugs. Audus et al., 
{Pharm. Res. (1990) 7(5):435-451) reviews epithelial and endotheUal cell models for drug 
transport and metabolism. 
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f;,l| UMARV o r THFi "■■Vf NTION 

The present tavemion reUtes to a method of screening and producing compound 
„,™es selected for absorption, and optionaily, one or more addUiona, propemes. Novel 
pr<.uced by the method ofthe invention a,so a. provrded. The medto^^ ts read,. 
5 adapted for high-throughput screening and production of compound .ibranesop.un.ed o 

JrpUon. The method and libraries of the invention can be utilized for preparahon of a 
medicament tor use in the treatment of a mammal. 

The method tnvolves sc^ng a first compound library or portion thereof which 
comprises a plurality of test samples containing isolated compounds and/or isolated mtxtures 
„ of Impounds per test sample. Screening U performed by: (i) generating an .„ v.vo 
absorption profile for each tea, sample fiom inihal d»e or amount and . v,.o b,oava.lab,l ty 
^Lprising permeability and solubility daia, and optionally dtssolution rate and transfer 
mechanism data, v^here an absorption profile is characterized by one or more of 
concentrauon, rate and extent of transfer of atest sample across a physio,og,c« fr^m 
= ,5 ote site of administration to a selected samphng site of a mammalian system of mterest,(n 

g selecUng compound having a desired absorption profile; (iii) and producmg a second 

S compound Ubrary comprising the selected compounds; and (iv, optionally repeattng «eps (0 

through(iii)oneormorerimes,whereacompoundUbraryselectedforabsorpUon,sobtamed. 

The present invention also provides a method for genentting an m vivo absorption 
20 profile. Tins method involves providing .>, ^ro bioavailability daU for a test s^ple of a 
compound library as input data to a computer-hnplemented phannacolcinetic too. (PK tool) o 
to invenhon. The PK too. includes as computer-readable component; (a) an mput/output 
system suit*,., for dau input and d«a o«pu.; (b) a aimulafion engine; and (c) a stmulanon 
model characterized byamulfi-compartmentphysiological model ofamammahansys^^^^ 

2S interest comprising a barri« to absorption that is based on . selected route of a— t^^ 
The input/output system, simulation engtae and simularion model are capable of «o,^ 
„geth« to carry on, the s.^ otreceivh,g as input data, initia. dose of a test sample the 
site of administr^ion and in n,ro bioavailability dau including one or more of penne*,hty. 
solubiUty, dissolution rate, and transfer mechamsm data, and generating as output da^ a 
30 staulati in Vvo absorption profile for each test sample that reflects rate, extent and/or tota, 
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concentration of the test sample at a given a sampling site located on the other side of the 
barrier to the site of administration. 

Subsequent rounds of screening according to the method of the invention provide new 
secondary compound libraries that are increasingly optimized for bioavailabihty. The 
libraries of the invention may be utilized to generate additional sub-libraries by screening. 
Accordingly, libraries produced by the method of the invention increase the chance of 
identifying compounds having desired in vivo pharmacological activity for a selected route of 
administration. 



nVFTNITIONS 

Absorption: Process by which a compound transfers across a physiological barrier as a 
fimction of time and initial concentration. Amount or concentration of the compound on the 
external and/or internal side of the barrier is a fimction of transfer rate and extent, and may 
range from zero to unity. 

BioavailabiUty: Fraction of an administered dose of a compound that reaches the sampling 
site and/or site of action. May range &om zero to unity. Can be assessed as a fimction of 
time. 

Compound: Chemical entity. 

Compound Library: A collection of two or more isolated compounds, pools of compounds, 
or combinations thereof. Examples include natural, synthetic and synthetic combinatorial 
compound libraries. May include computer-readable compound files. 

Computer Readable Medium: Medium for storing, retrieving and/or manipulating 
information using a computer. Includes optical, digital, magnetic mediums and the like; 
examples include portable computer diskette, CD-ROMs, hard drive on computer etc. 
Includes remote access mediums; examples include internet or intranet systems. Permits 
temporary or permanent data storage, access and manipulation. 
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Data: Experimentally collected and/or predicted variables. May include dependent and 
independent variables. 

Dissolution: Process by which a compound becomes dissolved in a solvent. 
Input/Output System: Provides a user interface between the user and a computer system. 

5 Permeability: Ability of a physiological barrier to permit passage of a particular substance. 
Refers to the concentration-dependent or concentration-independent rate of transport (flux), 
and collectively reflects the effects of characteristics such as molecular size, charge, partition 
coefficient and stability of a compound on transport. Permeability is substance and bamer 
specific. 

10 Physiologic Pharmacokinetic Model: Mathematical model describing movement and 
disposition of a compound in the body of a mammal or an anatomical part of the body based 
on phannacokinetics and physiology. 

Primary Compound Library: Compound library having compounds not yet screened and 
selected for (i) absorption, or (ii) absorption and one or more additional bioavailabihty 
15 properties. 

secondary Compound Library: Compound library derived from a primary compound 
Ubrary having compounds screened and selected for one or more particular properties. 

Simulation Engine: Computer-implemented instrument that simulates behavior of a system 
using an approximate mathematical model of the system. Combines mathematical model 
20 with user input variables to simulate or predict how the system behaves. May include log.c 
components, such as system control statements. 

SolubiUty: Property of being soluble; relative capability of being dissolved. 

Transport Mechanism: Hie mechanism by which a compound passes a physiological 
barrier of tissue or cells. Includes four basic categories of transport: passive paracellular, 
25 passive transcellular, carrier-mediated influx, and carrier-mediated efflux. 
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Figure 1 shows schematic of method of the invention to generate in vitro bioavailability data 
for screening a primary compound library by absorption parameters. 

Figure 2 shows schematic of method of the invention for selecting sampling site relative to 
administration site and barrier to absorption. 

Figure 3 shows schematic of method of the invention to generate an absorption compound 
library from in vitro bioavailability data. 

Figure 4 is a high level INPUT/PROCESS/OUTTUT diagram of the PK tool of the 
invention. 

Figure 5 is a high level flow chart and structure chart of the PK tool and method of the 
invention. 

Figure 6 is a graphical diagram illustrating generic compartment-flow simulation model and 
15 exemplary symbolic relationships among comparmients, flow regulators, converters and 
input links. 

Figure 7 is a key for Figure 6. 

Figure 8 is a graphical diagram illustrating generic pharmacokinetic first-order two- 
compartment open plasma model for intravenous injection. D is total drug. V is apparent 
volume of distribution, and C is drug concentration for either plasma (p) or Ussue (t). kl2 
and k21 represent first-order rate transfer constants for movement of drug from compartment 
1 to compartment 2 (kl2) and from compartment 2 to compartment 1 (k21). klO represents 
first-order rate transfer constant for movement (elimination) of drug from compartment 1 to 
compartment 0. 
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FigurH rr^Uca. compa«.=n.-flo. diagram iUu=.ating plasma simuladon mode, 
!; S aTd .xemplary reWonsMps amo„, compa«ne„., flow regula.0., co„ve«.. 

and input links. 

5 F,„re 10 show. schemaHc of a method ot *e i— fo, development of a„ im„a. 
physiologic-hased simulation model for PK tool and meftod of the invenhon. 

^r. n shows schematic of a method of the invenhon for development of a physiologic- 
based simulation model having selecdvely op,imi«d adiustttent parameters. 
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^,ur. 12 shows graph^cal co.part.ent-flow d.agran, iUustraUng the .ass-volume Gl tract 
simulation model of the invention linked to a training/validation plasma model. 

Figure 13 illustrates compartment, flow regulator and converter components for the mass- 
15 volume GI tract simulation model of the invention. 

Figure 14 illustrates structural relationship among compartment and flow regulator 
components for the mass-volume GI tract simulation model of the invention. 

20 Figure 15 illustrates stmctural relationship among flow regulator and converter components 
for the mass-volume GI tract simulation model of the invenuon. 

^r,«ntc fnr the mass-volume GI tract simulation model of 
Figure 16 illustrates converter components for the mass voiun 

the invention. 

" Figure .7 compares plasma oonceatration profiles dedved clinical smdies of 

gancyclovir arri simulation using volume GI tract simulation model of the mvenuon. 

IS compares plasma concentration profiles derived fiom clinical studies of 
30 ^cyclovirandshnulation„singmass-volun.eG..ractsimula«onmode,of.he,m,ennon 
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Hgure 19 *ows graphical compartm=n.-flow diagram illus^ing tire v,vo da,a analysis- 
processing IV/PO PK model (intravenous/oral administration) of the mvenhon. 

Flgare 20 shows schematic of method for development of initial integrated physiologic- 
5 based GI tract simulation model of PK tool and method of the invention. 

F,gu« 21 shoves graphical oompartment-now diagram illustrating the GI tract fluid transit 
model component of the PK tool and meUrod of the invention. 

0 Figure 22 shows graphical comparanen.-flow diagram illustrating the GI tract solubility- 
dissolution model component of the PK tool and method of the invention. 

Figure 23 shows graphical comparmren.-flow diagram illustrating the GI tract absorption 
model component of the PK tool and method of the invention. 

Figure 24 shows graphical compaflment-flow diagram illustrating integndon of dre GI tract 
flnid transit model, solubility-dissolution model, and absorption model components for one 
GI segment of the PK tool and method of die invention 

20 Figure 25 shows graphical compar^nent-Bow diagram illushradng integrated GI tract 
simulation model components (without converters or input link connectors) of the PK tool 

and method of the invention. 

Flgor. 26 shows graphical compartment-flow diagram illustrating integrated GI tract 
25 simulation model components (with convoters a.^ input linlt connectors) of the PK tool and 

method of the invention. 

Figure 27 shows schematic of method for development of selectively optimized adjusttnent 
p^eters and for optumzation of the inte^ed physiologic-based GI tract simulation 
30 model of PK tool and method of the invention. 
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Figure 28 shows schematic of method for selection of model parameters for utilization in a 
given physiologic-based GI tract simulation model of PK tool and method of the invention. 

Figure 29 shows schematic of method for regional (segmental) calculation/estimation of 
5 permeability from one or more user input values for permeability of a given GI tract 
region/segment. Regional permeability (Pe) correlation based on input of Pe value for 
duodenum is illustrated. 

Figure 30 shows graphical converter diagram illustrating volume, surface area, dose, time 
10 and pH parameters and calculations for integrated GI tract simulation model components of 
the PK tool and method of the invention. 

Figure 31 shows graphical converter diagram illustrating GI tract transit time parameters and 
calculations for integrated GI tract simulation model components of the PK tool and method 
15 of the invention. 

Figure 32 shows graphical converter diagram illustrating GI tract permeability parameters 
and calculations for integrated GI tract simulation model components of the PK tool and 
method of the invention. 

20 

Figure 33 shows graphical converter diagram illustrating GI tract solubility parameters and 
calculations for integrated GI tract simulation model components of the PK tool and method 
of the invention. 

25 Figure 34 shows graphical converter diagram illustrating GI tract control release formulation 
parameters and calculations for integrated GI tract simulation model components of the PK 
tool and method of the invention. 

Figure 35 shows graphical compartment-converter diagram illustrating GI tract concentration 
30 parameters and calculations for integrated GI tract simulation model components of the PK 
tool and method of the invention. 



12 



PCT/US99/2n51 

WO 00/16231 

Figure 36 shows graphical comparm^ent-converter diagram illustrating GI tract d.sso ut.on 
parameters and calculations for integrated Gl tract simulation model components of the PK 
tool and method of the invention. 

5 Figure 37 shows graphical compartment-converter diagram illustrating GI tract output 
calculations for absorption for integrated GI tract simulation model components of the PK 
tool and method of the invention. 

Figure 38 shows graphical convmcr disgran, illusMting GI «ct ou^u. calculations for 
10 soluble mass absorption rat. (flux) for integrated GI tract simuUtion model components of 
the PK tool and method of the invention. 

Figure 39 shows graphical compartment-flow-convexter diagram illustrating GI tract output 
calculations for cumulative dissolution rate and amount for mtegrated Gl tract simulaUon 
15 model components of the PK tool and method of the invention. 

Figure 40 shows graphical compartment-flow-converter diagram illustratmg GI tract output 
calculations for cumulative control release fomiulation rate and amount for integrated GI 
tract simulation model components of the PK tool and method of the invention. 

20 J ^ 

Figure 41 illustrates database and mlebase compartment, flow regulator and converter 

components for the integrated physiologic-based GI tract simulation model of the mvention. 

Figure 42 illustrates structural relationship among compartment and flow regulator 
25 componentsfortheintegratedphysiologic-basedGItractsimulationmodeloftheinvenUon. 

Figure 43 iUustrates structural relationship among flow regulator and converter components 
for the integrated physiologic-based GI tract simulation model of the invention. 

30 Figure 44 illustrates structural relationship among converter components for the integrated 
physiologic-based Gl tract simulation model of the invention. 
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Figure 45 is a high level INPUT/PROCESS/OUTPUT diagram of the PK tool of the 
invention as presented to a user of the carrying out a method of the invention, with inputs 
provided by the user and outputs provided by the PK tool. 

5 Figure 46 illustrates a flow chart and structure chart of a subsystem of the PK tool and 
n^ethod of the invention for selection of a physiological GI tract model from a model 
database and a parameter database. 

Figure 47 is a flow chart and structure chart of the system of the PK tool and method of the 
10 invention. 

Figure 48 is a flow chart and structure chart of a menu of the system of the PK tool and 
method of the invention. 

15 Figure 49 illustrates correlation of extent of absorption for fraction of the dose absorbed in 
portal vein (FDp). as predicted using physiologicbased GI tract simulation model and PK 
tool of the invention, to FDp derived fiom human clinical data for 12 compounds. 

Figure 50 illustrates correlation of mte of absorption for fraction of the dose absorbed in 
20 portal vein (FDp). as predicted using integrated physiologic-based GI tract simulation model 
and PK tool of the invention, to FDp derived from human clinical data for 12 compounds. 

Figure 51 compares plasma levels as predicted using integrated physiologic-based GI tract 
simulation model and PK tool of the invention, to plasma levels derived fiom human clinical 
25 data for a test compound. 

Figure 52 compares plasma levels as predicted using integrated physiologic-based GI tract 
simulation model and PK tool of the invention, to plasma levels derived from human clinical 
data for a test compound. 
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Fiaura 53 compares ptema levels as predicted using ■■..egraled physiologic-based GI m>c. 
simulation model and PK tool of the invennon, to plasma levels denved from human chmcal 

data for a test compound. 

Figure 54 shows high level INPUT/PROCESS/OUTPUT diagram of the PK tool of the 
invention for SAR/QSAR and CAD/CAE compound design and synthesis. 
Figure 55 shows high level flow and structure chart for screening method of the invention 
utilizing the PK tool and method of the invention. 

^^cf^ppT^nN OF SPFCTFTr FMBOPTMENTS 

The present invention relates to a method for screening compound libraries by 
absorption, or absorption and one or more additional properties. The invention also relates to 
compound libraries produced by the method of the invention. 

The method of the invention involves screening a first compound library or portion 
S thereof by absorption, where the compound Ubrary or portion thereof includes a plurahty of 

i 15 test samples containing isolated compounds and/or isolated mixtures of compounds per test 

sample. Screening of the first compound library or portion thereof is perfomied by: (1) 
generating an in vivo absorption profile ftom in vitro bioavailability data for each test sample, 
where the absorption profile is based on a selected route of administration and samphng site 
of a mammalian system of interest; and (2) selecting test samples having a desired absorption 
20 profile compared to others. 

•nxe in vivo absorption profile is characterized by absorption rate, extent of absorption, 
and/or concentration of a test sample relative to a selected site of administration and a 
selected sampling site of the mammalian system of interest, i.e.. rate and/or extent of transfer 
of a test sample from an external site (e.g.. apical) acn,ss a physiological bamer (e.g., 
epithelium) to an internal site (e.g.. basolateral) of that barrier. THis can include prediction of 
rate extent, and/or concentration of a test sample at the site of action when the samplmg site 
is t^e site of action. Transfer rate and/or extent are generated from initial dose data (e.g., 
amount) for the test sample and in vitro derived bioavailability data including permeabUity 
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and solubility data, and optionally dissolution rate and transport mechanism data (i.e., passive 
paracellular, passive transcellular, carrier-mediated influx, carrier-mediated efflux) for the 
test sample. Solubility and dissolution rate are interrelated and effect the ability of the 
compound to be solubilized at a rate sufficient for absorption to occur across a particular 
membrane. Penneability refers to the concentration-dependent or concentration-independent 
rate of transport (flux), and collectively reflects the effect of molecular size, charge, partition 
coefficient and stability of a compound on absorption for a particular physiological barrier, 
where the physiological barrier(s) depends on the selected route of administration. Molecular 
size, charge and partition coefficient determines in large part whether a compound is 
trar^orted via a paracellular or transcellular mechanism. Stability is a general feature that 
relates to whether the compound remains intact long enough to be absorbed. Together, in 
vitro derived solubility and permeability data, and optionally dissolution rate and transport 
mechanism data, are primary bioavailabihty factors utilized to prepare an absorption profile 
for a test sample of interest. 

An in vivo absorption profile may be generated by any number of pharmacokinetic 
techniques. The preferred method for generating an in vivo absorption profile is by providing 
initial dose and in vitro bioavailability data for each test sample as input data to a computer- 
implemented pharmacokinetic tool (PK tool) of the invention as described herein. The PK 
tool then generates as ou^ut a simulated in vivo absorption profile. This aspect of the 
invention provides a rapid and accurate way to predict in vivo absorption of a library of 
compounds from in vitro data, or absorption in one type of mammal (e.g.. rabbit) to a 
different type of mammal (e.g., human). This is important since in vitro absorption data 
cannot be used directly to predict absorption in vivo, nor can in vivo data firom one type of 
mammal be used directly to predict absorption in a second different type of mammal. 
Moreover, when a diverse set of compounds reside within a library to be screened, then 
conventional methods of utilizing in vitro data to predict in vivo absoiption, or in vivo 
absorption data from one type of mammal to the next wiU have an unacceptably high failure 
rate, i.e., false positives and false negatives relative to absorption. The PK tool of the 
invention also is readily adaptable for both high-throughput and high-resolution screening 
formats, and provides information necessary for ranking compounds by bioavailability 
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parameters comprising absorption where statistical correlation and other prediction schemes 
fail. 

The PK tool includes as computer-readable components, an input/output system, a 
simulation model, and a simulation engine. The input/output system may be any computer- 
implemented system suitable for data input and data output and operable interaction with the 
simulation engine and simulation model. The simulation engine includes a differential 
equation solver, and optionally a system control statement module. This includes various 
computer-readable algorithms for numerical iteration of mathematical equations over interval 
dt and for processing rules, scenarios, pattern matching and the like that direct the simulation. 
The simulation model corresponds to a physiologic-based multi-compartment model of a 
mammalian system of interest, where the mammalian system represents a barrier to 
absorption that is based on a selected route of administration, i.e., the location at which the 
compound is introduced to a mammal. More particularly, the physiologic-based simulation 
model of the PK tool of the invention is a mathematical model comprising as operably linked 
components: (i) differential equations for calculating rate, extent and concentration of a test 
sample for one or more physiological segments of the mammal system of interest; and (ii) 
initial parameter values for the differential equations corresponding to physiological 
parameters, and optionally one or more selectively optimized adjustment parameters, and 
optionally one or more regional correlation parameters, for one or more physiological 
segments of the mammal system of interest; and optionally (iii) control statement rules for 
one or more of absorption, permeability, solubility, dissolution, concentration, and 
mathematical error correction, for one or more physiological segments of the mammal system 
of interest. The simulation model also may include one or more smoothing functions that 
facilitate calculation of transitional parameter values occurring between one or more of the 
physiological segments. 

The differential equations of a selected simulation model of a mammalian system of 
interest describe the rate processes of absorption, and optionally other events, of that model, 
which in turn describe drug concentrations in the system as a function of time. (See, e.g., 
Shargel et al.. Applied Biopharmaceutics and Pharmacokinetics, Appelton & Lange, East 
Norwalk, Conneticut, 1993). Thus, the differential equations are selected for a particular 
model. 
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The initial physiological parameter values of a given simulation model can be 
generated de novo or obtained from existing sources including the literature. The selectively 
optimized adjustment parameter values of a given simulation model of the invention 
represent regression or stochastic analysis derived values that are used as constants for one or 
more independent parameters of the model. In particular, the selectively optimized 
adjustment parameter values are obtainable by using a stepwise fitting and selection process 
that employs regression- or stochastic-based curve-fitting algorithms to simultaneously 
estimate the change required in a value assigned to an initial absorption parameter of the 
model in order to change an output variable. The input variables utilized for fitting include a 
combination of /« vitro data (e.g., permeability, solubility) and in vivo pharmacokinetic data 
(e.g., fraction of dose absorbed, plasma levels) for a compound test set having compounds 
exhibiting a diverse range of m vivo absorption properties. Thus, the input variables used for 
regression- or stochastic-based fitting are derived from (a) a first data source corresponding 
to the mammalian system of interest (e.g., in vivo pharmacokinetic data from human for the 
compound test set), and (b) a second data source corresponding to a system other than the 
mammalian system of interest (e.g., in vitro solubility data and in vitro permeabiHty data 
from rabbit tissue for the compound test set). A fitted adjustment parameter value for a given 
independent parameter is then selected that, when supplied as a constant in the model, 
permits correlation of one or more of the input variables fi^m the first data source to one or 
more input variables from the second data source. The process is repeated one or more times 
for one or more additional independent parameters of the simulation model until deviation of 
the correlation is minimized. These "selectively optimized" adjustment parameters are then 
provided to a given simulation model as constants or ranges of constants or fimctions that 
modify the underlying equations of the model. The selectively optimized adjustment 
parameters facilitate accurate correlation of in vitro data derived from a particular type of 
assay corresponding to the second data source (e.g., Caco-2 cells, segment-specific rabbit 
intestinal tissue sections etc.) to in vivo absorption for a mammalian system of interest 
corresponding to the first data source (e.g., segment-specific portions of the human GI tract) 
for diverse test sample data sets. Selectively optimized adjustment parameters also can be 
utilized to facilitate accurate correlation of in vivo data derived from a first species of 
mammal (e.g., rabbit) to a second species of mammal (e.g., human). 
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For a simulation model representing two or more anatomical segments of a given 
mammalian system, the model will preferably include regional correlation parameters. The 
regional correlation parameters permit estimation of a selected parameter value for a first 
segment of the mammalian system from correlation using a value of the selected parameter 
for a second segment of the mammalian system. The regional correlation parameters 
represent a collection of empirically derived values or selectively optimized adjustment 
parameter values for various segments of the mammalian system of interest, for example, 
permeability values. The regional (i.e., segmental) correlation is performed by logic function 
of the model, which when activated utilizes a function/transformation algorithm to estimate 
the parameter value for the second segment from (1) the corresponding regional correlation 
parameters, and (2) a user provided input value for the same parameter, but for a different 
segment. The regional correlation logic function of the model is activated when a user does 
not supply an input value for a particular parameter. For example, when a user of the PK tool. 
suppUes a single permeability value as input into a GI tract simulation model of the invention, 
such as a permeability value derived from Caco-2 cells that corresponds to colon, then 
regional permeabiHty correlation is performed by the PK tool to estimate permeabihty in the 
other GI tract segments, such as duodenum, jejunum, and ileum. 

The control statement rules include various logic elements utilized for providing 
guidance as to how a given simulation is to proceed. For instance, a control statement rule 
would include "IF ... THEN" production rules. An example of a production rule would be 
"IF solubility of compound is zero THEN absorption is zero." The production rules are based 
on rules of thumb (heuristics) and the Uke, and may be generated by correlation of parameters 
and simulation runs. Rules can be added, modified or removed to change how a simulation 
model responds to incoming data. 

The input/output system, simulation engine and simulation model of tiie PK tool are 
capable of working together to carry out the steps of (1) receiving as input data, tiie initial 
amount of a test sample at the site of administration and in vitro bioavailabihty data including 
one or more of permeability and solubility data, and optionally dissolution rate and transfer 
mechanism data; and (2) applying tiie simulation engine and ttie simulation model to generate 
as output data a simulated in vivo absorption profile for each test sample tiiat reflects rate, 
extent and/or concentration of tiie test sample at a given sampling site for a selected route of 
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administration in a mammalian system of interest. This includes uni- and multi-dimensional 
output profiles that collectively reflect parameters of absorption, which can be directly or 
indirectly uUUzed for characterizing in vivo absorption, as well as one or more additional 
bioavailability parameters including distribution, metabolism, elimination, and optionally 
toxicity. 

The in vitro data corresponding to test samples of the library can be empirically 
derived ftom experimental assays (e.g., physiochemical, cell or tissue assays) or theoretical 
data predicted from one or more other bioavailability parameters derived from the assays 
(e.g., calculated estimation) and/or predicted from molecular structure information (e.g., 
structure-property) where structural information is available. The preferred data for high- 
throughput screening is empirically derived cell culture-based in vitro data. Preferred data 
for high-resolution screening is empirically derived tissue-based in vitro data. In vivo derived 
mammal (animal, human) data may be employed for model development, training and/or 
validation purposes, as well as for predicting absorption in a first species of mammal from in 
vivo data derived from a second species of mammal. 

Test samples of the first library can have known or unknown biological activity, and 
may be derived from compound libraries including nattiral and/or synthetic compounds and 
pools, as well as compound files. Libraries for high-throughput screening may range up to 
the maximal library size, and are preferably screened in blocks. The number of compounds 
or compound pools per block is detennined by the user, and typically range from 1,000 to 
100,000 compounds per block. Preferred libraries for high-resolution screening range from 
1,000 to 10,000 compounds, more preferably bam 100 to 5000 compounds, and even more 
preferably 50 to 1000 compounds. Of course the actual number of compounds per library 
and screening can vary depending on the intended end use, and may employ a combination of 
high-throughput and high-resolution screening ^proaches. 

The selected routes of administration include enteral (e.g., buccal or sublingual, oral 
(PO), rectal (PR)), parenteral (e.g., intravascular, intravenous bolus, intravenous infusion, 
intramuscular, subcutaneous injection), inhalation and transdermal (percutaneous). The 
preferred route of administration according to the method of the invention is oral 
administration. The selected route of administration determines the type and/or source of 
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assay or structure-property parameters employed for obtaining a set of in vitro bioavailability 
data utilized for generating a simulated in vivo absorption profile. That is, artificial, cell or 
tissue preparations and the like derived from or representative of a physiological barrier to 
absorption for a selected route of administration are chosen to generate the relevant in vitro 
bioavailability data for use as input into the PK tool. For instance, in vitro bioavailability 
data for simulating fate of a test sample following oral administration can be based on cell 
culture and/or tissue assays that employ biological preparations derived from or 
representative of the gastrointestinal tract of a mammal of interest, e.g.. gastrointestinal 
epithelial cell preparations for permeability and transfer mechanism data, and 
physiological/anatomical fluid and admixing conditions corresponding to the relevant 
portions of the gastrointestinal tract for solubility and dissolution rate assays. Assays for 
collecting in vitro bioavailability data for specialized physiological barriers such as the blood 
brain bairier may initially assume intravascular delivery and thus instantaneous absorption as 
a first step. In this situation an assay is selected to generate in vitro bioavailability data 
relative to the blood brain bairier, which include for instance cell culture and/or tissue assays 
that employ biological preparations derived from or representative of the interface between 
systemic blood and the endothelial cells of the microvessels of the brain for a mammal of 
interest, e.g., blood-brain-bairier microvessel endothelial cell preparations for permeability 
and transfer mechanism data, and physiological/anatomical fluid and admixing conditions 
corresponding to the relevant portions of the blood membrane bairier for solubility and 
dissolution rate assays. A series of assays may be employed to collect in vitro bioavailabihty 
data for two or more bairiers to absorption. As an example, oral, hepatic, systemic and blood 
brain bairier assays may be utilized to obtain in vitro bioavailability data for screening 
compound Ubraries for orally deUvered compounds that target brain tissue. 

The sampling site relates to the point at which absorption parameters are evaluated for 
a test sample of interest. This is the site at which rate, extent and/or concentration of a test 
sample is determined relative to a selected site of administration, and is separated from the 
site of administration by at least one physiological bairier to absorption. For generating 
simulated absorption profiles, the sampling site preferably is separated from the site of 
administration by an individual primary barrier to absorption, which can be utiUzed to 
evaluate additional absorption events by secondary bairiers to absorption so as to sequentially 
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and collectively reflect the summation of absorption events at other sampling sites of interest. 
As an example, the sampling site selected for oral delivery may be the portal vein where the 
primary barrier to absorption is the gastrointestinal lumenal membrane, or systemic blood 
where a secondary barrier to systemic absorption is the liver after the test sample passes from 
5 the portal vein through the liver to systemic circulation. Thus the type of physiological 
barrier(s) residing between a site of administration and a sampling site reflects the type of 
assay(s) employed for generating the desired bioavailability data for use as input data into the 
PK tool of the invention. 

As the selected route of administration determines the barrier(s) to absorption and the 
10 physiological parameters that affect absorption events following administration, it also 
determines the physiologic-based pharmacokinetic simulation model employed in the PK tool 
for generation of the simulated in vivo absorption profile. By way of example, if the 
proposed route of administration is oral, then a primary barrier to absorption is the lumenal 
membrane of the gastrointestinal tract, and a secondary event affecting systemic 
15 bioavailability is first pass metabolism by the liver. Thus, a given simulation model and its 
associated parameters for simulating the fate of a compound selected for oral delivery is 
chosen to represent this scenario. The model would include therefore relevant components of 
the gastrointestinal tract for administration and absorption (i.e., stomach, duodenum, jejunum, 
ileum, and colon) and a primary sampling site (i.e., portal vein) from which to evaluate a 
H 20 primary absorption event. In this instance a secondary barrier to absorption for oral delivery 
is the liver and a secondary sampling site is systemic blood/plasma. This basic approach to 
choosing a physiologic-based pharmacokinetic model also applies to models employed to 
simulate absorption by target organs and the like, where a physiological barrier to absorption 
is the tissue and/or membrane separating systemic blood from the target organ, such as the 
25 blood brain barrier. In this situation if oral delivery is selected as the preferred route of 
administration for a compound targeting brain tissue, then a gastrointestinal tract model and 
blood brain barrier model may be implemented separately and/or combined through a 
complementary plasma component of the models for screening purposes. 

The physiological models are selected from a repository of delivery route-specific 
30 models stored in a memory, a database, or created de novo. Physiological models of the 
invention include those corresponding to common routes of administration or barriers to 
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absorption, such as oral (Gl tract), ocular (eye), transdermal (skin), rectal, intravenous, rectal, 
subcutaneous, respiratory (nasal, lung), blood brain barrier and the like. For constructing a 
model de novo, the basic approach is to identify and isolate a primary barrier to a specific 
absorption event from secondary events so that each barrier to absorption can be tested and 
validated in isolation. This involves selecting a site of administration that is separated from a 
sampling site by a primary physiological barrier to absorption and then building a 
developmental physiological model that incorporates rate process relations and limitations to 
describe the isolated absorption event. If desired, the secondary events can be added 
sequentially so that each additional layer of complexity to the model can be tested and 
) vahdated in isolation from other components of the initial model. 

Test samples selected for their predicted/simulated in vivo absorption profile from the 
first library can then be utilized to generate a secondary compound library, which may be 
physically separate from the originating (first) compound library. The secondary library also 
may be defined by simply cataloging test samples in the originating library by a descriptor(s) 
5 related to their respective absorption profiles. In particular, for selecting compounds having a 
desired absorption profile compared to others, the profiles are compared and the compounds 
ranked in order of optimal to minimal concentration, rate and/or extent of absorption at a 
sampling site of interest, and/or one or more of absorption parameters from the group of 
permeabiUty, solubility, dissolution rate, and transport mechanism. Ranking profiles can then 
0 be utilized to select compounds having a desired absorption profile, which include optimally 
absorbed compounds, but also may include those that are poorly absorbed. For instance, the 
selection can be based on choosing compounds that fall within a user defined window of 
absorption rate, extent of absorption and/or concentration at a selected sampling site. The 
user defined window can be based on nmges of concentration, rate and/or extent of 
25 absorption relative to a control or set of standard compounds having known absorption 
profiles for the selected route of administration. An example of a desired absorption profile 
includes compounds that exhibit moderate to optimal rate and/or extent of absorption for a 
particular barrier to absorption as compared to a control. 

The secondary library can then be subjected to fiirther rounds of bioavailability 
30 screening, including additional rounds of more focused absorption screening, as well as other 
screens that characterize the test samples by metabolism, toxicity, biological activity and the 
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like This process can be repeated one or more times to obtain libraries containing 
compounds that are increasingly optimized for absoiption, and optionally absorption and one 
or more other properties. 

Secondary libraries produced by the method of the invention are unique in that they 
contain compounds possessing a desired absorption profile as the common functional 
denominator for a selected route of administration, and thus substantially retain route-specific 
structural and functional diversity of activities and bioavailability residing in the parent 
library. By way of contrast, secondary libraries produced by activity screening are likely to 
represent reduced activity and bioavailability diversity, since compounds are selected solely 
to interact with specific receptor(s). Activity screens tend to select compounds with similar 
molecular structure, which tends to result in similar properties and therefore absorption 
profiles. Absorption screening selects for molecular properties and thus in more likely to 
maintain structural diversity. For instance, depending on the diversity of the first library, a 
secondary library produced according to the method of the invention will contain a variable 
composite of compounds or mixtures of compounds having a user defmed absorption profile, 
compounds having low to high activity against a particular receptor, as well as compounds 
that exhibit no activity against that particular receptor(s). These secondary libraries also will 
contain optimal diversity with regard to route-specific structure-bioavailability information, 
which is extremely useful for subsequent structure-based compound design and iterative 
synthesis of analog libraries and the like. 

Thus an advantage of screening compound libraries by absorption profiles is that the 
structural and fimctional activity and bioavailability diversity residing in the parent library for 
a selected route of administration is retained in the secondary library, although no activity or 
structural information is required to define the content of the new Ubrary. Another advantage 
is that a majority of all compounds residing in the newly created library will exhibit a user 
selected absorption profile, and thus the chance of identifying well absorbed leads for 
secondary screening by activity and other properties, and ultimately drug development, is 
improved. An additional advantage is that the libraries are reduced to a more manageable 
size for secondary screening, while retaining optimal activity diversity within the library in 
terms of a composite of inherent biological activities for a selected route of administration. 
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The present invention is significant and counterintuitive in that biological activity 
(i.e., receptor-interaction activity) is not required for the process, or for obtaining libraries 
optimized to contain compounds having improved in vivo pharmacological activity for a 
selected route of administration. Furthermore, screening according to the present invention 
increases the utility of compound libraries in general through identification of new lead drug 
compounds from: (1) libraries that have not been screened for activity; (2) libraries 
previously screened for activity that contain compounds failing to pass activity screens; and 
(3) Ubraries previously screened for activity that contain compounds failing to reach the drug 
development stage for lack of desired biological activity or for failing to work in a preferred 
route of administration or formulation. As can be appreciated, the method and libraries 
produced by the method of the invention increase the chances of finding better absorbed leads 
for drug development for a selected route of administration. The method of the invention 
also permits early identification of possible routes of administration for a lead compound 
identified by screening of an absorption library of the invention by biological activity. 

Compound Library: 

Compounds employed in the method of the invention may be firom physical 
compound libraries that contain natural and/or synthetic compounds and pools. Compound 
files (computer-readable compound representations and theoretical "virtual" libraries) also 
may serve as a reservoir from which to obtain compounds for screening. Examples of natural 
compound libraries include those that contain compounds obtained from biological 
preparations, such as from micioorganisms (viruses, bacteria), algae, lower plants (fimgi), 
higher plants, lower animals, mammals and the like. Examples of synthetic libraries include 
those that contain compounds generated using various synthetic chemistry techniques, such 
as solid and/or solution phase chemistries. Synthetic libraries produced by combinatorial 
chemistries are of particular interest. Techniques and sources for obtaining the compound 
libraries are well known, and new sources and chemistries are being developed at a rapid 
pace. The methods of the invention are applicable for any of these libraries. 

Compounds of a physical Ubrary are typically stored in multi-vessel storage and/or 
testing units, such as multi-well microtiter plates, as Uquids or solids. In particular, 
compounds of a given library can reside in pools containing mixtures of two or more 
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compounds, as extracts from biological preparations, and/or as isolated individual compounds 
per storage and/or test vessel of a multi-vessel storage and/or testing unit. The units may be 
stacked or stored in separate locations. As can be appreciated, however, individual vessels of 
a library do not have to be stored in the same physical location; for instance, they can be 
assigned to a particular library although they reside in a different storage unit and/or location. 
The same applies to compound files. 

As an example, compound libraries may be represented as a machine-readable 
compound file. This includes compound files stored and/or accessible on a computer- 
readable medium. Examples include optical and magnetic mediums and the Uke. Machine- 
readable compound library files are particularly useful, for example, when large 
combinatorial libraries are screened according to the method of the invention, and then the 
absorption profile and/or ranking information is used to catalogue the larger library 
electronically without the need to go back and create a separate physical library that refiects 
the new library. The compound file can then be accessed to retrieve, add or change 
information in the file to generate new libraries and/or proceed through iterative rounds of 
new screening. 

Historically, natural products are the most common source of new pharmaceutical and 
lead drug candidates. A natural product library will contain extracts of various naturally 
occurring substances. Common sources of the extracts are microbial sources such as various 
fimgi, bacteria, or algae. Plant extracts are also a common component of natural product 
libraries. The natural product libraries are readily producible using standard methods. 
Natural product libraries also are commercially available. For instance, natural product 
libraries can be obtained from a variety of commercial venders, such as Pan Laboratories 
(Bothell, WA) and MycoSearch (NC). 

In contrast to natural product Ubraries, synthetic compound libraries are composed of 
chemicals that are not necessarily natural. Synthetic compound libraries suitable for the 
present invention include libraries constructed de novo or obtained commercially. Any 
number of methods for constructing synthetic compound libraries can be utilized. Synthetic 
compound libraries also are commercially available from various sources. Examples of 
commercial sources for synthetic libraries include Maybridge Chemical Co. (TreviUet. 
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Cornwall. UK). Comgenex (Princeton, NJ), Brandon Associates (Merrimack, NH), and 
Microsource (New Milford, CT). In addition, a "rare chemical" library can be obtained from 
Aldrich Chemical Company, Inc. (Milwaukee, WI). 

Compounds from both natural product and synthetic compound libraries are readily 
modified through conventional chemical, physical, and biochemical methods (Blondelle et 
al.,TIBTech (1996) 14:60). 

Of particular interest are synthetic compound libraries produced through 
combinatorial chemistry, referred to as combinatorial libraries. Combinatorial chemistry is a 
technique of creating large libraries of diverse compounds through systematic and repetitive 
use of chemical building blocks or templates. The combinatorial libraries may be based on 
any number of templates or core molecules modified by addition of constituents. For 
instance, combinatorial libraries may include isolated or mixtures of peptide, oligonucleotide, 
and drug-like small molecules, or combinations thereof, such as those generated by pin 
technology and split-pool methods. Examples include peptides, peptidomimetics. cyclic 
peptides, constrained peptides, small non-peptide organics, nucleic acids, chiral and non- 
chiral compounds, drug-like small molecule libraries and the like. 

Combinatorial libraries can be made de novo or obtained commercially. Virtually an 
unlimited number of techniques can be used to create a combinatorial library, such as 
solution and solid phase chemistries. An advantage of solution-phase synthesis is that it 
capitalizes on the vast range of solution chemistry available in the chemical literature. Solid- 
phase synthesis is usefiil for facile purification and easier automation. For instance, 
combinatorial chemistry libraries can be produced by semi- of fiilly-automatic equipment 
available from various sources following manufactures protocols. Examples include 
Hewlett-Packard (Palo Alto. CA), Perkin-Elmer Applied Biosystems (Foster City. CA) or 
ChemTech (Louisville. KT). These types of equipment are compatible with a wide variety of 
coupUng chemistries. Combinatorial libraries also can be made-to-order or purchased 
commercially. Examples of commercial venders include Affymax (Palo Alto, CA). ArQule 
(Medford. MA), Helios Pharmaceuticals (Louisville, KT). Gryphon Sciences (So. San 
Francisco, CA). 
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Additional compound libraries include compound file libraries. Compound files are 
databases containing one-dimensional, two-dimensional or three-dimensional descriptions of 
chemical compounds. These libraries can be created to describe any chemical compound. 
For example, there are several commercial sources of both two-dimensional (e.g., from 
Maybridge Chemical Company. Bark Information Services, or the National Cancer Institute) 
and three-dimensional (e.g.. the Cambridge Small Molecule Library) compound files. 
Alternatively, these compound files can be easily created based on known chemical formulas 
and the basic chemical knowledge, such as standard atoms, bond angles and lengths and the 
like. 

Virtual compounds are of interest as they can be exploited using computer methods to 
design, select and iteratively refine compound libraries of small molecule drugs from 
structx^e-bioavailability data, structure-activity relationship (SAR) data or three-dimensional 
structural or pharmacophore models having known and/or predicted bioavailability profiles. 
An important feature of the process is the initial computer generation of veiy large virtual 
libraries of synthetically accessible compounds. These compounds are designed to explore 
specific structural features suggested from an input SAR or stmctural model. Any number of 
suitable computer systems can be utilized for this purpose. Virtual libraries typically contain 
100.000 to 1,000.000 compounds where each compound has a validated chemical synthesis 
pathway, and is characterized by a set of molecular descriptors. Computer codes are then 
used to select library sub-sets (100-1000 compounds) for rounds of automated synthesis and 
testing using high-throughput absorption screening according to the method of the invention. 
Any number of approaches can be employed for coding. Test data resulting from absorption 
screening for a given round of synthesis are interpreted by a computer selector code that is 
able to optimize multiple objectives simultaneously, in order to refine the properties of 
molecules selected in fiirther rounds of synthesis and testing. Such a screening approach can 
be implemented through a hierarchical computer system (e.g.. cUent-server, mainframe, real 
time system) that tracts compounds from virtual conception through to testing and compound 
file database registration. This approach is extremely useful for creating new structure- 
bioavailability data from analysis of both positive and negative absorption profiles for a test 
sample, e.g., structural differences that contribute to good versus poor absorption. 
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Assays for Generating Bioavailability Data: 

In viti'o bioavailability data utilized to generate an absorption profile for a test sample 
include permeability and solubility parameters, and optionally transport mechanism and 
dissolution parameters. Bioavailability data can be generated de novo following any number 
of techniques, or obtained from public or existing sources where available. The 
bioavailability data can be derived from chemical, and/or biological assays as well as 
theoretical predictions. By way of example, the in vitro assays may employ artificial 
(synthetic) or naturally occurring biological preparations. This includes chemical, cell and/or 
tissue preparations. Assays for generating in vitro bioavailability data involve screening a 
plurality of test samples containing isolated compounds and/or isolated mixtures of 
compounds per test sample in an assay characterized by measurement of (1) permeability and 
optionally transport mechanism for a test sample; and (2) solubility and optionally dissolution 
for a test sample. Methods and materials for performing the assays are based on the selected 
route of administration, the associated barrier(s) to absorption and proposed sampling site(s). 
For instance, if oral deHvery is proposed for simulation and an initial sampling site is selected 
to be the portal vein (so as to isolate gastrointestinal absorption events from hepatic 
metabolism) then bioavailability data is collected from an in vitro assay that best 
approximates the luminal barrier and segmental physiology of the gastrointestinal tract. 

Examples of some common cell and tissue sources for permeability and transport 
mechanism assays for a selected route of administration are provided below in Table 1. 
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Tabid: Permeability and Transport^^ 

Cell Culture 



Route/Tissue 

bral/lntestinai 



Inhalation/Nasal 



Ocular/Comeal 



Oral-Buccal/Cheek 
TopicalAlransdennal 



iV/Hepatic 



IV/Blood Brain Earner 



"Caco-2 cells 
HT-29 cells 
T84 cells 

Intestinal epithelial cells (lEC) 
SV40 T Immortalized cells 
Organ culture/co-culture 
Primary culture 



SV40 T immortalized cells 
Primary culture 



RCEl cells 

Primary cultures 

SV40 T immortalized cells 



Primary cultures 

HaCat celir 
primary/co-cultures 



Hepatic carcinoma cell lines 
Primary cultures 
Co-cultures 

SV40 T immortalized cells 

Primary culture 
SV40 immortalized cells 



Bxan^ples of so.e conunon parameters for solubility and dissolution assays for 
given route of administration are provided below in Table 2. 
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Table 2: Solubilitj- and Dissolution Parameters. 

In vitro Parameters 



Oral 



Route/Anatomy/Physiology 

Gastrointestinal (GI) 



Buccal/Sublingual 



Rectal 



tract 



Stomach 

Duodenum 

Jejunum 

Ileum 

Colon 



Parenteral 



Mouth 
• Cheek 
Tongue 

Lower Gl tract 
Colon 
Rectum 



Skin 
Muscles 
Veins 



pH 

Temperature 

Concentration of test sample 
Volume 

Osmotic pressure 
Admixing conditions 
Physiologic 

Fluid/Buffer/solvent system 

Excipients 

Other Additives 

Test chamber composition 




Respiratory system 

Nose 
• Lungs 
Mouth 



■ . .hninues for collecting permeability and transport mechanism 
;„ and in v.vo ^^^^^^ ,^ art (Stewart et 

data using cell- and/or ttssue-based P'«P^^°" ^ 7„ ^35^51. Minth et al.. 

» n oos^ 12-693-699; Andus et al., Pharm. Res. 
al.. Pharm. Res. (1995) 12.693 . unv461-473). For instance, in vitro 
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wooo/1623' ^eaWlitv and transport mechamsm data typically are 

a..a,s. '"viv„-V=^-°"-"""°™:r; ra.«. hamster, .o. and monkey to 
oomiuoted in anitnal models streh as . ^^^^^ i,^bmo^ metaboKsm, 

• .Haractertze bioavaiUWl.-y -^^^ „„ ,„,„e-based . vi,™ assays are 
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oavWe. of Ae h«r. and of .he blood and lymph vessels, and U,= serious cavihes of the body, 
originating ftom the mesodenn. 

One of ordinary skill in me art will recognize that cells and tissues can be obtained de 
novo from a sample of intereat, or from cxisdng sources. Public sources include cell «ad cell 
5 line repositories such as U,e American Type Cul»e Collectton (ATCC), U.e Belgian Cul.« 
CoUections of Microorgamsms (BCCM), or the Qerrn^ CoUection of Microorganisms and 
ecu cultures (DSM), among many others. The cells can be cultivated by standard techmques 
known in the art. 

Preferred assays for collecting permeability data utilize devices and methods that 
,0 meagre change in resistance or conductivity of a membrane system by ion flux. Any devtcc 
stable for such studies can be employed. These include voltage-clamp type devtces and 
methods that employ either cell cultures or precision tis^e slices. Diffttsion chamber 
sy««ns utilizing culfl«d cells grown on permeable supports to measure permeabthty are 
preferred More preferred devices are readily adapted for high-throughput and automated 
,5 screemng. Examples of such devices are known and exemplified in U.S. Patent No 
5 599 688- WO 96/13721; ^ WO 97/16717, These devices also can be adapted for 
transport mechanisms.. As can be appreciated, however, measur«nen. of 
reslmnc conductivity and/or ion flux is no. required .o detemune pemreabtUty of 
compounds. Many otiter «chniques are a«aiable and can be employed in the invention. For 
20 instance, permeability d«a also may be predicted using Uteoretical models to approxm.a.e 
ti.is parameKr. for example, fiom SAR/QSAR (e.g.. log P, molecular weight, H-bondmg, 
surface properties). 

Transport mechanism of a test sample of mteres. can be detemtined using ceU culti^es 
and/or tissue s«=tions ibllowing standanj »chni,«es. These assays typicaUy involve 
25 conti^ting cells or tissue witi. a compound of int«est and measuring uptake into ti« cells, or 
competing for uptake, compared to a known transport-specific substrate. These «,perimen.s 
can be pe,fonn«l a. short incubation times, so that ktaettc paramet«s can be measured that 
will accurately chanKterize ti>e transporter systems, and mimmize a,e effects of non- 
aat^ting passive functions. (Baley et al.. Advanced Dn.g Delivery Reviews (1996) 22.85- 
103)- Hidalgo e. .1., Advanced Drtig Delivery Reviews (1996) 22:53-66; Andus e. al., Pharm. 
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Res. (1990) 7(5):435-451). For high-throughput analyses, cell suspensions can be employed 
utilizing an automated method that measures gain or loss of radioactivity or fluorescence and 
the like such as described in WO 97/49987. 

In a preferred embodiment, transport mechanism is determined using high-throughout 
transporter screening cell lines and assays. In this aspect of the invention a cell line is 
selected and/or manipulated to over-express one or more transporter proteins, and/or 
enzymes. The cells are then used to rapidly identify the mechanism(s) by which a compound 
is transported across the physiological barrier of interest. Transporters of interest represent 
the basic categories of transport including uptake and efflux transporters. These transporters 
aid in the movement of materials in biological systems, into and out of cells and across 
cellular layers. Natural combination(s) of enzyme(s) and transporter(s) also can provide the 
basis of a high-throughput transport mechanism screening assay. For instance, certain 
enzymes or transporters require secondary enzymes or transporters to function in a noimal 
physiological mode, i.e., cytochrome P4503A is co-regulated with P-glycoprotein. These 
proteins share the same substrate and their genes are co-regulated. Thus multiple artificial 
combination(s) of transporter(s) and enzyme(s) can be employed for characterizing transport 
mechanism of a test sample of interest. Examples of possible combinations of a fransporter 
and/or enzyme in a host cell of interest include cell-tiransporter-enzyme, cell-transporter, cell- 
enzyme, cell-enzyme-enzyme, and cell-dransporter-transporter. Examples of transporters that 
can be used to transfect the host cell of interest include peptide transporters (PepTl), amino 
acid transporters, organic cation transporters (OCTl), organic anion transporters, nucleotide 
transporters (Nl. N2. N3, ES, EI), glucose transporters (SGLTl, GLUT 1 through GLUT 7). 
monocarboxylate transporters (MCTl), and multi-drug transporters (LRP, MDR, MRP, 
PGP). Examples of enzymes that can be used to transfect the host cell are Phase I and H 
enzymes, cytochrome P450, 3A, 2D and tiie like. 

Nucleic acid and/or amino acid sequences for transporters/enzymes can be identified 
in various genomic and protein related databases. Examples of publicly accessible databases 
include as GenBank (Benson et al., Nucleic Acids Res (1998)26(1): 1-7; USA National Center 
for Biotechnology Infoimation. National Library of Medicine, National Institutes of Health, 
Bethesda, MD, USA), TIGR Database (The Institute for Genomic Research, Rockville, MD, 
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USA), Protein Data Bank (Brookhaven National Laboratory, USA), and the ExPASy and 
Swiss-Protein database (Swiss Institute of Bioinformatics, Geneve, Switzerland). 

Any number of known techniques can be used to prepare nucleic acid encoding a 
transporter(s) and/or enzyme(s) of interest. To express a target protein in a host cell the 

5 nucleotide sequence coding for the polypeptide is inserted into an appropriate expression 
vector, i.e., a vector that contains the necessary elements for the transcription and translation 
of the inserted coding sequence. The host cell line can be stably or transiently transfected by 
methods known in the art. Examples of transient transfection methods include calcium 
phosphate, electroploration, lipofectamine, and DEAE dextran. A cell Une can be stably 

1 0 transfected using methods known in the art such as calcium phosphate. In addition, the host 
cell can be infected with a retrovirus containing a target protein of interest, resulting in stable 
expression of the desired target protein. Host cells that express the target gene product can be 
identified by standard techniques. These include, but are not limited to, detection of the 
protein as measured by immunoprecipitation and Western blot analysis or by measuring a 

15 specific biological response. 

For synthesis in a ceU, a target transporter/enzyme protein can be generated by 
standard techniques. Cells that naturally express a target protein can be employed. 
Transfection and transformation of a host cell with DNA encoding a protein of interest also 
can be used. For example, a polymerase chain reaction (PGR) based strategy may be used to 

20 clone a target DNA sequence encoding all or part of a target membrane polypeptide of 
interest. (See, e.g., "PGR Cloning Protocols: From Molecular Cloning to Genetic 
Engineering," B.A. White, ed., Humana Press, Methods in Molecular Biology, Vol. 67, 
1997). For example, PGR can be used for cloning through differential and subtractive 
approaches to cDNA analysis, performing and optimizing long-distance PGR, cloning 

25 unknown neighboring DNA, and using PGR to create and screen libraries. PGR also can be 
used to introduce site-specific and random mutations into DNA encoding a target protein of 
interest. 

For general cloning purposes, complementary and/or degenerate oligonucleotides 
corresponding to conserved motifs of the target membrane polypeptide may be designed to 



35 



PCT/US99/21151 

WO 00/16231 

serve as primers in a cDNA and/or PGR reaction. Templates for primer design can be 
obtained from any number of sources. For example, sequences, including expressed 
sequence tags (ESTs) can be obtained from various databases, such as GenBank, TIGR, 
ExPASy and Swiss-Protein databanks. Homology comparisons performed using any one of a 
number of alignment readily available programs that employ search engines to find the best 
primers in a sequence based on various algorithms. Any number of commercially available 
sequence analysis packages, such as Lasergene, GeneWorks, DNASIS, Gene Jockey II. Gene 
Construction Kit, MacPlasmap. Plasmid ARTIST, Protein Predictor, DNA/RNA Builder, and 
Quanta. (See, e.g., "Sequence Data Analysis Guidebook," Simon R. Swindell, ed., Humana 
Press. 1996). The information can be used to design degenerate primers, nested/multiplex 
primers, site-directed mutagenesis, restriction enzyme sites etc. Primers can be designed 
from homology information, and computer programs can be used for primer design as well. 
Examples include "Primer Premier 4.0" for automatic primer selection (CloneTech. Inc.). 
The amplified cDNA and/or PGR firagment may be used to isolate fiilMength clones by 
radioactive or non-radioactive labeling of the amplified fragment and screening a library. 

Alternatively, transporter/enzyme DNA cloned from one source may be utilized to 
obtain a corresponding DNA sequence from other sources. Specifically, a genomic and/or 
cDNA library constructed from DNA and/or RNA prepared from a cell known or expected to 
express the target transporter/enzyme may be used to transform a eukaiyotic or prokaryotic 
host cell that is deficient in the putative gene. Transformation of a recombinant plasmid 
coding for the protein into a deficient host cell would be expected to provide the cell with a 
complement product corresponding to the protein of interest. In some cases, a host cell can 
be selected to express a particular phenotype associated with the target polypeptide and thus 
may be selected by this property. For a review of cloning strategies which may be used, see 
e.g., Sambrook et al., 1989. Molecular Cloning, A Laboratory Manual, Cold Springs Harbor 
Press, New York; and Ausubel et al.. 1989. Current Protocols in Molecular Biology. Green 
PubUshing Associates and WUey Interscience. New York. 

To express a target transporter/enzyme in a host cell the nucleotide sequence coding 
for the protein, or a fimctional equivalent for modular assembly as described above, is 
inserted into an appropriate expression vector, i.e.. a vector which contains the necessary 
elements for the transcription and translation of the inserted coding sequence. Host cells 
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containing the coding sequence and that express the target gene product may be identified by 
standard techniques. For example, these include but are not limited to DNA-DNA or DNA- 
RNA hybridization; the presence or absence of "marker" gene functions; assessing the level 
of transcription as measured by the expression of mRNA transcripts in the host cell; and 
5 detection of the gene product as measured by immunoassay or by its biological activity. 

Once a clone producing the target transporter/enzyme is identified, the clone may be 
expanded and used to over express the protein(s). If desired, the proteins may be purified 
using techniques well-known in the art including, but not limited to immunoaffinity 
purification, chromatographic methods including high performance liquid chromatography or 
10 cation exchange chromatography, affinity chromatography based on affinity of the 
polypeptide for a particular ligand, immunoaffinity purification using antibodies and the like. 
The purified proteins can then be bound to an artificial membrane matrix and utilized for 
assessing interaction of compounds to the transporter/enzyme of interest. 

Some commonly used host cell systems for expression of transport proteins and 
15 enzymes include E. coli, Xenopus oocytes, baculovirus, vaccinia, and yeast, as well as many 
higher eukaryotes including transgenic cells in culture and in whole animals and plants. (Sec. 
e.g., G.W. Gould, "Membrane Protein Expression Systems: A User's Guide," Portland Press, 
1994. Rocky S. Tuan, ed.; and "Recombinant Gene Expression Protocols," Humana Press, 
1996). For example, yeast expression systems are well known and can be used to express and 
20 recover target transporter/enzyme systems of interest following standard protocols. (See. e.g., 
Nekrasova et al, Eur. J. Biochem. (1996) 235:28-37; Gene Expression Technology Methods 
in Enzymology 185: (1990); Molecular Biology and Genetic Engineering of Yeasts CRC 
Press, Inc. (1992); Herescovics et al.. FASEB (1993) 7:540-550; Larriba, G. Yeast (1993) 
9:441-463; Buckholz, R.G., Curr Opinion Biotech (1993) 4:538-542; Mackett. M, 
25 "Expression of Membrane Proteins in Yeast Membrane Protein Expression Systems: A Users 
Guide," pp. 177-218, Portland Press. (1995). 

For high-resolution screening and validation, tissue-based assays may be employed to 
characterize transport mechanisms. For example, of the cytochrome P450 superfamily, 
CYP3A enzymes represent the most abundant isoforms in the liver and they are responsible 
30 for the metabolism of compounds of diverse chemical structure. The uptake of a compound 
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into hepatocytes can be mediated by passive or carrier processes. Once in the parenchymal 
cell of the liver, the drug can be metabohzed or bind to intracellular proteins. The drug or its 
metabolite(s) may return to the circulation or exit from the hepatocyte into the bile 
canaliculus, again by passive or carrier-mediated transport, before secretion in bile. 
Experimental systems have been devised to study these processes in isolation. Examples of 
such systems include isolated perfused rat liver (IPRL), and bile duct camiulated (BDC) rat 
models. (Chan et al., DDT (1996) 1:461-473). 

Tissue from transgenic animals designed to express particular transport properties in 
one or more particular tissues also may be utilized to characterize transport mechanisms. In 
this aspect of the invention, an animal can be genetically manipulated to express or not 
express one or more specific proteins in a tissue of interest, e.g. transporter protein in 
duodenum tissue. Tissue from the genetically engineered animal can then be used to examine 
transport mechanisms in a tissue-based assay. Transgenic, animal methodologies are well 
known (Gordon et al.. Hum. Cell (1993) 6(3):161-169; and Jaenisch, R., Science (1998) 
240:1468-1474). 

Artificially engineered tissue also can be used for permeability assays, such as tissues 
generated ex vivo for use as skin grafts, transplants, and the like. Such tissues can be 
obtained using standard techniques. See, for example, U.S. Patent Nos. 5.759,830; 
5,770,193; and 5,770.417. 

Solubility and dissolution data can be obtained in an in vitro assay by testing each 
sample of interest in an appropriate physiologic fluid^uffer system that best approximates 
the particular physiological system selected as the barrier to absoiption. A solubility profile is 
a plot of solubility of a test sample at various physiological conditions. As an example, the 
natural pH environment of the gastrointestinal tract varies from acidic in the stomach to 
sUghtly alkaline in the small intestine and fluid composition for each segment may vary as 
well. The solubility profile provides an estimation of the completeness of dissolution of a test 
sample in a particular physiological compartment or anatomical entity. In this instance, a 
panel of test wells each having different pHs and physiological fluid composition can be 
employed to generate a solubility profile for each test sample. Solubility and dissolution data 
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can also be predicted using theoretical models to approximate these values, for example, from 
SAR/QSAR information. 

In vitro dissolution assays measure the rate and extent of dissolution of a test sample 
in an aqueous solution. Various parameters are considered when performing a dissolution 
assay and are well known in the art. These parameters include size of the experimental 
vessel, amount of agitation and nature of the stirrer, temperature and nature of the dissolution 
medium. pH, viscosity, and design of the dissolution apparatus. Standard methods known in 
the art for measuring dissolution include rotating basket, paddle, rotating bottle, flow-through 
dissolution, intrinsic dissolution, and peristalsis methods. These methods can be adapted and 
used as a guide for high-throughput solubility and dissolution screening. 

For high-throughput collection of solubility and dissolution data, automated methods 
of solid and liquid handling are employed. This method involves addition of samples to a 
multi-well or multi-tube/plate system. The data associated with these tubes/plates, such as 
physiologic fluidA)uffer system, volume, concentration, pH and tube/plate maps, is 
transferred into an inventory system. The inventory system generates codes containing 
updated information pertaining to the aliquoting. diluting, or pooling methods applied to the 
original tubes/plates. Tasks created.in the database are then carried out physically in coded 
tubes/plates. Aliquots are then distributed to designated screen sites. After testing, the 
solubility profiles are generated and ported to a database for access and analysis. 



Assays for Screening Secondary Absorption libraries: 

Secondary libraries selected for absorption also can be characterized by one or more 
additional properties including, but not limited to. metabolism, distribution, elimination, 
toxicity, and biological activity. As with absorption, assays to characterize the relevant data 
are based on the selected route of administration. MetaboUsm or biotiansformation refers to 
the biochemical transformation of a compound to another chemical form. The 
biotransformation process typically results in a metabolite that is more polar (water-soluble) 
than the original parent molecule. 
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Most tissues have some metabolizing capacity but the liver is by far the most 
important organ, on the basis of size if not always concentration of target compound 
metabolizing enzyme. Phase I reactions are defined as those that introduce a fimctional 
group to the molecule and phase 11 reactions are those that conjugate those function groups 
with endogenous moieties. 

Since metabolism is a drug clearance process, metabolism of a compound contributes 
to elimination of the compound. Thus, compounds selected for absorption can be screened 
for metabolism in order to consider disposition of a drug after or concurrent with 
administration using standard techniques known in the art. (See, e.g., Sakuma & Kamataki, 
Drug metabolism research in the development of innovative drugs. In: Drug News & 
Perspectives (1994) 7 (2):82-86). 

Metabolism assays for high-throughput screening preferably are cell-based (cells and 
cellular preparations), whereas high resolution screening can employ both cell and tissue- 
based assays. In particular, test samples from compound libraries can be screened in cell and 
tissue preparations derived from various species and organs. Although liver is the most 
frequently used source of cells and tissue, other human and non-human organs, including 
kidney, skin, intestines, lung, and blood, are available and can be used to assess extra-hepatic 
metabolism. Examples of cell and tissue preparations include subceUular fractions (e.g., liver 
S9 and microsomes), hepatocytes (e.g.. coUagenase perfiision, suspended, cultured), renal 
proximal tubules and papillary ceUs, re-aggregate brain cells, bone marrow cell cultures, 
blood cells, cardiomyocytes, and established cell lines as well as precision-cut tissue slices. 

Examples of in vitro metabolism assays suitable for high-throughput screening 
include assays characterized by cytochrome P450 form-specific metabolism. These involve 
assaying a test compound by P450 induction and/or competition studies with form-specific 
competing substrates (e.g., P450 inhibitors), such as P450 enzymes CYPIA, 3A, 2A6, 2C9, 
2C19, 2D6, and 2E1. Cells expressing single or combinations of these or other metabolizing 
enzymes also may be used alone or in combination with cell-based permeability assays. A 
high-throughput cell-based metabolism assay can include cytochrome P450 induction 
screens, other metabolism marker enzymes and the like, such as with measurement of DNA 
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or protein levels. Suitable cells for metabolism assays include hepatocytes in primary 
culture. Computer-implemented systems for predictmg metabolism also may be employed. 

Absorption libraries also may be characterized by additional distribution and 
elimination events. In this aspect of the invention, in vitro assays are performed to assess 
protein binding to a test compound, since protein binding can affect compound distribution 
and elimination. In general, it is free compound that diffuses into cells and tissues. Binding 
can be classified as restrictive or permissive with regard to elimination, or quantitatively 
defined in terms of affinity. Affinity of the binding is defined as low or high when reversible, 
or more unusually when irreversible binding occurs. The biological half-life of a test 
compound will increase due to its interaction with a protem. Usually, the higher the affinity 
the lower the elimination that may be observed. Albumin is by far the most firequent 
contributors to plasma protein binding since it comprises about one half of the total plasma 
proteins. The al-Acid glycoprotein also plays an important role in the protein binding of a 
compound since it has an affinity for bases (many drugs are weak bases). It is an acute phase 
reactant and its concentration rises in inflammatory processes, malignant disease and stress. 
Lipoproteins (HDL, LDL or VLDL) bind drugs that are highly liposoluble and a fairly 
specific ligand-protein interaction occurs between certain steroids and gamma globulins. 
Thus, in vitro protein binding assays that employ one or more of albumin, al-acid 
glyc<!protein, lipoprotein, steroid and gamma globulins may be utilized to collect distribution 
and elimination data that can be utilized to further characterize an absorption library. 

Similarly, toxicity of a test compound may also be assayed and used to characterize 
compounds of an absorption library. Any number of techniques in the art may be employed 
for this purpose. Preferred methods are in vitro. Examples include determination of toxicity 
mechanisms, determination of cytotoxic potentials in cell and tissues of target organs, 
estimation of therapeutic indices from in vitro data, cytotoxicity screening of closely related 
drug compounds in cells from the same mammal or from different species, detection and 
quantification of peroxisome proliferation, screening of agents to prevent or reverse 
cytotoxicity, and specialized studies on target cells using co-incubation systems, e.g., red 
blood cells and hepatocytes. 
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Toxicity assays may utilize any technique that provides a toxicity parameter as an 
endpoint. For high-throughput screening, cell based assays are preferred. This includes gene 
expression (e.g., protein or nucleic acid based) enzymatic activity, and morphology screens 
and the like. Examples of cell-based assays include in vitro peroxisome proliferation studies, 
which can be used to assay palmitoyl CoA-oxidation in primary hepatocyte culture, with or 
without concurrent measurement of DNA or protein levels. Cytotoxicity assays in primary 
cultures also can be utilized, and include screening for cytotoxicity in hepatocytes or renal 
proximal tubules, enzyme release (lactate dehydrogenase), and MTT conversion 
(mitochondrial function) following standard techniques. Computer-implemented SAR/QSAR 
models for predicting toxicity also may be employed, such as when structural information is 
available. 

Absorption libraries produced according to the method of the invention also may be 
examined for activity hits using any technique suitable for such purpose. Examples include 
screening of isolated receptors or use of cellular preparations that contain a receptor target of 
interest (i.e.. compound/ligand-receptor interaction^inding). These include reporter gene 
assays, binding assays, cellular proliferation assays and the like. (See, e.g.. Wallace, R.W., 
and Goldman, M.E.. Bioassay Design and Implementation, High Throughput Screening, 
(1997) p. 279-328, Ed. Devlin, P). Activity assays also may use SAR/QSAR models. 

Screening of secondary absorption libraries by one or more of such additional 
properties can be performed concurrently or following the initial absorption screen of a 
primary compound library. 

PKTool and System: 

The PK tool of the invention is utilized to generate a simulated in vivo absorption 
profile from in vitro solubility and permeabiUty data, and optionally in vitro dissolution rate 
and transport mechanism data for a test sample of a compound library. The PK tool includes 
as computer-readable components, an input/output system suitable for data input and data 
output, a simulation engine having a numerical-based differential equation solver, and a 
physiologic-based simulation model comprising a pharmacokinetic model of the mammahan 
system to be simulated. In vitro bioavailability data is provided through the input/output 
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system, and then the simulation engine and simulation model are applied to facilitate a 
simulation run so as to generate a user selected in vivo absorption profile for the test sample. 
Together, the simulation engine and simulation model are employed to simulate the fate of a 
test sample in the system under investigation. 

The PK tool is based on a compartment-fiow simulation model system. The 
compartment-flow model employs compartments, flow regulators, and converters that 
collectively regulate flow among the compartments. The model components are represented 
by a series of differential equations which when run through the simulation engine are solved 
at each time increment dt based on the initial underlying values of the equations, the input 
values supphed by the user, and calculations performed by various subsystems of the model 
when activated in a particular scaiario. 

The PK tool optionally comprises a repository of different pharmacokmetic models 
and initial parameter values for a given model. The repository preferably resides in a 
database of the PK tool, and/or is accessible through an acquisition model. The PK tool also 
may include one or more curve-fitting algorithms for generation of absorption parameters and 
constants for correlation of in vitro data to in vivo data, or in vivo data from one species of a 
mammal to in vivo data of a second species of mammal based on a selected route of 
administration. The curve-fitting algorithms include regression-based and stochastic-based 
algorithms. 

1. Input/Output System 

With regard to the components of the PK tool, the input/output system provides a user 
interface between the user and other components of the PK tool of the invention. The 
input/output system may be any suitable interface between user and computer system, for 
input and output of data and other information, and for operable interaction with a simulation 
engine and a simulation model. For instance, the input/output system may provide direct 
input from measuring equipment. The input/output system preferably provides an interface 
for a standalone computer or integrated multi-component computer system having a data 
processor, a memory, and a display. Input mto the method and PK tool of the invention is in 
vitro bioavailabiUty data derived from an assay corresponding to a selected route of 
administration and mammalian system of interest. For example, the user enters the initial 
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parameter values for a test sample, such as dose, permeabiUty, solubility, and the like, and 
then optionally indicates the transport mechanism, e.g., passive transcellular, passive 
paracellular, carrier-mediated influx, or carrier-mediated efflux. When transport mechanism 
is not indicated, the PK tool can be designed to employ a default transport mechanism, such 
as passive transcellular. When set to the paracellular mechanism, the absorption of the 
compound is adjusted to compensate for the lower surface area available for absorption via 
the paracellular pathway. The model also may incorporate an operation by which the 
mechanism of absorption can be predicted using the permeability, solubility, molecular 
stmcture or other information. This allows the model to automatically compensate for 
paracellular and/or other absorption mechanisms without requiring prior input and knowledge 
from the user. Depending on the objective, the user also may specify the pH, delivery system 
rate such as controlled release rate or formulation release rate (delivery system referred to 
herein as "formulation"), dosing schedule, and simulation run time, as well as physiologic 
system specific parameters such as GI transit time when a'GI tract model is employed. If 
values for these parameters are not entered, the PK tool provides default values. 

Data may be entered numerically, as a mathematical expression or as a graph that 
represents a physiological or pharmacokinetic parameter, or alpha such as transcellular, 
paracellular, passive, active, etc. An advantage of entering data as a graph is that it removes 
any requirement to define the mathematical relationship between a dependent and an 
independent variable. The interface output displays and/or compares parameters related to 
absorption, such as graphs or tables corresponding to rate of absorption, extent of absorption, 
and concentration profiles, and the like. Output of the method and PK tool is utilized to 
profile and rank the test sample by one or more selected absorption parameters. 

The absorption parameters of a profile include concentration, rate and/or extent of 
absorption of a test sample. As can be appreciated, absorption parameters can be represented 
in multiple different ways that relate time, mass, volume, concentration variables, firaction of 
the dose absorbed and the like. Examples include rate "dD/dt" and "dc/dt" (e.g.. mass/time- 
mg/hr. concentration/time-jig/ml/hi), concenti-ation "C (e.g., mass/volume-jig/ml), area 
under the curve "AUG" (e.g., concentration • time, ^^g • hr/ml), and extent/fraction of the 
dose absorbed 'T" (e.g., no units, 0 to 1). Other examples include tiie maximum 
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concentration (C„«x), which is the maximum concentration reached during the residence of a 
compound at a selected sampling site; time to maximum concentration (T™,^), which is the 
time after administration when the maximum concentration is reached; and half-life Hm), 
which is the time where the concentration reaches V2 its maximum at a selected sampling site. 
Other examples of output include individual simulated parameters such as permeability, 
solubility, dissolution, and the like for individual segments, as well as cumulative values for 
these and/or other parameters. 

2. Simulation Engine 

The simulation engine comprises a differential equation solver. The simulation 
engine also may include a system control statement module when control statement rules 
such as IF...THEN type production rules are employed. The differential equation solver uses 
standard numerical methods to solve the system of equations that comprise a given 
simulation model. These include algorithms such as Euler's and Runge-Kutta methods. 
Such simulation algorithms and simulation approaches are well known (See, e.g., Acton, F.S., 
Numerical Methods that Work, New York, Harper & Row (1970); Burden et al.. Numerical 
Analysis, Boston, MA, Prindle, Weber & Schmidt (1981); Gerald et al.. Applied Numerical 
Analysis, Reading, MA, AddisomWesley Publishing Co.. (1984); McCormick et al., 
Numerical Methods in Fortran. Englewood Cliffs. NJ. Prentice Hall. (1964); and Benku, T., 
The Runge-Kutta Methods. BYTE Magazine, April 1986. pp. 191-210). 

Many different numerical schemes exist for the evaluation of the differential 
equations. There are literally lOO's of schemes that currently exist, including those 
incorporated into public commercially available computer applications, private industrial 
computer applications, private individually owned and written computer applications, manual 
hand-calculated procedures, and published procedures. With the use of computers as tools to 
evaluate the differential equations, new schemes are developed annually. The majority of the 
numerical schemes are incorporated into computer appUcations to allow quick evaluation of 
the differential equations. 

Computer application or programs described as simulation engines or differential 
equation solver programs can be either interpretive or compiled. A compiled program is one 
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that has been converted and written in computer language (such as C-f+, or the like) and are 
comprehendible only to computers. The components of an interpretive program are %vritten in 
characters and a language that can be read and understood by people. Both types of programs 
require a numerical scheme to evaluate the differential equations of the model. Speed and run 
5 time are the main advantages of using a compiled rather than a interpretive program. 

A preferred simulation engine permits concurrent model building and simulation. An 
example is the program STELLA® (High Performance Systems, Inc.). STELLA® is an 
interpretive program that can use three different numerical schemes to evaluate the differential 
equations: Euler's method, Runge-Kutta 2. or Runge-Kutta 4. The Kinetica™ program 
10 (InnaPhase, Inc.) is another differential equation solving program that can evaluate the 
equations of the model. By translating the model ftom a STELLA® readable format to a 
Kinetica™ readable format, physiological simulations can be constructed using Kinetica™. 
S which has various fitting algorithms. This procedure can be utilized when the adjustment 

parameters are being optimized in a stepwise fashion. 

15 3. Simulation Model 

The simulation model is a mathematical model of a multi-compaitment physiological 
model of a mammalian system (e.g., GI tract) that corresponds to the selected route of 
administration (e.g.. oral). A given physiological model is represented by series of 
differential equations that describe rate process interactions among anatomical segments for 
20 the physiological system under investigation. The individual segments or compartments are 
represented mathematically as a one, two and/or three compartment kinetic system. The 
segments are linked in a stepwise fashion so as to form an integrated physiological model 
describing absorption of a compound relative to the anatomical segments and at least one 
sampUng site for assessing an absorption event in isolation. For a model simulating oral 
25 delivery, anatomical segments of the GI tract are provided, which can include tiie stomach, 
duodenum, jejunmn. ileum and colon. A sampling site for the GI tract may be the portal vein 
and/or plasma. The rectum and colon would be applicable for modeling a rectal route of 
delivery. Segments and sampling site for buccal or sublingual delivery routes can include tiie 
mouth, cheek/tongue tissue and plasma. For ocular routes, this can include aqueous humor. 
30 conjunctival sac, tear duct, nasal cavity and plasma. For tiie lung routes, tiiis can include 
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respiratory bronchioles zone and plasma. For delivery via the nose, this can include nasal 
cavity and plasma. For the topical and transdermal routes, this can include epidermal, 
" dennal, subcutaneous tissue, muscle and plasma. Other systems adhere to these bas,c 

designs. 

5 Of course compartments representing a particular anatomical segment can be added or 

removed depending on the model's intended end use, such as when an isolated segment is 
examined or when it is desired to account for parameters affecting bioavailabibty at 
additional sampling sites. For example, compartments can be added to account for both pre- 
or post-absorptive protein binding or complex fom^ation to account for reversible association 

10 of a compound to the proteins (albumin and al-acid glycoprotein) of blood, or more usually 
plasma. Other compartments that may be added would include those that account for phase I 
and/or phase II hepatic metabolism. Formulation compartments that account for variable 
compound fomiulations also can be added, such as time-release, extended release or 
otherwise controlled release formulations. Another example is inclusion of kidney 

1 5 compartments to account for renal clearance. 

The compartments can be modified by factors that influence absorption such as mass, 
volume, surface area, concentration, permeability, solubility, fluid secretion/absorption, fluid 
transit, mass transit and the like, depending on the physiological system under investigation. 
As a rule of thumb, compartment modifiers relate to input variables. For instance, where 
20 transport mechanism and dissolution rate are variables considered for generating an 
absorption profile, then the physiological model will include compartments and parameters 
that account for these variables. 

When represented as a compartment-flow simulation model, the anatomical segments 
of a physiological model typically include one or more central and peripheral compartments 
25 that reversibly communicate through a flow regulator. A central compartment represents the 
interior or mucosal side of an anatomical segment. A peripheral compartment represents the 
blood side of the segment. The central and peripheral compartments are comiected by a flow 
regulator representing a physiological barrier through which material from the central 
compartment "flows" or is transferred to the peripheral compartment at some empmcally 
30 defined or calculated transfer rate "kl2" applied by a converter, which allows calculation of 
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parameters using compartment values. Transfers ("flows") between compartments can be 
zero order, second order and/or mixed order processes. As an example, a first order transfer 

■ from central compartment 1 to peripheral compartment 2 can be defmed by a finite difference 
equation connecting inputs (e.g., rate constant kl2 and amount of compound in central 

5 compartment = amount + dt*(-elimination - kl2+ k21)) to the flow controller between the 
compartments (e.g., kl2) and setting it as the product of the two variables. Thus the 
underlying equations of the model are utilized to calculate the amount of a compound in each 
compartment, and standard differential equations interrelate the system of compartments and 
their equations. This permits the model to simulate movement of a compound through each 

10 ■ compartment according to the calculated rates at each time increment (dt). Since all 
movement between compartments is in units of mass, the blood side and transferred test 
compound concentration is calculated from the amount of compound in the blood side 
(peripheral compartment) and volume of the mucosal side (central compartment). A model 
cycle is entered through the input/output user interface as incremental pulses (to simulate 

15 ramp, plug flow/lag tiroes) or as a fixed time range to initiate and effectuate cycling of a test 
compound of interest. 

The basic structure of a physiological model and mathematical representation of its 
interrelated anatomical segments can be constructed using any number of techniques. The 
preferred techniques employ graphical-oriented compartment-flow model development 
20 computer programs such as STELLA®, KINETICA™ and the like. Many such programs are 
available, and most employ graphical user interfaces for model building and manipulation. In 
essence, symbols used by the programs for elements of the model are arranged by the user to 
assemble a diagram of the system or process to be modeled. Each factor in the model may be 
programmed as a numerical constant, a linear or non-linear relationship between two 
25 parameters or as a logic statement. The model development program then generates the 
differential equations corresponding to the user constmcted model. For example, STELLA® 
employs five basic graphic tools that are linked to create the basic structure of a model: (1) 
stocks; (2) flows; (3) converters; (4) input links; and (5) infinite stocks (See, e.g., Peterson et 
al., STELLA® H, Technical Documentation. High Performance Systems, Inc., (1993)). 
30 Stock are boxes that represent a reservoir or compartment. Flows or flow regulators control 
variables capable of altering the state of compartment variables, and can be both uni- and bi- 
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directional in terms of flow regulation. Thus, the flow/flow regulators regulate movement 
into and out of compartments. Converters modify flow regulators or other converters. 
Converters ftmction to hold or calculate parameter variable values that can be used as 
constants or variables which describe equations, inputs and/or outputs. Converters allow 
calculation of parameters using compartment values. Input links serve as the internal 
communication or connective "wiring" for the model. The input links direct action between 
compartments, flow regulators, and converters. In calculus parlance, flows represent time 
derivatives; stocks are the integrals (or accumulations) of flows over time; and converters 
contain the micro-logic of flows. The stocks are represented as finite difference equations 
having the following form: Stock(t) = Stock(t-dt) + (Flow)*dt. Rewriting this equation with 
timescripts and substituting t for dt: Stock, = Stock..A, + At*(Flow). Re-arranging terms: 
(Stock, - Stockt-AtVAt = Flow, where "Flow" is the change in the variable "Stock" over the 
time interval "t." In the limit as At goes to zero, the difference equation becomes the 
differential equation: d(Stock)/dt = Flow. Expressing this in integral notation: Stock = 1 Flow 
dt. For higher-order equations, the higher-order differentials are expressed as a series of first- 
order equations. Thus, computer programs such as STELLA® can be utilized to generate 
physiologic-based multi-compartment models as compartment-flow models using graphical 
tools and supplying the relevant differential equations of pharmacokinetics for the given 
physiologic system under investigation. An example of iconic tools and description, as well 
as graphically depicted compartment-flow models generated using STELLA® and their 
relation to a conventional pharmacokinetic IV model is illustrated in Figure 6-9. 

The model components may include variable descriptors. Variable descriptors for 
STELLA®, for example, include a broad assortment of mathematical, statistical, and built in 
logic functions such as boolean and time functions, as well as user-defined constants or 
graphical relationships. This includes control statements, e.g., AND, OR, IF ... THEN ... 
ELSE, delay and pulsing, that allow for development of a set of production rules that the 
program uses to control the model. Variable descriptors are inserted into the "converters" 
and connected using "input links." This makes it is possible to develop complex nile sets to 
control flow through the model. The amount of time required to complete one model cycle is 
accomplished by inputting a total run time and a time increment (dt). The STELLA® 
program then calculates the value of every parameter in the model at each successive time 



49 



PCTAJS99/21151 

WO 00/16231 

increment using Runge-Kutta or Euler's simulation techniques. The preferred simulation 
technique is Runge-Kutta. Once a model is built, it can be modified and fiirther refined, or 
adapted or reconstructed by other methods, including manually, by compilmg, or translated to 
other computer languages and the like depending on its intended end use. 

5 A preferred method of the invention for constructing a physiological model is 

depicted in Figure 10. This method employs a two-pronged approach that utilizes a training 
set of standards and test compounds having a wide range of dosing requirements and a wide 
range of permeability, solubility, transport mechanisms and dissolution rates to refine the rate 
process relations and generate the initial values for the underlying equations of the model. 
□ 10 The first prong employs the training/validation set of compounds to generate in vivo 
S pharmacokinetic data (e.g.. human plasma profiles). The second prong utilizes the 

m training/validation set of compounds to generate in vitro permeability, solubility, transport 

Fi mechanism and dissolution rate data that is employed to_ perform a simulation with the 

developmental physiological model. The in vivo pharmacokinetic data is then compared to 
the simulated in vivo data to determine how well a developmental model can predict the 
actual in vivo values from in vitro data. The developmental model is adjusted until it is 
capable of predicting in vivo absorption for the training set from in vitro data input. Then the 
model can then be validated using the same basic approach and to assess model performance. 

In particular, three primary sets of data are generated from the training set for the 
20 comparison. The first set of data is empirically derived in vivo plasma data from animals or 
humans. The second set of data is obtained from conversion of the in vivo plasma data to a 
form corresponding to the primary sampling site of the developmental physiological model. 
The third set of data is empirically derived in vitro bioavailability data including 
permeability, solubiUty, dissolution rate and transport mechanism data. The raw data points 
25 are preferably collected and statistically analyzed to provide the best fit data. The best fit 
data may be obtained by any number of curve-fitting approaches, including standard 
regression techniques. 

The in vivo plasma data is utilized to judge how well a developmental simulation 
model is able to predict absorption of the training set of compounds relative to the 
30 empirically derived in vivo plasma values. Plasma data also is utilized to calculate absorption 
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at the relevant primary sampling site of the developmental physiological model. For 
instance, in order to use in vivo plasma data in a developmental physiological model, the 
plasma data must first be converted to data corresponding to the primary sampling site of the 
model. If plasma is the primary sampling site then no conversion is needed. However, if 
5 plasma is not the primary sampling site, then a pharmacokinetic training/validation model 
relating the primary sampUng site and the in vivo plasma data is utilized. For example, when 
the developmental model is of the gastrointestinal tract, the portal vein can be selected as a 
primary sampling site and plasma selected as a secondary sampling site. Thus in this instance 
the in vivo plasma data is converted to portal vein data so that the parameters affecting 
10 secondary bioavailability events are separated imm the primary absorption event resulting 
n from passage of the test sample across the gastrointestinal lumen. This is accomplished by 

adding a plasma-portal vein conversion/validation model that relates in vivo plasma data to 
m portal vein data. This plasma-portal vein conversion/validation model can be separate or 

0^ integrated with the developmental model. In most cases, the plasma-portal vein model is 

S 15 based on a standard central-peripheral pharmacokinetic compartment approach for data 
conversion. The third set of data, the in vitro derived data, is utilized to run the 
p developmental model, and the simulated absorption profile from this data set is compared to 

5^ the in vivo derived plasma and simulated sampling site data. The developmental 

[| physiological model is modified until the simulated absorption profiles are in agreement with 

O 20 the in vivo derived plasma and simulated sampling site data. 

As the number of parameters for evaluation increase it becomes more important to 
isolate and test each component of the model building process by validation using a standard 
validation set of compounds. The validation set of compounds should contain a diverse set of 
compounds that represent a broad range of absorption profiles for which both in vitro 

25 penneabiUty, solubility, dissolution rate, and transport mechanism data, and in vivo plasma 
data is available. Statistical criteria such as sum of squares of the deviations between 
experimental data and calculated values obtained from the developmental physiological 
model are used to determine how well the model fits the data. If the developmental 
physiological model does not predict a good fit for the data, then the model is adjusted by 

30 isolating or including additional rate processes by an iterative approach. 
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Parameter values utilized in the underlying equations of a given physiological model 
may be provided in a database for ready access and manipulation by the PK tool of the 
invention. The database may include values for physiological parameters, such as rate 
constants and various other values employed in the PK tool. The rate constants correspond to 
time-dependent numerical constants describing rate processes (e.g., kl2 and k21). The 
physiological parameters include rate constants, permeability, solubility, transport 
mechanism and dissolution rate variables, and the like, as well as pH, volume, surface area, 
transit times, transit rates, and the like, that are based on the physiology of a given anatomical 
segment represented in a selected physiological model. 

The database also may include adjustment parameter values and/or regional 
correlation parameter values. To account for differences between in vitro and in vivo 
conditions, as well as differences between in vivo conditions of different types of mammals, 
adjustment parameters that modify one or more of the underlying equations of a given 
simulation model can be utilized to significantly improve predictabiUty. The adjustment 
parameters include constants or ranges of constants that are utilized to correlate in vitro input 
parameter values derived from a particular in vitro assay system (e.g., rabbit intestinal tissue, 
Caco-2 cells) to a corresponding in vivo parameter value employed in the underlying 
equations of a selected physiological model (e.g., human Gl tract). This aspect of the 
invention permits modification of existing physiologic-based pharmacokinetic models as well 
as development of new ones so as to enable their application for diverse compound data sets. 
The adjustment parameters are obtainable from iterative rounds of simulation and 
simultaneous "adjustment" of one or more empirically derived absorption parameters (e.g.. 
physiological parameters for different anatomical segments) until the in vitro data from a 
given type of assay (e.g., Caco-2 cell data) can be used in the model to accurately predict in 
vivo absorption in the system of interest (e.g., human GI). In particular, the adjustment 
parameters are obtained by a stepwise selective optimization process that employs a curve- 
fitting algorithm that estimates the change required in a value assigned to an initial absorption 
parameter of a developmental physiological model in order to change an output variable 
corresponding to the simulated rate, extent and/or concend-ation of a test sample at a selected 
site of administration for a mammalian system of interest. The curve-fitting algorithm can be 
regression- or stochastic-based. For example, linear or non-linear regression may be 
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employed for curve fitting, where non-linear regression is preferred. Stepwise optimization 
of adjustment parameters preferably utilizes a concurrent approach in which in vivo 
pharmacokinetic data and in vitro data are utilized simultaneously for fitting to the model. A 
few parameters of the developmental physiological model arc adjusted at a time until the 

5 simulated absorption profiles generated by the physiological model for each of the 
training/validation compounds provides a good fit to empirically derived in vivo data. An 
example of this approach is depicted in Figure 11 and 27. Utilization of adjustment 
parameters permits predictabiUty of diverse data sets, where predictabihty ranges from a 
regression coefficient (r^) of greater than 0.40. 0.45, 0.50. 0.55. 0.60, 0.65, 0.60, 0.65, 0.70, 

10 or 0.75 for 80% of compounds in a compound test set having a diverse range of dose 
requirements and a diverse range of permeability, solubility and transport mechanisms. The 
preferred predictability ranges from a regression coefficient (r^) of greater than 0.60. with a 
regression coefficient (r^) of greater than 0.75 being more preferred, and greater than 0.80 
being most preferred. Adjustment parameters utilized for in vivo to in vivo prediction (e.g.. 

1 5 dog to human) employs the same basic approach. 

The regional correlation parameters of the PK tool include constants or ranges of 
constants that are utilized to estimate a selected parameter value of a first segment of the 
mammalian system under investigation when that value is not supplied by the user. The 
model performs this estimation by utilizing a function/transformation algorithm (e.g., 
20 utilizing polynomial, exponential, logarithm, or any other variety of transformation 
approaches) in which (1) regional correlation parameter values, and (2) one or more values 
for the parameter that is supplied by the user for a second segment of the mammalian system, 
are utilized to estimate the value for the first segment. The regional correlation parameters 
may be empirically derived values or adjustment parameter values for various segments of 

25 the mammalian system of interest such as for permeability. A preferred regional correlation 
approach employs a polynomial-based correlation. The polynomial is based on die particular 
parameter to be estimated. The regional correlation is performed by logic fimction of tiie 
model, which when activated utilizes the function/transformation algorithm to perform the 
estimation. The regional correlation logic fimction of tiie model is activated when a value is 

30 missing for the selected parameter. The estimated value(s) are then utilized as input variables 
for the particular parameter in question. The model then proceeds by employing tiie 
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estimated value for subsequent simulation. Various regional correlation parameters can be 
used, such as permeability, solubility, dissolution rate, transport mechanism and the like. The 
preferred correlation parameters are for permeability. This permits the PK tool of the 
invention to predict absorption of a test sample from minimal input permeability values, such 
5 as when the simulation model is a GI tract simulation model and when cell-based assays are 
employed to provide permeability data corresponding to a given GI segment (e.g., Caco-2 
cells and colon). 

Since the parameter values are specific for a given physiological model (e.g., GI 
model-parameters. Ocular model-parameters, Blood-Brain-Barrier-parameters, etc.), 
10 parameter values are chosen accordingly. These values are obtainable de novo from 
experiments or from the literature. The preferred values are based on a diverse collection of 
training/validation compounds for which in vivo pharmacokinetic data is available. 

The various physiological models also may reside in a database, in part or in whole, 
and may be provided in the database with or without the initial parameter values. The 
15 database will preferably provide the differential equations of the model in a compartment- 
flow data structure that is readily portable as well as executable by the simulation engine. 

An integrated physiological model corresponding to the GI tract of a mammal 
constructed using STELLA® and the above-described methodology is illustrated in Figures 
25-26 and 30-40. An example of information provided by the database is illustrated in 
20 Appendix 4 for the gastrointestinal model depicted in Figures 25-26 and 30-40. 

A physiologic-based simulation model of the PK tool and method of the invention 
may optionally include a training/validation model. This aspect of the invention can be used 
for determining whether the model is specific and accurate with respect to compounds of 
known membrane transport mechanism (e.g., passive transcellular, passive paracellular, 
25 transporter involved for absorption and secretion) and/or with respect to known drug 
solubility/dissolution rate limitations. 

A validation model can be linked to the physiological model of the invention as 
iUustrated in Figure 12. The linked system is then run to access the specificity and accuracy 
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computed values for rate and extent of absorption. These values are then compared to 
empirically measured plasma values. If computed values fall outside of an acceptable range 
the model can be reevaluated for these compounds and adjusmients made to the model. 



5 Absorption Profile Ranking/Cataloging: 

Absorption profiles generated for test samples of a compound library are compared to 
select those samples having a desired absorption profile. The selection process essentially 
involves ranking the compounds by one or more absorption parameters such as rate, extent 
and/or concentration, or parameters derived therefrom. By "ranking" is intended assigmnent 
10 of a feature to a compound or pool of compounds that distinguishes it fi-om another in a 
hierarchical manner. 

The manner in which ranking is recorded is not a limitation of the invention. For 
example, ranking may be any recorded by any descriptor that represents the desired 
absorption parameter, provided the descriptor is ultimately interpretable by a machine and/or 

15 a human. Bar codes, numbers, letters, symbols, scalars and the like are examples. A ranking 
value may be encrypted and later decoded if desired. In some instances, the ranking 
information can be stored in one form or unit of an absorption parameter, and keys utiUzed to 
convert the ranking unit to the next, such as for converting in vitro data to a value tiiat has 
meaning in the context of in vivo data. In other instances the bioavailability data may be 

20 recorded in otherwise abstract pieces, that when combined or processed by a converter, such 
as a mathematical algoritiim, yields a value in a new context. Ranking information is 
preferably stored in a database, i.e., location(s) where data can be permanently and/or 
temporarily stored, accessed and/or updated etc. 

Quality of the data, whether derived fiom the literatiire or otiier sources, is an 
25 important consideration in reliability of the ranking. For example, data in an evaluation test 
set can be assigned or scored with a quality grade, such as A, B, C etc., based on data source 
and quality. Data also may be converted into scalar units for ranking purposes. The scalar 
units can be qualitative and/or quantitative. They also may be assigned to various ranges, 
where values falling wifliin a particular scalar range indicate an extent, presence and/or 
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absence of a particular bioavailability property. For instance, a scalar value in the range of 1- 
10 may be used, where values 1-3 indicate poor bioavailability, 4-6 indicate marginal 
bioavailability, and 7-10 indicate increasing bioavailability. The scalar can be specific or 
general. An example of a specific scalar unit system of ranking is "Scalar Unit = Specific PK 
parameter." General scaling is exemplified by qualitative endpoints (i.e., endpoints simply 
reported as positive (absorption) or negative (no absorption)). Statistical methods can be 
used in the ranking process if desired. 

Generally the selected compounds also are catalogued using a descriptor, which 
reflects the absorption and/or ranking profile. By "cataloging" is intended assignment of a 
descriptor for indexing, filing and/or retrieval fi-om a database. Cataloging permits the 
information to be provided to a compound archive of the originating library, as well as 
separate libraries that contain the selected compounds. This includes cataloging and 
organizing compounds of the library according to a selected absorption profile or specific 
absorption parameter and cross-referencing by features other'than absorption. These features 
include, but are not limited to, for example, distribution, metabolism, elimination, toxicity 
and biological activity. Other features include transport mechanism (e.g., passive diffusion, 
active transport etc.), molecular size (e.g., molecular weight), polarity, charge (e.g.. pKa), 
method of preparation (e.g., synthesis, biosynthesis, extraction, etc.), structure (e.g., mass 
spectrophotometry. X-ray crystallography, NMR, etc.), and/or applications (e.g., 
pharmacologic class, fimctional materials, additives, catalysts, etc.). Altematively, the 
absorption profile/parameter information can be added to existing data files that characterize 
a test sample of a parent library or portion thereof. 

Cataloging is extremely valuable in terms of library management and data analyses of 
libraries. For instance, cataloging faciUtates grouping and organization of compounds in 
high-throughput screening compatible, multi-well units containing single or pools of 
compounds per storage vessel. Data maps outlining the position of each compound or pool of 
compounds and its coiresponding absorption catalog descriptor also may accompany storage 
units. The data maps are represented in a machine-readable format that can be provided to a 
database. 
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High-Throughput Screening: 

The method of the invention U particularly suited for high-throughput screening of . 
compound libranes. High-throughput screening provides a very rapid approach for carrymg 
7a large number of experiments .n a short penod of t.me through the combmat.n of 
. biochemical testing An additional aspect of high-throughput screemng .s the 
5 automation and biochemical tesnng. ♦ „f , tn he greened ' 

„ , for the assavs thus decreasing the amount of a to-be-screeneo 
use of small volumes for the assays, uiu* o .,. • *u dv 

1 of me — Wy P^vides a very ^id approach for canyns ou, a arge „u.b« 
rf«oav.U«U.ypredicSo„si„a.ho«periodof..e.SARandQSAKu,fo™aUonalso.ay 

0 be utilized for Wgh-lhroughput screening. 

A high-O^ughpm screening meUrod can be b^ed on an W,.o assay for collecting 
permeabiUty, solubiUty, dissoluUon and tran^o- mechanism data. The assays car, be based 
ZLX a variety of methods inCnding spectrophotonretric or opbca, ^t,n 
of radioacdve or nuorescen, tnarlccrs, conducHvity, hgh. absorbency, or o*er method 

,5 Itg a molecule based on size, ch^e, afSmty and the like. Any assay method th t .s 
^ ZT^ amenable to — n can be used as the basts of a high-throughput 
:^g rlethod. For «tample, an high-throughput screening method can measure gam or 
3:iioacti,i.y,g.norlossoffluor««nce.achangeintheresis.anceorcon„^^ 

a membrane (WO ,«.372„, cell monoUyers (WO ,7,«7,7, or a cell s„ (WO 
,0 „/49987) High-throughput solid and/or IKpid handling deuces may be employed fcr 
sllmty and dlludon screens. An advanuge of high-.hn>«ihpu. screening -s .he abrbty 
,„ assay a large collection of compounds in an extremely shonperiod of ttme. 

The method of the invention can be inte^ combinatoHal chemistry, 
.^g and instnnnentadon «> ctcate novel assay formats. Uad optimizahon n«y 
« Zloy .rltional single-eompound synthesis or parallel synthesU of dtscrete compotmds^ 
Z fcr parallel synthesis and screening is array chemistry, a system mvoMng 

rZf pamllc. 1-ons for establishing stnrcmre*ioav^«.ity reUtionshtps. Synt^es, 
S^L purtScation e<.ipment can be a«omated and integrate, wi. the PK tool and 
method oftheinv«.tion. This affords faster scree^ngoflargernumbersof compounds. 
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,n particular, compounds may be screened individually or as sets or arrays of sets. 
For instance, a compound set may cont^n „m,ue backbone cbcm,s»y with diverse s,de 
groups attached. These arrays are created by combinatonal reaction of a senal number of 
funcUonally identical bu, structurally diverse building Mocks to create a single compound. 
5 which is Chemically analyzed by HPLC and mass spectrometry. The oompomtds are then 

logically anarrged in spatially addressable multi-well microtiter plates (e.g.. 96 and/or 386 
well mic^titer plates) with a single compound per well. This format yields multtplex »ay 
where oter properties such as chemical structure, mass and the like are stored in a database 
that ,s searchable via a unique catalog descriptor. These arrays may be assembled mto a 
10 larger set of compounds referred to as an array of sets. 

Several strategies can be employed to manage dte high-throughput absorption 
^ng of large chemical compound sets against multiple cell types, tissues, and 
physiological conditions. The single compound per assay per well is the most direct. The 
advantages are that no deconvolution is required and minimal potential for maskmg exists. 
,5 Single compound per assay fits particularly well with the multiplex screening am>y, where a 
prim^ assay provides extensive ab»rptio„ and strucmre4,ioavailability data, the negahve 
assay data also adding value for subsequent lead optimization. 

Alternatively, pooling of compounds per assay well can be employed to quickly and 
efBcieml, assay large compound sets. The primary disadvantage is the need for subsequent 
20 deconvolution of positive readouts, the potential for masking of one compound's ab«,rp.,on 
profile by others, and thus information content of the compound set is parttally lost^ 
Nevertheless, me pcohng approach is very useftrl for rapid high-titroughput screemng of 
larger libraries, where the compounds are screened in blocks containing discrete pools or 
mixnnes of compounds per test sample, where about 1.000 to 100,000 compounds are 
25 represented per block, and about 3 ti, 10 compounds per diser«e pool. Permeabihty and 
solubility d«a are generated using a robotic high-titroughput system. For example, 
compounds can be stored in carousels with robotic retrieval systems mi conveyors dehvenng 
^les to a dispensing are. Individual samples are identiftable by bar codes. The syst«n 
can be configured to hok) a variety of libraries in different phases ranging from solid archives 
.0 liquid libraries in micro-tubes or micro-pUtes and. being modular in de^gn. can be readtly 
adapt^l accoriing to the differing needs of a user and/or piece of testing equrpment. If 



30 



58 



PCT/US99/21151 

WO 00/16231 

de^ir^, d,ssolu..o„ ra,e and transport mechanisn^s preferably screened in subsequent 
rounds of deconvolution. 

PermeaWmy and solubility studies on mixtures of compounds yield hybrid or 
conglomerate value;unless the compounds ^ pn>vided in stoichiometric concenh^tons and 

5 ,„an.ita,ive,y »acted throng the testing process, for instance via LCMS, so as to mnttrntze 
■ masldng effects. Conglomerate values can be separated further by collectmg the relevan 
data for each pool under different gradient tesring eondiUons, such as concentratton, pH. and 
physiological fluid/solvent system gradient conditions. This approach generates pemteab.l,-* 
and solubility profiles representing r^es of values that ate re^tly supplied as input ,n,o the 

1 0 PK tool of the invention. 

Discrete compomtd pools selected by a hybrid absorption ptcfile for a given route of 
administration can menbedeconvolutedby any technique icnowninthe art. smaller pools o, 

individual compounds separated ftom a test pool can be subjected to fUrOter subsequent 
n,„„ds of more focused absotpUon screening accordmg to the methods of the invenhon. 
15 If desired, multiple parallel high-throughput bioavaiUbiUty screens in addition to 

ahsorpdon can be uUli^d to sereen large compound arrays so as to maxtmrze 
stnicture/bioavailability infcrmation content and «unaround time. Tlte compound hhrar,« 
^ed for bioavailability can he integmled into data management programs, and 
multiplexil* in primary and secondary assays. 

20 

As can be appreciated, the method of the invention does not require knowledge of 
biological activity «> cre«e hWes optimi^ for bioavaUability. However, acdvity h.U m 
a^ese libraries are likely U. =thibit hnpn.ved administration route-spedflc bioavatlabtl, y « 
v(» and thus improved »w. phannacological activity compared to those selected by 
25 «:Hvi,y alone. Accordingly, the method of the invention pmvides a receptor-mdependen, 
approach to c«aloging and deigning compound Ubraries with optimally diverse properttes, 
I well as selection and desi^ of compomuls for lead drug development and ophmtzatton^ 
The method is readily adapted for a«oma.ed high-throughput screemng and ranlnng of 
„«! Ubraries (pristine), previously ^ Ubraries (screened), focused hbranes 
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• and selection (focused), or combmauons thereof. Libranes 
pared down by J l,lasethechanceof,den^^^^^^ 

produced by the methods of the invent ^elected route of administration 

having improved pharmacological activity v,vo for a selected 
before or early in drug development. 

~-,.fc nf thi<! invention. These examples 
The following examples illustrate vanous aspects of this 

do not Umit the scope of this invention. 



0 



15 



20 



A 1 f«r «redictinE oral absorption of a compound 
i„ . n^an^al from « v.,™ (e.g.. ^ ^^^^^^^^ 

using a combination of m yitro data and in vivo data. 

,ti.»i model development tool with a graphical user 

smXA* "7'"^; ana a, each s.g= of tt-e «d, as 

«e„ as calculation of flow b«wc«. P ^^^^^^ 3, 
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Example 2: Compound Data Sets 

compound dau sc. for deve,op„,e„,, and .hu. building, ,es,ing, .-ng ^d 
, H „ „L models were obtained from various sources including the hterature and cell, 
rr alTldlman tests as desenbed be.n. The data se. included relevant 
«cal parameters related to absorption of a compound including 0. trac rela ed 
Z PH. initial volnmes, surface area, average transit time, volume transfer rates. 

dissolution, permeability, sotabiUty etc.). 

Bau s«s were selected for compounds that permtted development and isolated 
, ^ and validation for each stage of the build. Compounds suttab.e for this purpose were 
H as follows. For «,e mass, volume ^ int=,.ated mass-volume s,mulat.on mode s 
lidate compound was chosen based on the premise that the best candtdate compound » 
.oae, development would not a dnrg that is highly correlated P— 
between cell, tissue, animal and humans, bu, one ^ is poorly correlated That .^a 
5 compound predicted to have high total absorption h, humans based on pre-clmical smd.es. 
r« exhibited poor ab.on,tion in humans when tes.^ in clinical trials was chos^. 
I^tionally a compound was selected « is no. subject to pre-ahsorptive or hepat, 
:lolism so as to isolate absorption components of the models 6om P-*-n. ^ 3^ 
^etaboUc fac«,rs. Gancydovir „-<,3-dihydroxy-2-pmpoxymethyl)guamne. ,n„„ sodtum 
,0 ^.pHPGJorCytoveneJwassuitablefcrthispurpos. AW,. significant antma, an human 
rJal J was publicly av^-able for Gancydovir Oacobson « An.n..r. .a< A,en. 
W a^ra„. vol. 31, No. 8. p. 1«1-1254 (1987). New Drug /^phcatton Ibr 
Gancydovir Sodium (Syntex. Inc. USA, obtaht^d ftom the Food ^ Drug Attamtstiati . 

25 CImiccI Therapies. (19W) 17:425^32 (1995)). 

For development »,d testing of the integrated GI model, a set of training and testing 
lead drug compounds in vahous sUges of human clinical testing were selected. Thts ^ se, 
included compounds having diverse dosage re,uiremen.s and ranges of permeab-hty. 
solubility, dissolution and transport mechamsms. as shown below in Table 3. 



61 



PCT/US99/21151 



WO 00/16231 



Table 3. 

Comoound Test Set 




al 
a2 
a3 
a4 
a5 
a6 
alO 

Pl 
P2 

P3 

P5 

P6 



4-4-HH- 



-HH-++ 



4- 
+ 



paracellular 
unclassified 
transcellular 
paracellular 
transcellular 
transcellular 
transcellular 
transcellular 
paracellular 
unclassified 
unclassified 




E«unpl.3= Exp«im«.«..D.t.C.U«fl.n«.dPr«essU,g 



? „I for ttainiae »d vaUdadon. exp«m,«.Ul conditions «ere specific 

enough to ensure proper a development 
vart-ions in individn- ^ ,ep«ise 

..„dcdindi™«^.i2:-^^^ , 

""° p„. dissoluSon. d«a poinfs «c« fit » dcenni- 

rrrr - da» ^a-vs. .nd p— ^npio^ed a 

T jTlV^ mode, and weighed .=as..s,ua«s regression analysis (See F>^.. 

,3 :r^r;>""----'-=^^'"^ 
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JpO doses (=™r funC™ input), a hepatic ccmparmten,, as well as an IV dose andfi. 
X e,— eo.pa™,ent. The p,as.a san,p,e is .a.en fro. the een.., eon,p.«n.en., 
and the FDp sample from the hepatic compartment. 
A. Humm/n vivo Data- Oral (PO) 

Plasma levels following oral administration (PO) in humans were used to detennine 
,camol::fcompoundinput.othehepaticvein(n.p,asafrmeUonof«me.Plasma^eve. 

!r L in humans following oral administration of drug solution or suspenston after ^ 
o^^r^t were used as a daU set. If no solutions or ^e^ions were admuus^-. 
3a.ed dosage .rm data were used. ThoPO pro.es i^ludedmch^^^^^^^^^ 
each paUen. enrolled ,n the s»dy from the Sme of administrahon through 24 hours 32 
11 after administradon. along with dosage. If mtadple dose regimens were adm,mstered, 
plasma profiles for all doses were used. 

B. Hum.. /«vi«. Data- Intravenous AdmiBistnitloDaV) 

Plasma leveU following intravenous administration (IV) in humans w^ used to 
d«ermine the amount of drug input to the hepatic v«n (FDp) as a taction of hnte. I V 
p^h.cludedindividualdatapoin.3foreachpatien.enrol,edinthes»d,^m.e^^ 

admit^tion through 24 hours to 32 hours aft. admi^s»hon, along wth the dose, 
multiple dosageregimens were admimster»d.plaamaprofil« for all d^eswereused, 

0 C. /nvtfro Permeability Data 

;„ ^ penneahility da., was used to calculate drug fluxes across various regions of 
^ intestinal mucos. This included rabhi. inteshna. tissue from one or more of duoden,^ 
Smnm. i.»n and colon, ^ C.c.2 cells. T.e mechanism of transport, such as passtv 
C'ular or para^Hul.. canier-mediated .h»n>tion, carrier-medi^ secretion, or mtxed 

„ ^sm, 1 determined for several test compounds and penneahiltnes or eac 
:^smandasse«^aslist^inT.M.4. Protocols for permeahiUt, assays are descnhed 

in Example 4. 
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mechanism permeabilities and parameters for each Gl region. 



wo 00/16231 
Table 4: Transport 

Mechanism 

Passive transcellular 



SX«.io„ » BL w,U,o« inhibition 




Carrier-mediated secretion 



AP to BL and BL to AP 
without inhibition 



Pe 

Kn,, Pc,andPn„orPeat 
entire concentration 
range 

p„,Pe,andPm,orPeat 
entire concentration 
range 



D. Solubility Data 

Solubili.i« of test compounds «s a fanction of pH were determined &om pH 1 5 to 
8 2 in increments of 0.1 pH units. ProtocoU describing condittons for solubfty 
de^rmination are found in M *■ AUematively. solubility a. each pH urn. ftom L5 o 
8 0 was used, with a minimun. of 5 data points a. pH 1.5, 6.0, 6.5, 7.0. and 7.5^ mse 
soiubiUties were used to clcul«e the amount of soluble compound available for absorpuon 
across the intestinal mucosal barrier. 
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E. Dissolution Data 

, „ 5 Protocols d«cnbi„g condUion. for d,sso.u.ion de.=m„na«on fo».d 

'TV r dUso.u.o„ of powd.«. co^po^a, and — y. 
Example 4. a,. fonnulMed dosage fonn «ed to collect oral plaama 

— ™ tltl'^were t^ wl* — - " ^ ' 

C Jfi:— r« — - — — °' 

the intestine. 

Example 4: Protocols for Data Collection 

H helow are detailed protocols utilized for collecting and calculating data 

provided for development of the simulation models. 
A. In vitro Permeability Protocols 
1. Diffusion Chambers 

Peccability data is detccined »si„g intestinal tisane in verHca. difii^ion ch»bc. 
Pem.eab,Uty ^o^.^^^e. or 9.nm round tissne d,ffiis,on 

" Jdonor and receiver buffers witbin .c chancer a. contrnuonsl, 

mixed throughout the experiment. 
,Q 2. Mathematical Calculations 

Effective permeability (Pe) is calculated using Equation 1. 

_ V dC 

(Eq. 1) 



'-ACT* dt 



1- /«,n A ic the surface area available for 
1 ^ tViP receiver chamber (ml), A is me suriau^ 
„bcre V is the votoe of the r«^v^ „ „ ^ fc, , „,„ round and Uw-volume 
fusion (1.78 c,n2 for 8X24 nun ch^nbers, OM^ 
25 chambers), C is the donor concentrabon, ».d dOd. .s calcula 
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^^sion line of tt>e cotecled rece^r co„centta.io„ (see Sampling) v. ume plot- Two 
conditions ™s. be satisfied for .his equaUon .o apply: (D s,.* condiUons in *e recewer 
chamber i e. .he accumulated concen«ion, mus. be virnaally zero when compared ,o .he 
donor conccn.ra.ion; and (2) U-e donor concen.aUon mus, be co„s.an. (C.) ftroughou. .he 

experiment. 

Tl»e parameters for carrier-mediated absorption and secretion are calculated usmg 
Equation 3. 

Pe = 



Co 

1 + Km 



(Eq. 2) 



where Pc is the carrier-mediated permeability. Pm is the passive pem^eability. Km ,s the 
10 affmity of the drug for the earner, and Co is the donor concentration. Pc, Pm. and Km are 
calculated using non-linear regression, Pe is calculated using Equation 1. and Co is gtven as 
part of the experimental conditions. To obtain valid parameter values, Pe is detemimed for a 
sufficient number of Co's to determine Km using Equation 2 (a minimum of 6 Co's .s 
recommended ranging between the analytical Umit and the solubility limit). If Pe values are 
1 5 provided, the variability of the mean as well as the number of experiments performed for 
each concentration are provided to allow accurate regression analysis. 



3. Experimental Conditions 
a. Buffers 



Experiments are performed in appropriate, non-cytotoxic, physiological sahne iso- 
20 osmotic buffers at pH 7.4 (basolatera^serosal side) or pH 6.5 (apicaymucosal side). 
Preferred buffers are Ringer's buffer (pH 7.4), Ringer's with glucose (pH 7.4), MES rmger s 
buffers (pH 6.5), or MES Ringer's with glucose (pH 6.5) (Table 5). 
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Table 5: Formulas for Ringer's buffer and Wger's v».tb gLcose butter. 



Ringer's buffer 
(mM) 

5 
1.15 
0.3 
25 
l.l 
1.25 
qs iso-osmotic 



Chemical 

"kcT 

Na2HP04 
Na2HP04 
NaHCOj 
MgS04 
CaS04 
NaCI 
MES 

Glucose — 

pH adjusted with 1 N HCl or 1 N NaOH 



Ringer's with 
glucose (mM) 

5 

1.15 
0.3 
25 
1.1 
1.25 
qs iso-osmotic 



25 



MES Ringers 
Buffer (mM) 

5 



1.1 
1.25 
qs iso-osmotic 
25 



MES Ringer's 
With glucose 
(mM) 

5 



1.1 
1.25 
qs iso-osmotic 
25 
25 



10 



b. Sampling 

samples are col,.«cd from U>e receiver chamber begirding o,.e steady s,a,e has b«n 
achievl7dco„.iuuing,bra..eas„Omi„u..s.Four.osU«samplesa..l.ec. 

r^low accural «— on of dOd. (E,u...o. 1). The volume r^nove^ .^n, ^ 
Iter Chamber a. each Ume poln. is replaced »i«, buffer con.aining no drug .0 mamum 
ZZ volume in .he receiver chamber. The di,u.o„ of .he receiver co„cen»«on 
r lotion of buffer is corrected during da. analysis and Pe calculaUon. The cone— 
:;Te correct, by: (1, adding d,e mass removed a. each samphng .me .o d.e m. 
IL^ fron, *e receiver chamber a, all prior sampUng .mes, by summmg calcul«^a^.ass 
absorbed and adding » mass for sample calcuMon; and (2, u.„g E,u...n J (preferred). 

/OX n-l 



1 . p 



(Eq.3) 



v^here d,e co«c.ed receiver chamber c„ncen««on U obutaed by dividing .he co„e«ed 
.5 J^c concen^a^on by E,uaa.n 3 (l/X), S .s .he volume of sample wiU...awn V ,s ^ 
Zer Chamber volume, . is U,e se,uen«a, sample number. i.e.. for «,e fi^ samp. 

A c«n,nie time k-3 for the third sample time, etc., and (5 is the 
time, k=2 for the second sample time, k j ror 

corresponding number from Pascal's triangle (Table 6). 
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Table 6: Pascal's Triangle for determining p coefficients, 

4'" term 5'" term 6 ' term 




Donor concentration (Co) is determined by sampling the donor buffer containing the 
test compound with subsequent analysis directly from the donor chamber, or from a stock 
5 solution of donor buffer provided binding and absorption to the interior of the chambers does 
not occur. 

c. Intestinal Tissue 

Rabbit intestinal tissue is used for permeability experiments. During mounting of 
tissue onto chambers, intestinal muscles are stripped off the mucosa and discarded. Care 
10 should be taken to ensure integrity of the tissue. A minimum of three chambers are used to 
determine Pe values for each region, concentration and compound. The mean Pe and 
Standard Error of the Mean are provided for each study. 

d. Cell monolayers 

Caco-2 cell monolayer Pe is determined in diffusion chambers similar to NaviCyte 
15 Snapwell™ diffusion chambers and follow all procedures described above except the 
recommended buffers are Ringer's with glucose or MES Ringer's with glucose as listed m 
Table 6. 

Caco-2 cells are grown using DMEM media supplemented with 10% FBS, 5% PCN- 
STEP and 1% NEAA under 95-100% humidity and 5% CO3 at 3rC. Cells are grown m 
20 flasks and the culture split at 85-95o/o confluence. Snapwells- are seeded at 65,000 cell/cm 
and used in the permeability experiment within 21-28 days post seeding to allow for 
differentiation. 
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4. Determination of absorption mechanism 

Absorption mechanism for a compound is determined by one of the following 
methods Determination of Pe in both the apical-basal (AB) to basal-lateral (BL) and BL to 
AB directions using Equation 1, or determination of Pe in the AB to BL direction at 
5 concentrations, (a) close to the analytical limit, and (b) close to the solubility limit. 

Similar Pe values in both the AB to BL and BL to AB indicate a passively absorbed 
compound and no fiirther studies are required. AB to BL Pe greater than BL to AB indicates 
caiTier-mediated absorption and Pe must be detennined for 5 additional Co in the AB to BL 
direction. BL to AB Pe greater than AB to BL indicates carrier mediated secretion and Pe 

10 determined for 5 additional Go's in the BL to AB direction. 

Similar Pe values at low and high concentrations indicate a passively absorbed 
compound, and no further studies are required. Low concentration Pe higher than high 
concentration P. indicates carrier-mediated absorption and Pe is detennined for 5 additional 
Go's in the AB to BL direction. High concentration Pe higher than low concentration Pe may 

15 indicate carrier-mediated secretion. BL to AB Pe is then detemiined at the low concentration 
and the mechanism determined as described above. 

B. Solubility determination 

Solubility of a compound is determined using an accurate and scientifically sound 
method similar to the Phase Rule and Phase-solubility analysis as described in Remington's: 
20 The Science and Practice of Pharmacy, 19* edition, Chapter 16. 

The solubility is determined at pH 1.5 using Simulated Gastric Fluid (USP XXH) 
minus pepsin. Solubility at pH 6.0, 6.5, 7.0, and 7.5 is determined in Simulated Intestinal 
Fluid (USB xxn) minus pancreatin. Parameters are for data collection are carefully 
monitored by ensuring purity of the test compound and accuracy of the Simulated 
25 Gastrointestinal fluids. A temperature of 37°C is maintained accurately during the course of 
the determination. Complete saturation and accurate analysis of saturated solutions are 
employed. 
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C. Dissolution determination 

The dissolution rates are determined using the equipment, apparatus, and methods 
described in USP XXII, <711> dissolution. The dissolution rate at pH 1.5 is determined in 
Simulated Gastric Fluid (USP XXII) minus pancreatin. Concentrations are collected and 

5 analyzed for drug compound from the vessel for a sufficient time (6 hours, preferable) to 
allow the initial slope of the concentration v. time curve to be detemined. The slope 
(dissolution rate) is determined using the initial linear portion of the concentration v. time 
plot if non-sink conditions exist. Under sink conditions, the entire plot are used to calculate 
the slope. The slope is reported as the dissolution rate. Explanations of the dissolution rate, 

10 sink and non-sink conditions, and equations for calculation are given in Remington's: tiie 
Science and Practice of Pharmacy, 19th edition, Chapter 34. 

If a formulated dosage form is used for dissolution testing, the dissolution protocols 
described are used to determine tiie dissolution rate for drug compound from the formulated 
dosage form. 
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Example 5: Standards and Protocols for Evaluating Permeability Data Collection 

This example provides detailed protocols for controlling the quality of permeability 
data coUection described in Examples 3 and 4. Compounds listed in Table 7 are used as 
standards for monitoring permeability data collection and quality. The compounds were 
chosen to represent each intestinal transport mechanism (passive transcellular, passive 
paracellular, carrier-mediate influx, or carrier-mediated efflux). 

Table 7: Permeability Standards 



Transport mechanism 


Compounds 


Passive Paracellular 


maxmitol 


Passive Transcellular 


hydrocortisone 


Carrier-mediated Influx 


D-glucose 


Carrier-mediated Efflux 


etoposide 



Mannitol, hydrocortisone, D-glucose, and etoposide also were chosen since they are 
25 widely used as markers for intestinal transport across rabbit tissue and other systems with 
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well characterized Pe values. These compounds also are available commercially as either 
3H-labeled or 14C-labeled. 

Permeability data for standards is compared to the values for rabbit Hsted in Table 8 
(or other standard values) using basic statistical analyses. If the data is significantly different 
5 (p-value>0.05) for any of the standard compounds, data collection is repeated. 

Table 8: Transport Characteristics of PermeabUity Standards* 

Compound (donor 

concentration) Duodenum Jejunum Ileum Colon 

mannitol(lmM)5 1.73 xlO'^ 3.MxlO^ 4.02 xiU^ 5.53x10- 
hydrocortisone 3.00x10' 1.31x10-* 2.91x 10 3.85x 10 

Hu^eOO 4.55x10-^ 1.02x10-^ 1.45x10-^ 9.28x10-* 

mM)5 

etoposide(lOO^iM) ^ ^ 

*Note: permeability values are representative of ranges. Other values or extended ranges 

may be used. 

A. Experimental Conditions 

10 Protocols, conditions and calculations for permeabiUty evaluation of standards are as 

described in Example 4, with the foUowing modifications. 

Permeability experiments are performed using Ringer's buffer at pH 7.4 on both the 
apical/mucosal and basolateral/serosal sides. Ringer's buffer is as described above excepting 
that glucose is substituted with mannitol when Pe values for glucose are being measured. 

15 Samples are collected from the receiver chamber beginning 30 minutes after 

experiment initiation and continuing every 15 minutes until 6 samples have been collected 
(105 minutes). One-half ml is removed from each receiver chamber at each time point and 
compound concentration detemMned. The volume removed from the receiver chamber is 
replaced with buffer containing no drug to maintain constant volume in the receiver chamber. 

20 The dilution of the receiver concentration due to the addition of buffer should be corrected 
during data analysis and Pe calculation. The concentration is corrected by using Equation 4. 
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(Eq. 4) 



Where the corrected receiver chamber concentration is obtained by dividing the collected 
sample concentration by Equation 4 (1/X), S is the volume of sample withdrawn, V is the 
receiver chamber volume, k is the sequential sample number, i.e. k=l for the first sample 
time, k=2 for the second sample time, k=3 for the third sample time, etc., and p is the 
corresponding number bom the modified Pascal's triangle below (Table 9). Note: Since the 
sample intervals are not even (i.e. the 1st interval is 30 minutes, all others 15 minutes) 
Equation 4 as well as the p coefficients are modified from those listed in Example 4. 

Table 9: Modified Pascal's Triangle for determing p coefficients 



Sample 


1st term 


2nd term 


3rd term 


4th term 


5'" term 


6th term 


1 


2 












2 


3 


2 










3 


4 


5 


2 








4 


5 


9 


7 


2 






5 


6 


14 


16 


9 


2 




6 


7 


20 


30 


27 


11 


2 



The donor concentration Co is determined by sampling 0.02 ml of the donor buffer 
containing drug (with subsequent analysis) directly from the donor chamber. Potential 
binding of drugs to the chambers also is monitored. Donor samples (0.02 ml) are taken at 
experiment initiation and at experiment conclusion. If a significant decrease in drug 

15 concentration has occurred (>10%) the experiment is repeated using procedures which 
compensate for the drug loss in the donor chamber. It is recommended that the donor 
chamber solution be removed and replaced with fresh donor buffer containing drug at 
appropriate intervals. The intervals and volumes to be used are determined using sound 
scientific judgment. Adequate data is collected to show the donor drug concentration has 

20 remained constant throughout the experiment. 

For tissue-based permeability assays, during mounting of tissue onto chambers, 
intestinal muscles should be stripped off the mucosa and discarded. Care should be taken to 
ensure integrity of the tissue. 
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Animals donating tissue are euthanized immediately prior to expenment mitiation. 
The small intestine is excised from the animal and kept in ice cold Ringer's buffer pH 7.4 
until mounted in diffusion chambers. As soon as possible after excision, the tissue is cut into 
an appropriately sized piece and placed over the diffusion chamber pins with the mucosal 
5 side down. The muscle layers are carefully stripped away using forceps. After the tissue is 
mounted the two half chambers are placed together and the donor and receiver sides filled 
with the appropriate pre-warmed (37°C) buffer. If NaviCyte chambers are used, the gas hfl 
system is comiected with 95% 0,15% CO. flowing at ~5 - 15 ml/min (depending upon 
chamber volume) into each half chamber to maintain pH and mixing. Sampling begins 30 
10 minutes after connection of the gas lift system. 

The mean Pe and Standard Error of the Mean are determined for each study. 
Permeabilities from at least 6 chambers from 3 different animals are used in calculating the 
mean and Standard Error of the Mean. 

In addition the Pe of radiolabeled mannitol is determined simultaneously with the 
15 standard compourid as a marker of intestinal integrity. Mannitol Pe values may be 
determined by concurrent diffiision using a donor buffer containing mannitol and the standard 
drug compound, or by continuing the experiment for 60 minutes after the last standard 
compound sample is collected using donor buffer containing mannitol and fresh receiver 
buffer containing no compounds. 
20 Special experimental conditions are followed for certain standard compounds. This 

includes such conditions as a proton gradient, a sodium gradient, presence of glucose, etc. 
These conditions are Usted in Table 10 and are substituted or added to the general conditions 
listed above. 



Table 10: Experimental Conditions 



Standard Compound 



Donor 
Concentration 

ImM 
lOmM 
0.01 ^M 
100 jiM 



Special Conditions 
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Example 6: Physiologic-Based Mass Simulation Model 
A. Design 

A multi-compartment physiologic-based mass simulation model (the "mass-model") 
was designed to integrate mass-flow relationships among GI compartments representing the 
stomach, duodenum, jejunum, ileum, and colon, and thus throughout the Gl tract, and to 
characterize drug movement in units of mass into peripheral compartments. Converters that 
interrelated transfer rates and associated rate constants (kt), which in turn were modified by 
various factors including pH, solubility profiles, compartment surface area and drug 
permeability were incorporated to account for drug movement among compartments. A 
plasma kinetics model also was included for validation purposes and for correlating clinical 
plasma data to the mass model. Converters also were used for unit conversion. 

Gancyclovir was chosen to develop and test the mass model. Gancyclovir exhibits no 
in vivo biotransformation and is pooriy absorbed. Thus, the mass model assumes no 
metabolism or protein binding. Additionally, dissolution rate and delivery system were not 
used in the mass model as modifying parameters of drug absorption, i.e.. drug assumed to be 
completely dissolved in the stomach and solubilized according to its solubility profile. 

Surface area values for each compartment of the mass model represented a 
"functional surface area." as opposed to an absolute value. A functional surface area was 
utilized since (1) fluids entering the gastrointestinal compartments do not cover the surfaces 
of the compartment instantaneously, but rather over a time course; and (2) solubilized drug 
within the fluid is not ideally presented to aU absorptive areas. Functional surface areas for 
each compartment were calculated by solving Equation 5 for the area using various data 
inputs fi-om the literature. 

p»A»Sp=aM/at 

(Eq. 5) 
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Where P is the permeability coefficient, A is the surface area of the membrane, Sp is the 
solubihty of the drug in the relevant segment of the intestine, and 9M/9t is drug flux, where 
flux dWdt is determined from the permeability of the drug in the particular intestinal 
compartment, the surface area covered by drug solution and the solubility of at the pH of the 
5 intestinal compartment. 

For example, several studies have been conducted comparing permeability of various 
compounds (Rubas et ai, Pharmaceutical Research, Vol. 10, No. 1 (1993)). Mamiitol, 
which has similar physicochemical properties to Gancyclovir, also has similar permeabiUty 
characteristics and a bioavailability of approximately 10% in humans when it is orally 

10 administered. For mannitol, permeability is well characterized. Thus, data obtained from the 
Uterature related to permeability in each compartment, pH-dependent solubility and mass 
concentration relationships was used to solve Equation 5 for area. Thus, it was this area, and 
not the theoretical total surface area of each compartment, that was used as the functional 
area of a compartment, which represented a good approximation of in vivo surface area 

15 relationships for initial model building. 

Permeability values were obtained from published in vitro cell diffiision experiments 
and were accounted for by converters that modified luminal and peripheral flow (K12) for 
each compartment. For solubility, a solubiUty curve was used based on experimental data 
available in the literature. pH was then isolated in a separate converter to modify the 
20 solubility curve for the particular compartment. In contrast, for validation purposes, an 
absolute solubility value was used and pH was entered as 1 to isolate that converter from the 
validation model. 

Absorption "transfer" rates among each two compartment sub-system were collected 
into a separate flow representing total absorption rate, which in turn was collected into a 
25 compartment representing the total amount of drug absorbed for each GI tract compartment, 
namely, stomach, duodenum, jejunum, ileum, and colon. Absorption rates among stomach, 
duodenum, jejunum, ileum, and colon modules were connected by flows modified by the 
associated rate constants between each GI segment. 
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For validation purposes, a plasma kinetics model was integrated with the mass-flow 
compartments by linking the total absorption rate to a flow representing the absorption rate 
constant, which in turn fed into the central plasma compartment. A standard two- 
compartment plasma kinetics model (Ramsay, European Journal of Pharmaceutics and 
Biopharmaceitucs, Vol. 37, No. 3 (1991)) was used for this purpose. (See Figures 9 and 10) 
The plasma kinetics model incorporated first order transfers between the blood compartment 
and peripheral compartment. Two flows were used and set up as first order systems and thus 
different rate constants were applied in each direction. Compartment values were represented 
as mass units. Blood volume was input in a converter, which modified a converter for 
concentration along with the mass compartment. An elimination rate constant was also 
obtained form the literature in a first order process. In addition, while most drugs are given 
in milligram doses, plasma concentrations are reported in microgram or nanogram per 
milUliter. This is done since compounds are distributed rapidly into a large volume after 
entering the blood resulting in a concentration of drug in systemic circulation that is quite low 
with respect to the concentration at the site of administration. Accordingly, an additional 
converter was added to convert milligram units to nanogram or microgram units expected for 
concentrations of the test compound based on human bioavailability data. A compartment 
also was added to collect elimination data. 



B. Mass Model Parameters 

Parameters and associated values of the mass model include pH, solubility, 
permeability, and intestinal transit, and are illustrated in Table 11. 

Table 11: Mass Model ParametersA^aiues 



Parameter 


Value 


Dose 


lOOOmg 


dt 


0.125 


RunTime 


24hrs 


ka assumed (mass transit) 


2.8 or 3 




Stomach 




Area 


50 cm" 


Solubility 


31 mg/ml 


Permeability 


1.1 X 10"" cm/sec 


Duodenum 





76 



in 

W 
g1 



W 



wo 00/16231 



PCT/US99/21151 



Area 


125 cm' 


SoluDility 


3.65 mg/ml 


rermeaDiiuy 


1.1 X 10"" cm/sec 


Jfejunuiiin 




Area 


182 cm" 


DOiuoiuty 


3.65 mg/ml 


renneaDiuiy 


2.17 xlO*" cm/sec 






Area 


102 cm" 


oOiuDiiiiy 


3.65 mg/ml 


Permeability 


4.06x 10"" cm/sec 


Colom 




Area 


138 cm" 


Solubility 


3.65 mg/ml 


Permeability 


3,80 10*" cm/sec 




Plasma Kiaetics 






0.839 


k2i 


0.670 




0.161 


Fluid Volume 


76,800 ml 



m 

in 
o 



The mass model also was tested by inputting values derived fipom the literature 
(Gibaldi et al. Pharmacokinetics, pp. 284-288, Marcell Dekker (1975)) into the plasma 
kinetics model. These values are shown in Table 12. 



TfflMe 12: Values for Plasma Kinetic Modole 



Dose 


Ig 


1505a 


2.718 h" 


1505b 


0.254 h■' 


k2i 


0.3737h■■ 


ki2 


0.7509h-' 


kio 


1.3474h-' 


v„ 


20.1241 
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Example 7: Testing and Validation Mass Model 

The mass model was tested using parameters shown in Table 11 with an initial dose 
of 1000 mg over a time course of 24 hours. AUG, Q„ax. T„.x, and T,/2 were simulated using 
various doses (New Drug Application for Gancyclovir Sodium, Syntex (USA), (obtained 
from the FDA under the Freedom of Information Act (FIA)) and compared to human clinical 
data obtained for Gancyclovir. Bioavailability simulated by the mass model for Gancyclovir 
was approximately 6%. Compared to human clinical data, obtained for two Phase I clinical 
studies (designated here as ICM 1505 and 1505b). bioavailability of fasted patients in clinical 
trials typically ranged from 3-20%. The mass model also was tested using a plasma kinetics 
validation model illustrated in Figure 9. 

Figure 17 shows the area under the concentration time curve for a 1000 mg dose of 
Gancyclovir, Tmax = 1.4 hrs, Cmax = .51 ng/ral., using the mass model, as compared to 
clinical study data of ICM 1505 and 1505b. The results demonstrate that the mass model 
underestimated plasma concentration during the post-absorptive period. Table 13 shows 
comparison of some values between cUnical studies and those predicted by the mass model. 
The clinical studies also used a 70Kg body weight for normalization of concentrations. 



Table 13: Comparison of Mass Model to Clinical data 



Parameter 


Mass Model 


Clinical 1505a 


Clinical 1505b 


Cmax (mcg/ml) 


0.51 


0.55 


0.59 


Tmax (hrs) 


1.40 


1.43 


1.43 



Example 8: Physiologic-Based Volume Simulation Model 
A. Design 

A physiologic-based simulation model for incorporating fluid volume flux and GI 
transit (the •'volume model") was developed for integration with the mass model to account 
for changes in absorption resulting from fluid absorption/secretion and transit, and thus 
appartmt drug concentration. The volume model was constructed so that fluid enters a 
compartment and was absorbed by a first order process based on an absorption rate for that 
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Of fluid between compartments was dependent on a zero or first order fluid 



B. Volume Model Parameters 

As a starting point for the volume model, values were obtained from literature that 
described in general terms absoiption and secretion of fluid throughout the body (Change 
al Gastrointestinal, Hepatobiliary and Nutritional Physiology, Chapter 5. p. 92, Lippincott- 
Raven (1996)). Values representing total intake of fluid per day and total secretion of fluid 
per day were modeled into the system normaUzed linearly to increments of dt for the model. 
To permit for changes in dt for the model, the values were entered as pulses. Values used in 
the volume model are shown in Table 14. 

Table 14: Volume Model ParametersA^alues 



Source 
Intake/Secretion 


ml/24brs 


ml/0.1 hrs 


Stomach 


6500 


27.08 


~~ Orally 


2000 


8.33 


Salivary 
Glands 


1500 


6.25 


Stomach 


2500 


10.42 


Duodenum 


2000 


8.33 


Bile 


500 


2.08 


Pancreas 


1500 


6.25 


Jejimimfi/Ileum 


1000 


4.17 


Jejunum 


641 


2.67 


Ileum 


359 


1.50 


Colon 


0 


0 


Total 


9000 


337.57.5 


AbsorptioD 






Duodenimi 


2598 


10.82 


Jejunum 


3783 


15.76 


Ileum 


2120 


8.83 


Colon 


400 


1.67 


Total 


8900 


37.09 


Note- Values for coirpartments based ou %lolal uilcilmal aicu 



Where data was only available for a series of compartments, values were assigned to 
each compartment based on the percentage of the total area for that series (e.g. secretions for 
and ileum and absoiption for parts of the small intestine). The model was set as two 



15 jejunum 
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flows between the blood (serosal) side of the compartment and the compartment itself. Each 
flow represented the rate constant for secretion and fluid absorption. 

For development purposes, absorption and stomach secretion were assumed to be zero 
order when using values fiom Table 14 for both flows. Also, daily volume for fluid entry 
into the stomach was entered as a pulse according to the dt values shown in Table 14. Thus, 
total intake and secretions of fluid was modeled as a pulse occurring every 6 minutes 
throughout a 24 hour period. Initial volume in the stomach also was set up as a pulse of the 
total oral intake, salivary excretion, and stomach secretion over each dt increment. 



Example 9: Testing and Validation of Volume Model 

To test movement of fluid between compartments the volume model was modified to 
approximate zero order fluid transit or emptying and isolated 'ftom the mass component of the 
model. Initial values of 1000 ml and 250 ml were used for testing. 



Example 10: Physiologic-Based Mass-Volume Simulation Model 
A. Design 

A physiologic-based simulation model integrating the mass and volume models (the 
"mass-volume model") was constructed to integrate complex mass and fluid flow 
relationships. The integrated mass-volume model also included compartments to characterize 
drug movement into peripheral compartments. A plasma kinetics model for 
training/validation purposes also was included. The basic design for the integrated mass- 
volume model, Unked to the plasma kinetics model shown in Figure 9, is illustrated in 
Figure 12. 

Volume for a compartment was added as a product to obtain the amount of drug 
solubilized at a time increment volume. Additionally, an "IF ... THEN ... ELSE" control 
statement was added to prevent the equation from indicating that more drug was solubilized 
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than dosed. Thus, the integrated mass-volume model shows the mass of drug in the stomach 
connected to the absorption rate constant as well as the volume compartment. 

Mass and fluid transit rate constants of 2.8 and 3 for the stomach were calculated from 
values obtained from the literature for Gancyclovir (Syntex, Clinical Studies ICM 1653 and 
1774, FDA NDA available data and Bachrach et al. Functional and Diagnostic Aspects of 
the Upper Digestive Tract, Digestive System. Part I. Upper Digestive Tract, Netter (1989)), 
and determined for each of the remaining compartments to approximate mass and fluid 
movement. 

B. Mass-Volume Model Parameters 

Parameters and associated values and equations were systematically varied or as 
described above for individual mass and volume models; an example of the equations and 
parameters employed in the mass-volume model are shown in Appendix 1. Dissolution rate 
and delivery system (controlled release device/formulation) were excluded from in the mass- 
volume model, and thus the model assumes a test compound is immediately in solution in the 
stomach. 

Example 11: Testing and Validation of Mass-Volume Model 

The mass-volume model was tested using the equations and parameters shown in 
Appendix 1. These parameters included the pulsed estimate of fluid absorption and 
gastrointestinal secretions, and rate constants extracted from the Uterature. Alternate sets of 
parameters for fluid absorption and secretions also were tested. For example, simple zero and 
first order rate constants of 1 or a sequential integer and various doses were evaluated for 
comparison to human clinical data. 

Figure 18 shows the area under the concentration time curve for a 1000 mg dose of 
Gancyclovir, Tmax = 1.1875 hrs, Cmax = .54 mcg/ml., using the mass-volume model of 
Kgure 12 with the estimated absorption and secretion rates, relationships, and values of 
Appendix 1, as compared to clinical study data of ICM 1505 and 1505b. The data is now 
less favorable for Tmax but more favorable for AUC compared to the mass model. These 
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results demonstrate that the mass model underestimated plasma concentration during the 
post-absorptive period, while the combined mass-volume model appeared to overestimate it. 

The mass-volume model was modified to incorporate simple zero and first order 
absorption and secretion. This model was then run using an initial volume of 250 ml and also 
4 administrations of 250 ml water as done during clinical studies. Results were similar to the 
results shown in Figure 18, but with slightly higher absorption. 

The mass-volume model also was run using the following combinations of data input: 
(1) doses of 500 mg, 750 mg, 1000 mg at qid, bid, and tid dosing; (2) initial volumes of 250 
ml, 500 ml, 1000 ml; (3) varying absorption and secretion rates based on differing 
assumptions for daily secretion and fluid intake; (4) vatying pH values in the various 
compartments; and (5) simulation of food intake and fasting conditions. Correlation was very 
good with some clinical data and less than optimal with others. Correlation with theoretical 
estimations also varied from very good to poor. 

Collectively, the mass-volume model represented an improvement over the individual 
mass and volume models in that it provided a better approximation of in vivo conditions. 
While the simpler mass-model correlated better with clinical data, the integrated mass- 
volume model was more sensitive to changes in the various input parameters, physiological 
conditions and underlying constants, and thus a more rigorous model of the GI tract. 



Example 12: Physiologic-Based GI tract Simulation Model 
A. Design 

The mass-volume model was selectively improved in a stepwise fashion to create an 
integrated physiologic-based simulation model of the GI tract of a mammal (the "GI model") 
capable of compound-independent prediction of oral absorption with a high level of accuracy. 
The model was developed to be flexible. That is, it was designed so that additional 
physiological factors that influence oral absorption could be identified and incorporated into 
the model as needed to improve the quality of the prediction for a diverse set of test 
compounds. Additionally, the GI model was developed to minimize input data requirements. 
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The basic approach involved generation, testing and integration of a GI transit model 
(Figure 21), a pH-dependent solubility and dissolution model (Figure 22), and an absorption 
model (Figure 23), as well as underlying equations and parameters, constants, calculated 
parameters, and rules by which a given simulation is to proceed. A controlled release device 
and formulation compartment also was included. A graphical compartment-flow model of 
the integrated GI model is illustrated in Figure 25 (without converters, ghost or connectors) 
and Figure 26 (with converters, ghost and connectors). Parameter inputs, calculations and 
outputs are illustrated in Figures 30-40. An abbreviation key for the GI model is provided as 
Appendix 3. 

The GI model also incorporated additional features to improve the predictive power 
and versatility of the simulation model. One feature was the development and incorporation 
of regression analysis derived adjustment parameters based on analysis and processing of 
human clinical data and in vitro data for a diverse set of compounds. The adjustment 
parameters were utilized as constants in the GI model, and thus modify underlying equations 
of the model. A second feature was development and incorporation of regional permeability 
correlation parameters and equations that permitted estimation of values for segments of the 
model that were missing user provided input values for corresponding parameters. This 
facilitates prediction of oral drug absorption when permeability values or other parameter for 
a given compoimd are provided for a to limited number of GI segments, for example, when 
cell-based input data, such permeability data derived from Caco-2 cells is used to provide 
permeability input data of colon. Another feature was development and incorporation of 
parameters and calculations to account for transport mechanism and thus transport-specific 
variations in compoimd absorption. Another feature was incorporation of the abiUty to 
isolate and evaluate specific regional absorption events related to dissolution and mass 
transit. Also, the GI model was developed to separate absorption into the portal vein (FDp) 
from hepatic metabolism, so as to account for individual primary barriers to absorption. 

B. GI Model Equations, Rules and Parameters 
L General Equations For GI Model: 

Various differential equations and rules utilized for the GI tract model are provided 
below. For the equations, adjustment parameters are designated by the letter Z. 
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Transit time: 

order tr ansit process 

^ = k^[A] (Eq. 6) 

dt 

dA/dt = rate of transit (or absorption), krr = rate of constant, A = amount (compound or 
water) in proximal compartment. 



H^te consi ^n^ r.aiculation 



TTadj = adjusted transit time 

TTadj = (TTp 'Zrr ' Vserrr) (Eq- 8) 

TTp = physiological transit time, Zrr ^ transit time adjustment parameter, Userir = User 
controlled adjustment to transit time. 

Krr is a regionally dependent parameter, i.e. different rate constants are used for each region 
of the GI tract. 

Fluid volume absorption/resorption: 

^ = U^] (Eq.9) 
at 

dA/dt = rate of absortpion, kvA = rate constant, A = amount of fluid (water) in the 
compartment 

ZvA = volume absorption adjustment parameter, kemp is determined emperically to match 
human fluid absorption in vivo. 

Dissolution and Solubility: 
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dt 

A = Amount dissolved, ko = User supplied dissolution rate constant, Zd = Dissolution rate 
adjustment parameter, Sadj = solubility, C = concentration 

5 Snliihilitv (regionally dependent) 

iP"n-P">i-\) 

Sadj = Solubility. S„ = user supplied solubility {S.-Ss) , pH„ = user supplied pH values 
{pH,...pH5} corresponding to user supplied solubilities, pH = pH value appropriate to region 
of the system, such as GI tract, n is selected such that pHn > pH, and pH„., < pH. If any of 
10 pH,...pH5 are equal to pH, the corresponding Sn is used as the solubility. 

rnncentration (regionally dependent) 

r = ^m. (Eq. 13) 

V 

C = concentration of soluble drug, V = volume of fluid 
Flux/Absorption: 

J-^Paixj.SAadj.C (Eq.l4) 



15 



J = flux, Padj = Adjusted permeability, SAadj = Adjusted surface area available for 
absorption, C = concentration 



^ADJ ~ 



f^Uz, -3600*^^2^2^ (Eq.l5) 



Zeff = Efflux transport adjustment parameter. Pm = passive membrane permeability, Zf - 
20 passive permeability or flux adjustment parameter. Zact = active permeability adjustment 
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parameter, Pc = active carrier permeability, C = concentration, K^, = Michaelis-Menten 
kinetic parameter. 

RRginnal Pen TTrjjhflitv nnrrelation 

Any regional permeability, Pm, can be calculated using any number of other provided 
permeabilities. 



[\2 



(Eq. 16) 



Pa = permeability calculated using the regional correlation, Pb, permeabihty provided by the 
user. A, B, and C = correlation coefficients fitted to determine correlation. 

By way of example, rules utilized for a GI tract model of the PK tool and method of 
the invention include the following general processes. 

2. General Processes For Rule Generation: 

1 . GI transit. The transit of drug compound and fluid volume are somewhat controlled 
and the transit of formulations and/or controlled release devices is much more strictly 
controlled. 

2. Controlled Drug Release. The release of drug from the dosage form must be 
controlled such that drug is released into the correct intestinal region at the 
appropriate time. 

3. Dissolution. A comparison between the concentration and the solubility must be 
made to determine if additional insoluble compound will dissolve, or if compound 
already dissolved must precipitate to insoluble drug due to solubility limitations. 

4. Absorption. Mathematically, absorption may occur when physiologically it is 
impossible, e.g. when the volume in the colon becomes low enough that any dissolved 
drug will be within fluid contained in other solid waste also present in the colon and 
therefore unavailable for absorption. IF...THEN production rules control these 
situations. 
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5. Permeability calculations. To estimate unprovided permeability values from provided 
permeability values logical evaluations must be made to determine the correct 
equations necessary to make the correlations. 

6. Concentration calculations. The concentration in the intestine cannot exceed the 
solubility for that particular region. If it does, an incorrect flux will be calculated. 
IF...THEN production rules are used to ensure the correct concentration is used in the 
flux calculation. 

7. Mathematical anomalies. At certain times during the simulation (especially early and 
late in the simulation) some compartments, flow regulators, or converters used in 
other calculations may have a value of 0 which will result in a computational error, 
e.g. division by 0. Production of rules are used to identify these situations and avoid 
the errors. 

The following table lists the specific processes, conditions, results that control 
statement rules, e.g., IF...THEN production rules, are used to control. Generally, separate rules 
used for each region of the GI tract and are combined into one line in the table. 
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Table 15: Rules for Physiologic-Based GI tract Simulation Model 



Process 


Condition 


Result in True 


Rcsuh if False 


Comments 


GI Transit of drug 
compound or fluid 
volume 


Time < 4 hours 


No transit to 
waste 


Transit to waste 
by fu-st order 
process 


Applies to GI 
regions using 
different values 

lUl LUC k^UIlUiUUii. 


GI Transit of 
formulations or 
controlled release 
devices 


Time, cumulative 
physiol. transit 
time 


no transit to next 
compartment 


Immediate transit 
to next 
compartment 


The rate constant 
for first order 
transit is set 
exceeamgiy large 
to provide near 
instantaneous 
transit. 


Controlled release 


Time to reach GI 
region < Time < 
Time to exit GI 
region 


Drug is released 
from dosage form 
to GI region 


No drug release 
into that GI region 


Drug is released 
according to user 
proviucu rcicodc 
profile. 


Dissolution 


Soluble 
drug/volume 
(concentration) < 
Solubility 


Dmg moves from 
insoluble to 
soluble 
compartment 
according to 
dissolution rate 


Drug moves from 
soluble to 
insoluble 
compartment 
according to 
precipitation rate 


Precipitation rate 
is set to provide 
near instantaneous 
precipitation 
without causing 
"overshoot". 


Absorption 


Volume<lxlO"^ 
ml AND Mass < 1 
xlO**mg 


No absorption, i.e. 
concentration = 0 


Absorption by 
flux equation 




Permeability 
Calculations 


Duodenum, 
Jejunum, and 
Ileum 

Permeabilities all 
provided 


Use provided 
Permeabilities 


Estimate 
unprovided 
penneabilities 
from provided 
permeabilities 


1 or2 

permeabilities can 
be used to 
calculate 
unprovided 
permeabilities 


Concentration 
Calculation 


Concentration < 
Solubility 


Concentration 
used in flux 
equation 


SolubiUty used in 
flux equation 




Mathematical 
anomalies 


Volume = 0 


Dissolution rate = 
0 


Dissolution rate 
calculated by 
Noycs-Whitney 
equation 


Dissolution given 
as an example. 
Similar 

condications are 
provided for 
concentration 
calculations and 
other processes. 
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Exemplary equations, rules, parameters and initial values for the graphical 
compartment-flow model and various sub-models of the integrated GI model illustrated in 
Figures 21-26 and 30-40 are provided in Appendix 4, as related to the abbreviation key 
provided as Appendix 3. Various aspects of the physiological, adjustment and regional 
correlation parameters employed in the GI model and their development are described in 
further detail below. 

1 . Physiological Parameters 

Physiological parameters of the GI model included physiological ranges reported in 
the literature (Table 17) as well as specific values utilized in the model and compiled for 
each of five regions of the gastrointestinal tract (stomach, duodenum, jejunum, ileum and 
colon)(Table 16). These included values related to pH, transit time, surface area, and volume 
parameters. 

Table 16: Physiological Parameters Employed In GI Model 





pH" 


Initial 
Volumes 
(ml) 


Surface 
Area 


Average 
Transit 
time 
(hr)^ 


Volume 
Transfer 
Rates (t9o) 
(hr'r 


New Water 
Absorption 
Rates* 
(hr')** 


Stomach 


1.5 


100 


NA 


0.5 


4.6 


0 


Duodenum 


6.0 


0 


150 


0.225 


10.8 


0 


Jejunum 


6.5 


0 


1000 


1.5 


1.54 


1.75 


Ileum 


7.0 


0 


1000 


1.5 


1.54 


1.75 


Colon 


6.5 


0 


850 


24 


0.094 


0.1 



♦Water absorption rate parameters were set so that cumulative water absorption from each 
region using the GI model were in agreement with values listed in Table 17. 
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Table 17: Physiological Parameters Employed In GI Model 
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Initial 
Volumes 
(ml) 


Surface 
Area 
(cm^) ' 


Average 
Transit 
time 
(hr)' 


Volume 
Transfer 
Rates (t9o) 
(hr')' 


New Water 
Absorption 
Rates 
(hrV 


Stomach 


1.0-2.5 


100 


NA 


0.5-3.0 


0.8-4.6 


0 


Duodenum 


4.0-6.4 


0 


147-168 


0.20-0.25 


9.2-11.5 


0 


Jejunum 


4.4-6.4 


0 


913.5-1044 


1.0-2.0 


1.15-2.3 


4.0-4.5 


Ileum 


6.8-7.4 


0 


913.5-1044 


1.2-1.5 


1.54-1.9 


2.4-2.7 


Colon 


5.5-7.0 


0 


763-872 


18-36 


0.064-0.13 


1.4-1.6 



Lui et al. J Pharm Sci 1986;75(3):271.4; Youngberg et al. Dig Dis Sci 
1987-32(5)-472-80; Channan et al. J Pharm Sci 1997;86(3):269-82; Langguth 
et al'siopharm Dmg Dispos 1994;15(9):719-46; Kararli TT. Biopharm Drug 
Dispos 1995;16(5):351-80; 

b) Wagner JG. J Pharm Sci 1961;50(5):59-87; Ho NF. Park JY, Ni PF, et al. 
Crouthamel W, Sarapu AC, editors. Animal Models For Oral Drug Delivery 
In Man: In Situ And In vivo Approaches. Washington, D.C. American 
Pharmaceutical Association. 1983; 2. Advancing quantitative and mechanistic 
approaches in interfacing gastrointestinal drug absorption studies in amraals 
and humans, p. 27-106; 

c) Ho et al. Crouthamel W, Sarapu AC, editors. Animal Models For Oral Drug 
Delivery In Man: In Situ And In vivo Approaches. Washington, D.C. 
American Pharmaceutical Association. 1983; 2. Advancing quantitative and 
mechanistic approaches in interfacing gastrointestinal drug absorption studies 
in animals and humans, p. 27-106; Oberle et al. Journal of Pharmacokinetics 
& Biopharmaceutics 1987;15:529-44; Davis SS. S TP Pharma 1986;22:1015- 
22; Davis et al. Gut 1986;27:886-92; 

d) Tumberg LA. Digestion (1973) 9:357-81. 



2. Adjustment Parameters 

Differences between in vitro and in vivo conditions, as well as differences between in 
vivo conditions for one species of mammal and a second hamper accurate prediction of 
absorption using a simulation approach. For example, in vitro dissolution rate may or may 
not be comparable to dissolution rates existing in vivo, or, the permeability in rabbits may or 
may not be comparable to the permeability in humans. 
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To compensate for such differences, a set of selectively optimized adjustment 
parameters were developed. These parameters were designed to be utilized as constants that 
modify the underlying equations of specific compartments of the Gl model to pemiit 
automatic correlation of input data to output data as well as facilitate accurate prediction of 
oral absorption for a diverse set of compounds. For example, the differential equation 
utilized to calculate fluid volume absorption/resorption employs a rate constant obtained from 
an equation that is modified by a volume absorption adjustment parameter Zva (see Eq. 11). 
Listed below (Table 18) are examples of parameters that can be used to adjust parameters 
and equations as well as those which can be added or removed to a given model if necessary. 
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Table 18: Adjustment Parameters 



PCT/LIS99/2U51 



Compartment 


Segment 


Regional fluid absorption 


stomach 

duodenum 

jejunum 

ileum 

colon 


Flux/Penneability 


duodenum 
jejunum 
ileum 
colon 


Active/Carrier mediated 
Transport (absorption) 


duodenum 
jejunum 
ileum 
colon 


Compound Efflux (secretion) 


duodenum" 
jejunum 
ileum 
colon 


Transfer rates 


stomach to duodenum 
duodenum to jejunum 
jejunum to ileum 
ileum to colon 
colon to waste 


Surface Area 


duodenum 
jejunum 
ileum 
colon 



The adjustment parameters were developed and optimized using a stepwise selective 
optimization process. Initial adjustment parameters were developed for correlation between 
humans and rabbit as follows. Two primary sets of data were used: 1) FDp and best fit 
plasma profiles fiom in vivo clinical pharmacokinetic (PK) data, and 2) simulated FDp and 
plasma profiles generated fmva the GI model. The FDp and best fit plasma profiles from in 
vivo PK data was obtained by analyzing and processing IV and PO data from humans for the 
test set of compounds described in Example 2 using a regression-based curve fitting 
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algorithm to determine the best fit curve that matched the actual cUnical plasma profiles. The 
second set of data was generated using a developmental GI model. 

In vitro data (pemeability, solubility, dissolution rate, and dose) were used as inputs 
into the GI model with the adjustment parameters set to some initial value previously 
determined to provide reasonably predictable values for FDp. The GI model was used to 
provide FDp data for each test compound. The FDp data generated from the GI model also 
was used as input data into an IV/PO PK model, such as the one shown in Figure 19, to 
determine plasma profiles. 

The PO input to the IV/PO PK model of Figure 19 used for fitting clinical data is an 
error function and shown in Equation 17. 



Where D is the dose of drug delivered to the intestine, t is time in minutes, t50 is the time for 
50% of the drug to be absorbed, and Pe is a parameter (Peclet number) related to the slope of 
the linear portion of the absorption curve. 

When fitting the data, all available in vivo PK data (multiple intravenous (IV) dosing 
and multiple oral (PO) dosing) was analyzed simultaneously using the IV/PO PK model of 
Figure 19. The data were weighted by 1/Standard Error of the Mean (SEM) or by 
1 /Concentration^. 

The initial adjustment parameter values were determined empirically. Using a limited 
set of compounds and corresponding in vitro data from rabbit tissue, the adjustment 
parameters were manually varied to obtain FDp values that were reasonably consistent with 
the actual PK data. After the initial values were determined, the GI model developed using 
STELLA® was converted to a program file readable by a program having fitting algorithm, 
such as KINETICA™. The initial adjustment parameters were then simultaneously fit using 



1-1 



JO 




(1) 



(Eq. 17) 
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non-linear regression analysis in a stepwise manner to determine optimized values (i.e., best 
fit values) for the adjustment parameters. Within each step, a few parameters were selected 
for optimization by simultaneous fitting. The fitting was approached using an iterative 
process, where selected adjustment parameters were varied systematically such that the 
deviation of the GI model determined absorption from the actual PK determined absorption 
was minimized. Once the deviation was reduced to a satisfactory level, few more parameters 
were then selected and optimized. The process was continued until all parameters were 
successfully optimized. The new parameters were then placed into the GI model and the FDp 
determined for each compound which is compared to the PK FDp to establish the goodness 
of fit. This process was repeated until an acceptable goodness of fit was established. Using 
this approach, adjustment parameters were developed to correlate, for example, in vitro 
solubility, dissolution, dose and permeability in rabbits to in vivo human absorption. 
Although FDp was employed as the reference for deviation, the actual measurement of 
absorption can be evaluated using any number of parameters, such as plasma levels, 
absorption constants, or others. Moreover, it will be appreciated that many sets of adjustment 
parameters may be developed and established. For instance, others sets of adjustment 
parameters may be established to correlate dog, rat, monkey or other species permeability 
data to human, dog, rat. rabbit, monkey, or other animal in vivo absorption. 

3. Regional Permeability Correlation Parameters 

Since Pe in all intestinal regions may not be available, for instance when cell 
monolayer data is employed to determine Pe in colon, a correlation was developed that 
provides a reasonable prediction of unknown Pe values in the other intestinal regions. 

An objective was to establish a correlation between regional permeabilities that 
allowed prediction of peimeabiUty in the duodenum, jejunum or ileum using known 
permeabilities in one or two of the other regions. 

Correlation development involved obtention of regional permeability values in 
intestinal tissue from the literature and experimentally using methods consistent with the 
experimental protocols as described in Examples 4-5. 
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The regional correlation parameters are estimated using a polynomial equation 
developed for this purpose (Equation 17). Any regional permeability, can be calculated 
using any number of other provided permeabilities. 

The regional correlation parameter function was incorporated into the GI model using 
a logic function module. A control statement was utilized to regulate activation of the 
regional correlation parameter estimation function when a user provides less than the total 
number of permeability values for the segments of the GI tract. 

The following (Table 19) shows correlations that were established along with the 
corresponding correlation coefficient. Correlations were accomplished by data 
transformation and fitting to a non-linear fiinction. 



Table 19: Results of Regional Correlation 





Variable 


Correlation CoefTicient 


Dependent 


Independent 




Duodenum 


Jejunum 


0.870 


Duodenum 


Ileum 


0.906 


Jejunxim 


Duodenum 


0.858 


Jejunum 


Ileum 


0.914 


Ileum 


Duodenum 


0.855 


Ileum 


Jejunum 


0.894 



As an example of the capability of the correlation, two of the above correlations were 
evaluated by estimating the penneability in the duodenum and jejunum using ileum Pe 
values. The compounds chosen were those for which complete Pe values were available. 

The error and % error for the permeability calculations were determined by 
comparing predicted values to the known permeabilities (Table 20). 
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The above results demonstrate that the regional correlation parameter function of the 
GI model was able to accurately predict Pe values for compounds within the initial data set 
(i.e., high r^). 



10 



15 



Example 13: Validation and Testing of GI Model 

To demonstrate that the physiological parameters of the model were operating in a 
logical manner consistent with expected behavior in vivo, the parameter were varied and the 
effect on output monitored. For example, a decrease in the surface area available for 
absorption should result in a decrease in the amount of compound absorbed. Thus, the 
physiological parameters of the model were varied by increasing and/or decreasing their 
values The effect of these variations on the rate, as measured by T50 (time for 50% 
absorption), and extent, as measured by FDp. were simulated. Table 21 shows the 
physiological parameters that were varied and the expected effect on FDp and T50. 
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Table 21 : Physiological Parameter Variations* 



Parameter 


Range evaluated 


Expected effect 


Surface Area or 


0.05 to 10 X Normal* 


Increase in: 
Surface Area or 


Increase tup 


Permeability 
GI Transit Time 


1x10'*^ to IslO'^ cm/s 
0.05 to 10 X Normal* 


Permeability 
Increase in: 
GI Transit Time 


Decrease T50 
Increase FDp 

Increase T50 


Dissolution Rate 


0.05 to 10 X Normal* 


Increase in: 
Dissolution Rate 


Increase FDp 
Decreased T50 


Solubility 


1 to lOOmg/ml 


Increase in: 
Solubility 


Increase FDp 
Decrease T50 

.^1- ...» 4Ua Mm 



chosen was varied, all other parameters were held constant. 



5 All effects on FDp and T50 were as expected with the changes in the physiological 

parameters. While not all of the ranges were in the physiological range, the lower part of the 
range was included to assure that the model would limit to zero FDp ar. the various 
parameters approached zero. 

The basic structure of the GI model also was assessed by comparing its ability to 
10 predict, from dose and in vitro solubility and rabbit tissue permeabiUty data, the rate and 
extent of oral drug absorption in humans and dogs for several drugs, including atenolol, 
ganciclovir, verapamil, and napioxin. These compounds were chosen for their well known 
and diverse in vivo absorption properties and interspecies absorption variabihty. By changing 
the physiological parameter values of the simulation model so that they corresponded to the 
15 species under investigation, but not changing the model structure, i.e., compartment, flow 
regulator, converter relationships, efficacy of the model structure could be evaluated. Initial 
parameter values for dog and human were derived from the literature. Adjustment 
parameters were used to build the correlation between the in vitro data and in vivo absorption. 
For all four drugs, the GI model accurately predicted FDp for both dog and human. 

20 To assess the basic power of the GI model for predicting oral drug absorption, the 

model was tested by simulating FDp as a function of time so as to separate absorption across 
intestinal tissue from first pass metabolism and drug concentration in systemic circulation. 
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Accordingly, methods were developed and used to determine FDp from clinical plasma data 
so that transport across the intestinal tissue could be determined. This was accomplished by 
simultaneously fitting clinical pharmacokinetic data (PO and IV) to the two compartment 
open rV/PO PK model illustrated as a compartment-flow model in Figure 19. Elimination 
was from the central compartment. Input from oral doses was into a pre-systemic 
compartment (for metabolism) which was in equilibrium with the central compartment. FDp 
was determined simultaneously for each oral dose. Clinical pharmacokinetic data fitted to the 
IV/PO PK model demonstrated the ability of the model to accurately determine blood levels 
in the central compartment. 

The fitted clinical FDp data for a test set of compounds was then compared to FDp 
predicted by the GI model using both experimental in vitro values for permeability as input as 
well as estimated permeability values calculated by the model using the regional permeability 
correlation function. The permeability source of the test compounds are shown in Table 22 
below. 
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Table 22: Permeability Source of Test Compounds 



Compound 


rermeauiiiiy 
source* 


ocl 


Experimental 


oc2 


Experimental 


oc3 


Experimental 


oc4 


Experimental 


oc5 


Estimation 


oc6 


Experimental 


oclO 


Estimation 


Pl 


Estimation 


P2 


Estimation 


P3 


Estimation 


P5 


Estimation 


P6 


Estimation 



segments were submitted. Estimation - permeability 
values were calculated using regional permeability 
correlation parameters. 



Figures 49-53 are illustrative of the results of these tests. The physiological model 
was found to accurately predict FDp for the test set of compounds. The accuracy of the 
prediction is based on both rate and extent of absorption. Correlation of FDp extent between 
the clinical data and as predicted by the model for the test set of compounds yielded a 
collective regression coefficient (r^) of greater than 0.92. 



Example 14: Smoothing Functions for GI Model 

In the in vivo physiological situation. permeabiUty and pH do not change at distinct 
points or places within the GI tract (with the exception of the gastro-duodenal junction). For 
example, permeabiUty of a drug in the duodenum may be measured at 1.5 x lO"' cm/s and 2.5 
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X 10-* cm/s in the jejunum, but there is no distinct point in the intestine where such an abrupt 
change exist. Since the Gl model simulates five regions or segments of the Gl tract, and each 
segment utilizes its own set of initial permeability and pH values, an abrupt change, as 
opposed to an incremental transition, is simulated for a dosage form or dissolved drug as it 
passes distally through the segmented Gl tract. 

To account for this phenomenon, and to generate a Gl model that is as physiologically 
accurate as possible, smoothing functions were incorporated into the model. Pairs of 
exponential functions were used to adjust the pemieability and pH values in each segment of 
the intestine. The functions were developed to be time/position dependent using the mean 
cumulative transit time as cues for adjustment. For example, prior to the cumulative transit 
time to reach the ileum (CrrI), the ileum permeability will be equal to the user provided or 
regional correlation estimated jejunum permeability. As time approaches C-rrl, the ileum 
permeability will correspond to the exact average of the jejunum and ileum permeability at 
that point. Immediately after CttI, the ileum permeability continues to gradually 
decrease/increase exponentially until it reaches the user provided, or estimated, ileum 
permeability. 

Two exponential functions were used in combination to effectively smooth the 
permeability and pH values. The Gl model was adapted to employ Equation 19 as the time 
approaches a mean cumulative transit time (Crr), and Equation 20 immediately after Crr. 

p = A-ke(") (Eq.l9) 

P = B + ke-°<-'^ (Eq.20) 

Where A = permeability or pH in the previous intestinal region or segment, B = permeability 
or pH in the latter region, k is defined in Equation 21, a = a constant used to determine the 
steepness of the transition between regions and is inversely proportional to the transit time of 
the region, t = time, and TT = cumulative transit time to exit the previous region. 

k = 0.5(A-B)/e°" (Eq.21) 
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These smoothing functions were utihzed to adjust permeability and pH at junctions of 
stomach/duodenum, duodenum/jejunum, jejunumyileum, and ileum/colon. 
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Appendix 1: Abbreviation Key for Mass-Volume M< 

Abbreviation 

Kf sd = associated rate constant for stoma ch and duodenum 
Ka dj = associated rate constant for duodenum and jejunum 
Ka ji = associated rate constant forjejunu m and ileum 
Ka ie = associated rate constant for ileum and colon 
Ka CO = associated rate constant for colon a nd excretion 
SD trans = transfer rate between stomach and duodenum 
DJ trans = transfer rate between duodenum and jejunum 
JL trans = transfer rate between jejunum and ileum 
IC trans = transfer rate between ileum and colon 
Waste = transfer rate between colon and excretion 
pH s = pH stomach 
pH s2 = pH duodenum 

pH s3 = pH jejunum 

pH s4 = pH ileum 

pH s5 = pH colon 

sol profile = solubility profile for stomach 
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sol profile 2 = solubility profile for duodenum 



sol profile 3 = solubility profile for jejunum 



sol profile 4 = solubility profile for ileum 



sol profile 5 = solubility profile for colon 



stom 



ka = associated rate constant for stomach compartments 1 and 2 



duo ka = associated rate constant for duodenum compartments 1 and 2 



Jej ka = associated rate constant for jejunum compartments 1 and 2 



II ka = associated rate constant for ileum compartments 1 and 2 



Colon ka = associated rate constant for colon compartments 1 and 2 



SA stom = surface area of stomach 



SA duo = surface area of duodenum 



SA jej = surface area of jejunum 



SA il = surface area of ileum 



SA colon = surface area of colon 



Perm stom = permeability of stomach 



Perm duo = permeability of duodenum 



Perm jej = permeability of jejunum 



Perm il = permeability of ileum 



Perm colon = permeability of colon 
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Ka sd = associated rate construct for stomach fluid absorption 
Ka du = associated rate construct for duodeunm fluid absorpti 
Ka je = associated rate construct for jejunm fl uid absorption 
Ka il = associated rate construct for ileunm fluid absorption 
Ka CO = associated rate construct for colon fluid absorption 

Note: other abbreviations adhere to above descriptors and ar< 
explanatory 
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Appendix 2: Equations, Parameters and Values For Mass-Volume Model 

amt_plasraa(t) = amt_plasma(t - dt) + (trans_21 + ka - eliminaUon - trans_12) * dt 
INTT amt_plasma = 0 

INFLOWS: 

traiis_21 = k21*comp_2 
ka = tot abs rate 
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OUTFLOWS: 

elimination = amt_plasma*k_elim 
trans_12 = kl2*amt _plasma 

bIood_side_col(t) = blood_side_col(t - dt) + (colon_ka_5) * dt 
P 15 INTT blood_side_col = 0 

i SjcYj = IF Vol_colon*sol_profile_5 >=CoIon THEN 

Colon*SA_colon*perm_colon*3600 
20 Vol_colon*sol_profile_5*SA_colon*penn_colon*3600 
f J blood_side_dou(t) = blood_side_dou(t - dt) + (duo_ka) * dt 

INIT blood side_dou = 0 



3 INFLOWS: ^ ^ 

^ 25 duo ka = IF Vol_duod*sol_profile_2 >= duodenum THtN 

duodenum*SA_duo*penn_duo*3600 ELSE 
^ Vol_duod*sol_profile_2*SA_duo*pemi_duo*3600 
□ blood_side_il(t) = blood_side_il(t - dt) + GLka) * dt 

INTT blood_side_il = 0 

30 

{f^T^ff Vol_ileum*sol_profile_4 >=Ileum THEN Ileura*SA_Il*penn_Il*3600 ELSE 
Vol_ileum*sol_profile_4*SA_Il*perm_Il*3600 
blood_sideJej(t) = blood_sideJej(t - dt) + (Jej_ka) * dt 
35 INIT blood_side Jej = 0 

INFLOWS: *o . ■ •* • ■ ^-^^nn 

Jej ka = IF VolJej*sol_profile_3 >=Jejunum THEN Jejunum*SAjej*permjej *3600 

ELSEVolJej*sol_profile_3*SAjej*permJej*3600 
40 blood_side_sto(t) = blood_side_sto(t - dt) + (stom_ka) * dt 
INTT blood side sto = 0 



ELSE 



INFLOWS: 

stom ka = IF Vol_stom*sol_profile >= Stomach THEN 
45 StoiMch*SA_stom*penn_stom*3600 

Vol_stom*sol_profile*SA_stom*penn_stom*3600 
Colon(t) = Colon(t - dt) + aC_trans - Waste - colon_ka_5) * dt 
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INIT Colon = 0 
INFLOWS: 

IC_trans = ka_ic*Ileum 

5 

OUTFLOWS: 

Waste = ka_col*Colon , * , r-, c -.-r-„i«„ thfm 

colon_ka_5 = IF Vol_colon*sol_profile_5 >-Colon THEN 

Colon*SA_colon*penn_colon*3600 
10 Vol colon*sol_profile_5*SA_colon*perm_colon*3600 
comp_2(t) = conip_2(t - dt) + (trans_12 - trans_21) * dt 
INITcomp_2= 0 

INFLOWS: 
1 5 trans_12 = kl2*amt jslasma 

OUTFLOWS: 

trans 21 = k21*comp_2 

duodenum(t) = duodenum(t - dt) + (SD_trans - duo_ka - DJjrans) dt 
20 INIT duodenum = 0 

INFLOWS: 

SD_trans = if Stomach >0 then kf_sd*Stomach else 0 

25 OUTFLOWS: . , 

duo ka = IF Vol_duod*solj)rofile_2 >= duodenum THEN 
duodenum*SA_duo*penn_duo*3600 ^^SE 
Vol_duod*sol_profile_2*SA_duo*perm_duo*3600 
DJ_trans = ka_dj*duodenum 

30 excretion(t) = excretion(t - dt) + (vol_cw) * dt 
INIT excretion = 0 



INFLOWS: 

vol_cw = Vol_colon*ka_col 
35 excretion_2(t) = excretion_2(t - dt) + (Waste) * dt 
INIT excretion_2 = 0 

INFLOWS: 
Waste = ka_col*Colon 
40 neum(t) = Ileum(t - dt) + (JL_trans - IC_trans - Il_ka) * dt 
INIT Ileum = 0 

INFLOWS: 

JL_trans = kaji* Jejunum 

45 

OUTFLOWS: 

IC trans = ka_ic*Ileum 
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„ . IF Vo, i,eum.solj™me.4 >..l=un> THEN I,eun,.SAJ,.p™J«600 ELSE 
INIT Jejunum = 0 

5 

INFLOWS: 

DJ_trans = ka_dj*duodenum 
OUTFLOWS: 

l^rr-^FHIu"Clj,ro^ 3 >=Jejunum Jeju«um*SAjej*pennJej *3600 

Sosal.coKt) = semsal.coKt - dt) + (Adsorp_col - col_secretion) dt 
INIT serosal_col = 0 

15 

S^coi=PULSE(1.67,0,.l)H-0*Vol_colon*ka_co 
OUTFLOWS: 

::;^tu» -UsaLdouC. - ^ (Ad.o,p.Duo - duo.—) • d. 

INIT serosal_dou = 0 

25 Sj^'uo = PULSE(10.82.0..1H0*Vol_duod*ka_du 
OUTFLOWS: 

duo secretion = PULSE(10.82,0,.l) 

s^lS m(t) = serosaljiKt - dt) ^ (Adsorptjll - ile.secretion) dt 
30 INIT serosal JU = 0 

^oSLn^PULSE(8.83,0,.10)+0*VoUleuni*ka_il 

35 OUTFLOWS: 

lie secretion = PULSE(1. 50,0,1) ^ ^ * j. 

seSeKO = serosalJejCt - dt) -h (Adsorpjej - jej_secretion) dt 

INIT serosal jej = 0 

'° SS!ej=PULSE(.5.76.0,.mO.VoUe!-k.J. 

OUTFLOWS: 

« So« - * (^"-P-^""" - • 

INIT serosal_sto = 0 
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INFLOWS: 

Adsorp_Stom = 0*Vol_stoin*ka_sd 
OUTFLOWS: 

Stom Secretion = PULSE(16.67,0,.1) 
StomichW = StomachO - dt) (- SD.trans - stom.ka) dt 

INIT Stomach = 1000 
DT ITFLOWS ■ 

stoni_Ka - " *-,£ni\ t,Li>E, 
Stomach*SA_stoin*perm_stom*3600 

1 5 INIT total_drug_absorbed - 0 

^ab^'^te = stom ka+duo_ka+Jej_ka+ILka+colon ka_5 
T^EuL.tionii^= Total_EUmmation(t - dt) + (elunination) dt 
20 INITTotal_Elimination = 0 

INFLOWS: 

25 INIT Vol_colon = 0 
INFLOWS: 

vol_ij = Vol_ileum*ka_ic 
col_secretion = 0 

30 

OUTFLOWS: 

vol cw = Vol colon*ka_col 

35 INIT Vol_duod = 0 
INFLOWS: 

vol sd = kf.sd*Vol_stom 
duolsccretion = PULSE(10.82,0..1) 

40 

OUTFLOWS: 

voil di =Vol duod*ka_dj , ^ j., a 

45 lNrrVol_ileum = 0 
INFLOWS: 
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vol ii = VolJej*kaJi 
ile_secretion = PULSE(1.50,0,-l) 

5 S'^U = PULSE(8.83A.10K0*VoLileum*ka.il 

INTT Voljej = 0 

10 INFLOWS: 

veil di=Vol_duod*ka_dj 
jejJecretion = PULSE(2.67,0..1) 

OUTFLOWS: 

INlf Vol_stom = PULSE(8.33,0,.l) 
^m°s^retion = PULSB(l6.67,0,.l) 



25 



OUTFLOWS: 
vol_sd = kf_sd*Vol_stom 



Adiorp Stom = 0*Vol_stom*ka sd 
conc_prasma= (amtj)lasma/volume)*mg_to_ug 



kl2 = .839 
k21 = .67 
ka_co = 1 
30 ka_col = 3 
ka_dj = 3 
ka_du = 1 
ka_ic = 3 
kajl = 8.83 
35 ka Je = 1 
kaji = 3 
ka_sd = 1 
kf_sd = 2.8 
k_elim = .161 
40 mg_to_ug=1000 

penn_colon = 3.80e-6 
penn_duo = l.lOe-6 
penn_Il = 4.066-6 
perm Jej = 2.17e-6 . 
45 perm_stom=1.10e-6 
ph_s = 1.5 
ph_s_2 = 6.6 
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ph_s_3 = 6.6 
ph_s_4 = 7.5 
ph_s_5 = 6.6 
SA_colon= 138 
5 SA_duo = 125 
SA_I1=102 
SAjej = 182 
SA_stom = 50 
volume = 4* 19200 

10 sol proble 

= GRAPH(ph.s) nofl. 4.00). (3.50, 3.80), (4 00. 3 

fS: fss) M 65) S S («.50: 3.65), (7.00, 3.65), (7.50, 3.65), 
S'^S: 3'.6^t(S^S^ or5m'(«0, U.0), (.0.0, 23.5) 

soljirofile_2-OR^H&.M-^) (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), 

K ?5S' [50 : "6 , ;fo; IS S 3S '(6.0 3.65), (7.00, 3.65), (7.50, 3.65), 
(tw, 3.65);(?.50,'4.00);(9.00, 5.00). (9.50, 12.0), (10.0. 23.5) 

sol_piofilO = G'^H(ph sj> 4 00), (3.50, 3.80), (4.00, 3.65), 

life- 3 50} '(5 oo- 3 65?; * am 3:°°); poo. ^•'^>- 

m St 00; 3.6U (8 50°-^^l), (9.00, 5.00), (9.50. 12.0). (.0.0, 23.5) 

rU solj>roBte_4 = GRAPHJph sjl) „ „(, 4 (3.50, 3.8O), (4.00, 3.65), 

n To- f52' (5 00 6 ) 6 ; toS; 3.6?}, l^o! 3.65!, (7.00, 3.65), (7.50, 3.65), 
S 25 lt.o°o;3.«t(K?(loO,5.00),(9.5^ .2.0), (.0.0, 23.5) 

rij sol_profile_5 = GRAPH(ph s 5) 3 gO), (4.OO, 3.65), 

a ilfo' f5S' J5 00 3 fo; "a ; m 3S 6.S: 3;65?; J7.00, 3.65), (7.50, 3.65), 

g P, 3.iS:(8 50°'«0% 00, 5.00,^50, .2.0). (10.0, 23.5, 
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Appendix 3: Abbreviation Key For GI Model 

The legendAcey has been divided into sub-sections corresponding to the sub-sect.ons of the 
5 model diagram. 

u . fw.n 2 3 4 5 6)havebeenassignedto distinguish between intestinal 

10 1 - stomach 

2 - duodeniun 

3 - jejunum 

4 - ileum 

5 - colon 
15 6 -waste 

3, the jejunum, etc. 

Ghos« «e listed »nder ^ sub-section comatog A= original reservoir, flow regulalor. or 

converter. 

Abbreviations listed in italics are regionally dependent and set up as arrays to allow 
25 independent values for each intestinal region. 

In general ADJ as a prefix indicates a calculated par^eter value (ADJ = adjusted), while 
^Tas a ^iffix indicates an adjustment parameter (ADJ = adjustment). 
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Intestinal model 



f . m?mirf /^""^partmeDts 

VOLABS Fluid volume absorbed 

1% iSg conuined «i* a formulata or comroBed release 

Sublemassofdrug(notcontainedwithintheformulationor 

controlled release device) 
„QT Soluble mass of drug 

ABSORPTION Mass of drug absorbed 

45 r^fpilators 

j^g^g Rate of water absorption 
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VOL OUT 
CROUT 
CR INPUT 
MASS OUT 
DISS PRECIP 
SOL OUT 
FLUX 
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Fluid volume transit rate 

Formulation or controlled release device transit rate 

Drug release rate from formulation or controlled release device 

Insoluble drug mass transit rate 

Dissolution rate 

Soluble drug mass transit rate 

Absorption rate 
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ADJ PARMS (Adjustment Parameters) 
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VOL ADJ 
DISS ADJ 
TRANSIT ADJ 
SA ADJ 
FLUX ADJ 
EFFLUX ADJ 
CARRIER ADJ 



PARMS (Parameters) 

VOL FARM 
SURFACE AREA 
DOSE 

INIT VOLUME 
TIME IN HOURS 
pH 

PARACELLULAR 



Fluid volume absorption adjustment parameter 
Dissolution rate adjustment parameter 
Transit time adjustment parameter 
Surface area adjustment parameter 
Passive Absorption adjustment parameter 
Efflux or secretion adjustment parameter 
Active absorption adjustment parameter 



Fluid volume absorption rate constant 
Surface area available for absorption 
The administered dose of drug 
The administered volume of water or fluid 
A clock 

The physiological pH value 

A user controlled switch used to adjust absorption based on 
absorption mechanism 



TRANSIT TIME 

35 

TRANSFERS 
CUMUTT 

ADJ TRANSIT TIME 
40 USER TT INPUT 



GI transit rate constant 
Cumulative transit time 

Adjusted GI transit time incorporating adjustment parameter 

and user input _ . 

User controlled adjustments to the GI transit time 
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OUTPUT CALCULATIONS 

ABSORBED TOTAL Total mass of drug absorbed (sum of ABSORPTION 1...5) 
ABSORBED luim. potion or the dose absorbed into portal vem x 100 

FLUX TOTAL Total absorption rate (sum of FLUX 1...5) 
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rTTTwrnTqs Cumulative drug mass dissolved 

r^SlTe Cumulative drug mass released from formulation 

PERMEABILITY CALCULATION 

AnT PFRM Adjusted permeability incoiporating all transport mechanisms 

and relevant adjustment parameters 
^ Active or carrier-mediated absorptive permeability 

^ Constant from the Michaelis-Menten type permeability equation 

for active fransport , , , 

T,T:r^TOMAT Passive permeability after regional correlation calculation (same 

RbUiuiNAi. p^gg if regional correlation is not used) 

DAccPF Passive permeability entered by user 

PASS PE ^ ^^^^.^^ determimng the regional correlation 

A logical function used in determining the regional correlation 



SOLUBILITY CALCULATION 

USER oH User supplied pH value for which a solubility value is av^ii^^^^ 

SsER SOLUB User supplied solubility value corresponding to the USER pH 

Am QHT T m Solubility calculated (if necessary) at the appropriate pH value 

ADJ SOLUB ^.^^ the entered USER pH and USER SOLUB values 

30 CONTROLLED RELEASE CALCULATION 

P The instantaneous release rate fi^m the formulation 

PR DOSE The total dose contained with the formulation 

AT TIME The cumulative drug mass release profile 

35 CR AT LAST The cumulative drug mass release profile 

Note- CR AT TIME holds the value at the current time value (t). CR AT LAST holds the 
value at the immediately preceeding time value (t-1) 
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CONC CALCULATION 

CONCENTRATIONS The dissolved drug concentration 
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DISSOLUTION CALCULATION 

Precipitation rate constant 

ADJ DISS PKbCir ^.aicuiated concentration 



114 



PCT/US99/2n51 

Appendix 4: E<,ua.i.ns. Para-neurs and V,l»« F.r Gl MoW 
(7=*- AOJ FARMS 

^ O CARRIER_ADJ[C0MPS1 = 0 
O D1SS_ADJ[C0MP_11 = 1 
Q DISS_ADJ[C0MP_21 = 1 
Q DISS ADJlCOMP_31 = 1 
O DISS~ADJlCOMP_41 = 1 

O dissIadj(C0MP_5] = 1 

Q EFFLUX_ADJ[C0MPS] = 1 
O FLUX ADJ[C0MP_1] = 1 
O FLUX ADJIC0MP_2] = 10 
Q FLUX ADJtCOMP_31«8 
O FLUX_ADJIC0MP_4] = 2 
Q FLUX_ADJ(C0MP_51 = 1 
O SA ADJICOMP.II = 1 
O SA A0J(C0MP_21=1 

O SA ADJ[C0MP_31s1 

Q SA_ADJ(C0MP_4)»1 

O SA ADJ[C0MP_51 = 1 

O TRANSIT_ADJ(C0MP_11=1 

n TRANSIT AOJ(COMP_21=1 

O TRANSIT ADJIC0MP_3]=1 

O TRANSIT ADJlCOMP_4l=1 

O TRANSIT_ADJ[C0MP_51=1 

O VOL.ADJICOMPJI = ■• 

O V0L_ADJtC0MP_2] = 1 

Q VOL ADJ[C0MP_31 = 1 
O V0L~ADJ(C0MP_41=1 
O V0L_ADJIC0MP_51 = 1 
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- CONCENTRATIONSfCOMP 5] - if V0L_5 » 0.0 th«n 0 eta. if (V0L_5<1e^ AND SOL_5<1e-7) 
^ SnTerseTAcS.SOL^ t»^e" ADJ.S0LUBIC0MP_5] else 

lo^slT^^iot 4*S0L 3*SOL_2WOL_3WOL.UVOL_2WOL_4) 
51 CONTROL RELEASE CALCULATION 

' O CR"S?TC = (CR AT_T1ME-CR_AT.LASTr20-(CR_OOSEyi00) 

^ ^''-''l-^l^^'TfilT^?^^) (0 75 42 2) (1.00. 50.6). (1.25. 57.1). (1.50. 62.1M1.75. 
'° Lo nil' ?1 6W2^ 7?^) (2 75 7^9). (3.(^0 76.0). (3.25. 76.9). (3.50. 77.6). 

78° ; 78 5) [4 S S 4.5i ?9 1 . (4.75. ^9.3). (5.00. 79.5). (5.25. 79.6). (S^^J. 
Si t; .li^: 79?; [e m' 79 8) (6 25 79.8). (6.50. 79.9). (8.75. 79.9). (7.00. 79.9). (7.25. 79.9). 
n 80 0; ?7 75 K (8 S 80 0) (8.25. M.O). (8.50. 80.0). (8.75. 80.0). (9.00. 80.0). (975 
in m' K ?9 75 80 0) {10.(^. 80.0). (10.3. 80.0). (10.5. 80.0). (10.8. 80.0). (1 1.0. 80.0). 
01 3. 80 0)". 5. 80 O)'. (118. S) 0). (12.0. 80.0). (12-3. 80.0). (12.5. 80.0). (12.8. 80.0). (13.0. 
80.0)... 

® fn^'^^'Vllwo m'^Ts 31 5) (0 75. 42.2). (1.00. 50.6). (1.25. 57.1). (1.50. 62.1). (1.75. 
(0.00. Om. (0.25 7 7). 3 -^Uy . j' ^ 77 g). 

Ili's 78~)' ?4 S' 78t 4 S 78?: J* 79 1 ) (4.75. !r9.3). (5.00. 79.5). (5.25. 79.6). (5^«) 
7q Is 75 79?) (6 M 79 8) (6.25. 79.8). (6.50. 79.9)..(6.75. 79.9). (7.00. 79.9). (7.25. 79.9). 
It Jo 80 01' ?7 75 MO)' (loo, 80.0) (8.25. 80.0). (8.50. 80.0). (8.75. 80.0). (9.00. 80.0). (9^25 
K ffl S* So 0^ « 75 80 0) ml 80.0). (10.3. 80.0). (10.5. 80.0). (10.8. 80.0). (11.0. 80.0). 
(1°1°3. am. 5. Si O)! ?1 1 8'. io 0)! (12.0. io.O). (12.3. 80.0). (12.5. 80.0). (12.8. 80.0). (13.0. 
80.0)... 

©DISSOLUTION CALCULATION 
n ADJ DISS PREClPlCOMPJ]«ifVOL_1=0ttHin0el8eif 

OL UVOL 2*VOL_3*VOL_4WOL_5) 
O ADJ DISS "PRECIPIC0MP_21« if VOL_2=0 then 0 else « 

OL 1+VOL 2+VOL_3*VOL_4*VOL_5) 
n ADJ DISS PRECIPIC0MP_31 = if VOL_3=0 then 0 else if 

OL_1 ♦VOL_2*V0L_3*VOL_4*V0L_5) 
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r ADJ DlSS_PRECIPtCOMP.41»ifVOL_4=0than0a*i«if 
" (SOL ANOL 4<ADJ S0LUB(C0MP_41) then 

(DISSOUCOMP 4]*DISS.ADJ[COMP_4rMASS_4'{ADJ_SOLUB(COMP_4)-SOL_4/VOL.4)) 

((SOL 4/VOL 4>1aDJ SOLUB[COMP_4irPRECiP(COMP.41 

♦0-(MASS 1 +MASS J+MASS.3+MASS_4^MASS_5*S0L.1 ♦S0L_2+S0L_3+S0L.4+S0L.5*V 
OL UV0L_2WOL 3WOL_4ivOL_5) 
Q ADJ DISS_PRECIP[C0MP_51 = rf VOL_5=0 then 0 else if 
(SOL 5A/0L 5<ADJ S0LUB(C0MP.51) then 

DISSOUCOMP 51'DISS_ADJ(C0MP 5rMASS_5*(ADJ^S0LUB(C0MP.5]-S0L_5/V0L.5)) e*se 
((SOL 5A/0L sPaDJ SOLUB[COMP.5irPRECIPlCOMP_51 

lo-(MASS 1 ♦MASS J*MASS_3*MASS.4*MASS.5+S0L_1 ♦SOL.2*SOL_3*SOL_4*SOL.5W 

OL 1W0L 2+V0L.3+V0L_4+V0L_5) 
Q DISSOUCOMPJl = 1 
Q DISS0L(C0MP_21 = 1 
Q DISS0L(C0MP_31 = 1 
O DISS0L[C0MP_4] » 1 
3 DISS0L[C0MP_5) =« 1 
Cj PREC!P[COMPJ1 = 10 
Q PRECIP{C0MP.21 = 10 
Q PRECIPtCOMP_3] = 10 
Q PRECIP[COMP.4] = 10 
O PRECIPICOMP.S] a 10 
g INPUTS 

(Bl INTESTINAL MODEL 

□ ABS0RPTI0N_1(l) = ABSORPTIONJ(t - dt) ♦ (FLUXJ) • dt 
INITABSORPflONJ «0 

INFLOWS: 

FLUX 1 ^ 

CONCENTRATIONS[COMPjrADJ_PERM(COMP_irSURFACE.AREA[COMP_1] 

□ ABSORPnON_2(t) = ABS0RPTI0N_2(t - dt) ♦ (FLUX_2) * dt 
INIT ABSORpfiON_2 = 0 

INFLOWS: 
^ FLUX 2 B 

^ CONCENTRAT10NS(COMP_2rADJ_PERMlCOMP_2rSURFACE_AREA{COMP_2] 

□ ABSORPnON_3(t) - ABSORPnON_3(t - dt) ♦ (FLUX_3) • dt 
INIT ABSORPfiON_3 ■ 0 

INFLOWS: 

FLU)C 3 " 

^ CONCENTRATIONSICOMP.3rADJ_PERMICOMP.3rSURFACE_AREA[COMP.3) 

□ ABS0RPnON.4(t) = ABSORPTI0N.4(t-dt)*(FLUX_4)*dt 
INIT ABSORPnON_4 « 0 

INFLOWS: 
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^ FLUX 4 * 

CONCENTRATIONS(COMP_4rADJ_PERM(COMP_4rSURFACE_AREA{COMP_4) 

□ ABS0RPTI0N_5{t) = ABS0RPTI0N_5(t ■ dt) ♦ (FLUX_5) * dt 

INIT ABSORPflON_5 = 0 

INFLOWS: 
4» FLUX 5 = if time<32 then 

CONCEIfrRATIONS[COMP_5]*ADJ_PERM[COMP_5rSURFACE_AREA{COMP_5r(32-ti 

mey48*(VOL S/17.2)el8e0 

□ C REL 1(t) = C_REL_1(t - dt) ♦ {- CR_0UT_1 - CRJNPUT.I) * dt 

iriir cJrel_i = cr.dose 

OUTfLoWS: 

H» CR OUT 1 = IF TIME >= CUMU_TTtCOMPJ] THEN C_REL_1*10000 ELSE 0 
5 Cr'iNPUT 1 =ifTIME>CUMU TT[C0MP_1] then 0 else CR_RATE 

□ C REL 2(t) = C_REL_2(t- dt) ♦ (CR_0UT_1 - CR_0UT_2 • CRJNPUT_2) * dt 

inIt cJrel_2 = 0 

"^^°Cr1oUT_1 = IF TIME >= CUMU_TTIC0MP_1] THEN C_REL_riOOOO ELSE 0 
OUTFLOWS: 

^ CR OUT 2 = IF TIME >= CUMU_TTIC0MP_21 THEN C_REL_2*10000 ELSE 0 
^ Cr'iNPIJT 2 = ifTIME>CUMU fTIC0MP_21 then 0 else CR.RATE 

□ C_REL_3(t) « C_REL_3(t - dt) * (CR_0UT.2 - CR_0UT_3 - CRJNPUT_3) • dt 
INIT cJrEL_3 = 0 

INFLOWS: 

^ CR_0UT_2 = IFTIME>» CUMU_TTIC0MP_2] THEN C.REL_2*1 0000 ELSE 0 
OUTFLOWS: 

^ CR.OUT 3 = IF TIME >« CUMU_TTIC0MP_31 THEN C_REL_3*1 0000 ELSE 0 
^ CR INPtir 3»lfTIME>CUMU TTICOMP.S) then 0 else CR.RATE 

□ C ^ 4(t) = C.RiL_4(t-dt)*(CR_OlJr.3-CR_OUT_4-CRJNPUT_4)*dt 

INITcJrEI-_4»0 

'^^°Cr1oUT_3 = IF TIME >= CUMU_TTIC0MP_3I THEN C_REL_3«10O0O ELSE O 

Cr'oUT_4 s IF TIME >» CUMU_TnC0MP_4] THEN C.REL_4'10000 ELSE 0 
^ CR~INPUT 4 » if TIME>CUMU_TnC0MP_41 then 0 else CR_RATE 

□ C REL 5(tl = C_REL_5(t-dt)*(CR_OUT_4-CR_0UT_5-CRJNPUT_5)*dt 

inttcJrel.s^o 

"^^°Cr10UT_4 = IF TIME >- CUMU_TT[C0MP_41 THEN C.REL_4*10000 ELSE 0 

^'iTcT'oUT 5 = IF TIME >= CUMU_TTIC0MP_5I THEN C_REL_5*10000 ELSE 0 
^ CR JNPUT_5 = if TIME>CUMU_TTIC0MP_5] then 0 else CR.RATE 

□ C_REL_6(t) = C_REL_6(t - dt) ♦ (CR_0UT_5) ' dt 
INlT C llEI-_6 = 0 

"*'^°C?0UT_5 - IF TIME >- CUMU.TTICOMP.Sl THEN C.REL.nOOOO ELSE 0 
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3 MASS_1 (t) = MASS_1(t . dt) * (CRJNPUT.I - MASS_0UT_1 - DISS_PRECIP_1 ) * « 
INIT MASS_1 = DOSE 
INFLOWS: 

^ CRJNPUT_1 = if TIME>CUMU_TT{C0MP_11 then 0 else CR_RATE 
OUTFLOWS: 

=» MASS_0UT_1 = MASS_1TRANSFERS{C0MP_1] 
^ DISS PRECIP 1=ADJ DISS_PRECIP[C0MP_11 

□ MASS 2(t) = MASS_2(t - dt) * (MASS_0UT_1 * CR_INPUT_2 - MASS_0UT_2 - 
DISS_PRECIP_2) * dt 

INIT MASS_2 = 0 
INFLOWS: 

^ MASS 0UT_1 = MASS_1TRANSFERS[C0MP_1] 
^ CRJNPUT_2 = if T1ME>CUMU_TTIC0MP_2] then 0 else CR.RATE 
OUTFLOWS: 

^ MASS OUT 2 = MASS_2'TRANSFERS(COMP_21 
J DISS PREcTp_2=ADJ DISS_PRECIPtC0MP_2] 

□ MASS 3(t) = MASS_3(t - dt) + (CR_INPUT_3 ♦ MASS_0UT_2 - MASS_0UT_3 - 
DISS_PRECIP_3) • dt 

INIT MASS 3 s 0 

INFLOWS: ^ 

m OR iNPUT_3 = ifTIME>CUMU_TTtCOMP_3] then 0 else CR.RATE 

^ MASS_0UT_2 = MASS_2TRANSFERStCOMP_2] 
OUTFLOWSP 
^ MASS 0UT_3 = MASS_3*TRANSFERStCOMP_3] 
^ DISS PRECIP 3»ADJ_DISS_PRECIP[COMP_31 

□ MASS_4(t) ="MASS_4rt • dt) ♦ (CRJNPUT_4 ♦ MASS_0UT_3 - MASS_0UT_4 - 
0ISS_PRECIP_4)*dt 

INIT MASS 4 s 0 

INFLOWS: ^„ 
^ OR INPUT_4 « if TIME>CUMU_TTIC0MP_41 then 0 else CR.RATE 
^ MASS_0UT_3 » MASS_3TRANSFERStCOMP_31 

OUTFLOWSP 
^ MASS OUT 4 = MASS_4'TRANSFERS(C0MP_4) 
^ DISS PREcTp 4»ADJ_DISS_PRECIPIC0MP_41 

□ M ASS_5(t) -"mASS.S^I . dt) ♦ (CRJNPUT.5 ♦ MASS_0UT_4 - MASS_0UT_5 - 
DISS PREClP_5)'dt 

INIT MASS_5 ■ 0 

INFLOWS: ^„ 

OR INPUT 5 = ifTlME>CUMU_TTlCOMP_5]then0el8eCR_RATE 

^ MASS_0uf_4 = MASS_4*TRANSFERSIC0MP_4] 
OUTFLOWSP 

MASS OUT 5 = lftime>4thenMASS_5'TRANSFERS[COMP_51elseO 
S DISS PREcTp_5 = ADJ_DISS_PRECIPtCOMP_51 

□ MASS_6(t) sMASS.elt ■ dt) ♦ (MASS_0UT_5) * dt 
IN!TMASS_6»0 

INFLOWS: 
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H» MASS OUT 5 ■ if tifm>4 th«n MASS_5TRANSFERS(COMP_5) ciM 0 

□ S0L_1(t) » SOL.Kt' dt) ♦ (DISS_PRECIP_1 • S0L_0UT_1 - FLUX_1) • dt 
INIT S0L_1 = 0 

INFLOWS: 

^ DISS_PREC1P_1 - ADJ_DISS_PRECIP[C0MP_11 

OUTFLOWS: 
^ S0L_0UT_1 =S0L_1TRANSFERS[C0MP_1] 

C0NCENTRAT!0NS[C0MP_1]*ADJ_PERM[C0MP_irSURFACE>R^^ 

□ S0L_2{t) = S0L_2(t - dt) ♦ (S0L_0UT_1 + DISS_PRECIP_2 - S0L_0UT_2 - FLUX_2) * dt 
INITSOL_2 = 0" 

INFLOWS: 

^ SOL 0UT_1 = S0L_1'TRANSFERStC0MP_1] 
^ D!SS_PRECIP_2 = ADJ_DISS_PRECIPIC0MP_2] 

OUTFLOWS: 
^ S0L_0UT_2 = SOL_2TRANSFERSlCOMP_2] 

i» FLUX 2 s 

^ CONCENTRAT10NSICOMP_2rADJ_PERMtCOMP_2rSURFACE_AREAtCOMP_2) 

O S0L.3(t) « S0L_3(t - dt) ♦ (DISS_PRECIP_3 ♦ S0L_0UT.2 - S0L.0UT_3 - FLUX_3) • dl 
INITSOL_3 = 0 

INFLOWS: 

^ DISS PRECIP 3*ADJ_DISS_PRECIP(COMP_31 
^ SOl3uT_2 trsOL_2'TRANSFERS(COMP_21 

OUTFLOWS: 

SOL_OUT_3«SOL_3TRAN8FERS[COMP_31 

Hfc^ FLUX 3 * 

CONCENTRAT10NSICOMP_3rADJ_PERM[COMP_3rSURFACE_AREA(COMP_3J 

□ SOL 4(t) = S0L_4(t - dt) ♦ (DiSS_PRECIP_4 ♦ S0L_0UT_3 - S0L_0UT_4 - FLUX_4) * dl 
INITSOL_4 = 0~ 

INFLOWS: 

i» DISS PRECIP 4 = ADJ_DISS_PRECIPIC0MP_41 
^ SOLrOUT_3=~SOL_3TRANSFERSlCOMP_3] 

OUTFLOWS: 
^ SOL_OUT_4 = SOL_4TRANSFERS(COMP_4] 

C0NCENTRATI0NStC0MP_4rADJ_PE 

□ SOL 5(t) = SOL_5(t-dt) + (DISS_PRECIP_5*SOL_OUT.4-SOL_OUT_5-FLUX_5)-dt 
INIT S0L_5 = 0 

INFLOWS: 
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^ OISS^PRECIP.5 « ADJ_DISS^PRECIP(C0MP_5] 
^ S0L_0UT^4 = SOL_4*TRANSFERS(COMP_4] 
OUTFLOWS: 

S0L_0UT_5 = if time>4 then SOL_5*TRANSFERS[COMP^5I else 0 
y FLUX.5 = if time<32 then 

C0NCENTRATI0NS{C0MP_5rADJ^PERM(C0MP_5]'SURFACE AREA(COMP 5r(32-ti 
mey48*(VOL_5/1 7.2) else 0 

□ SOL 6(t) = S0L_6(t - dt) ♦ (S0L_0UT.5) * dt 
INIT SOL 6 = 0 

INFLOWS: 

^ S0L_0UT_5 = if time>4 then SOL.5*TRANSFERS(COMP_5] else 0 

□ VOL 1(t) = V0L 1(t.dt)+(-REABS_1 -VOL OUT 1)«dt 
INIT VOL 1 = InTt.VOLUME 

OUTFLOWS: 
^ REABSJ =V0L_1*V0L_PARM[C0MP_1J 
VOL.OUTJ =VOLJ*TRANSFERS[COMPJ) 

□ V0L_2(t) = V0L.2(t - dt) + (VOL.OUT^I - V0L_0UT_2 • REABS.2) * dt 
INIT V0L_2 = 0 

INFLOWS: 

^ V0L_0UT_1 = V0L_1*TRANSFERS(C0MPJ1 
OUTFLOWS: 
y V0L.0UT_2 = VOL.2*TRANSFERS(COMP.21 
^ REAis_2 = VOL_2*VOL_PARM(GOMP_2] 

□ V0L.3(t) = VOL.aO - dt) ♦ (V0L_0UT.2 - VOL.OUT.S - REABS_3) ' dt 
INIT VOL_3 = 0 

INFLOVVS: 

•gp V0L_0UT_2 = VOL.2TRANSFERS[COMP.2) 
OUTFLOWS: 
^ V0L.0UT_3 « VOL.3*TRANSFERStCOMP.3I 
y RE/^S.SsVOL.S-VOL.PARMpOMP.S] " 

□ V0L_4(t) - V0L_4(t - dt) ♦ (V0L_0UT_3 - V0L_0UT_4 - REABS.4) • dt 
INIT V0L_4 = 0 

INFLOWS: 

^ V0L_0UT.3 = VOL.3TRANSFERSICOMP.3] 
OUTFLOWS: 

V0L_0UT.4 H VOL_4*TRANSFERS[COMP_4] 
^ REAis.4»VOL_4*VOL_PARM[COMP.41 

□ VOL 5(t) = VOL 5(t - dt) ♦ (V0L_0UT_4 - V0L.0UT.5 - REABS.5) * dt 
INIT V0L_5 a 0 

INFLOWS: 

^ V0L_0UT_4 a VOL_4TRANSFERS(COMP_41 
OUTFLOWS: 
^ VOL_OUT_5 = VOL_5TRANSFERStCOMP_5] 
^ REABS.5 = VOL.5*VOL_PARM[COMP.5J 

□ VOL_6(t) = VOL_6(t-dt) + (VOL_OUT.5)*dt 
INIT VOL.e ■ 0 
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INFLOWS: 

^ V0L_0UT_5 a VOL^5TRANSFERS(COMP_51 
rj V0L_ABS_1 (t) = >7oL_ABS J (t - dt) * (R£ABS_1 ) * dt 
INIT VOlJaBS.I = 0 
INFLOWS: 

REABS_1 =V0L_1*V0L_PARM(C0MPJ1 

□ V0L.ABS.2(t) = V0L_ABS_2{t • dt) + (REABS_2) * dt 
INIT V0lJaBS_2 = 0 

INFLOWS: 

<^ REABS 2 = VOL_2*VOL_PARM(COMP_21 

□ V0L_ABS_3(t) = V0L_Ais_3(t - dt) ♦ (REABS_3) * dt 
INIT V0lJaBS_3 = 0 

INFLOWS: 

^ REABS 3«VOL_3*VOL_PARMtCOMP_31 

□ V0L_ABS_4(t)»V0L.Ais.4(t-dt) + (REABS_4)*dt 
INIT VOlJaBS.4 = 0 

INFLO\a/S: 

REABS 4 - VOL^4'VOL_PARM{COMP_4] 

□ V0L_ABS_5(t) = V0L.Ais_5(t - dt) + (REABS.5) * dt 
INIT V0lJaBS_5 « 0 

INFLOWS: 

^ REABS.5 = V0L_5^V0L.PARM(C0MP_51 
Q MULTI DOSE CALCULATION 
(a OUTPUT CALCULATIONS 

□ CR.ReleaseW = CR.Relea8e(t - dt) ♦ (CR.cumrate) * dt 
INlfcR_Release = 0 

INFLOWS: 

^ CR.cumrate « CRJNPUT.1K:RJNPUT.2+CRJNPUT.3+CRJNPUT_4+CR.INPUT_5 

□ CUM DISS{t) = CUM_OISS(t-dt)*(CUMM_DISS_RATE)*dt 
INIT CUM.DISS - 0 

INFLOWS: 

CUMM DISS RATE ^ 

DISS.PRECIP,1*DISS_PRECIP_2*DISS_PRECIP.30ISS_PRECIP.4+DISS_PRECIP 
O ABS0f4eD_TOTAL « ABSORPTION.2^ABSORPnON_3+ABSORPT10N_4>ABSORPTION_5 

O PDp% = ABSORBED TOTAL/DOSEMOO 
Q FLUX_TOTAL « FLUX.2+FLUX.3*FLUX_4*FLUX_5 
g) PARMS " 

O DOSE = 1000 
O INIT.VOLUME a 100 
Q PARACELLULAR « 1 
O pHIC0MP_11 = 1.5 
Q pH(COMP.21»5 
Q pH(C0MP_31 « 6.5 
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Q pH(COMP_41«7 

O pH[C0MP_5] = 6.5 

Q SURFACE AREA(C0MP_11 = if PARACELLULAR=0 then 50*SA_ADJ[COMPJ1 else 
50*SA ADJ[C0MP_1] 

O SURFACE AREA[C0MP_21 = if PARACELLULAR=0 then 1 50*SA_ADJ[COMP_2) else 
7 5*SA ADJtCOMP 2] 

O SURFACE AREA(C0MP_31 = if PARACELLULAR=0 then 1000*SA_ADJ[COMP_3) else 

50'SA ADJ[C0MP_31 ~ ^ . 

O SURFACE_AREA(C0MP_4] = if PARACELLULAR=0 then 1000*SA_ADJ[COMP_4] else 

50*SA ADJfCOMP 4] 

O SURFACE AREA(C0MP_5] = if PARACELLULAR=0 then 850*SA_ADJ[COMP_5] else 

42.5*SA_Abj(COMP_5] ~ 
O TIMEJN HOURS = TIME 
O VOL PARMlCOMPJ]»(rVOL_ADJICOMPJ] 
n VOL~PARM[COMP 21=0'VOL_ADJ(COMP_21 
O VOL>ARM[COMP 3] = 1 .75'VOL_ADJ(COMP_31 
O V0l"pARM(C0MP_41 = 1.75'VOL_ADJ[COMP_41 
O V0L>ARM[C0MP_5] =0.10*VOL_ADJ(COMP_51 
3 PERMEABILITY CALCULATION 
O ACT_PEIC0MPS1 = 10 . 

0. 

0. 

0. 

01 

^ !S"l*EFFLU?'JIo^ 1]))'REGIONAUCOMP_irFLUX_ADJtCOMP_1l*3800*(CARRIER_ 
DJ[COMP_irACT_PEICOMPjr3600/(1*(CONCEhmWTIONS[COMP_1J/(KrnCCOMP_1I)))rO 

^ {Sl*SLi?'ADicOM^ 21))*REGIONAMCOMP_2rFLUX_ADJICOMP_2]'3600+(CARRIER_ 
DJ(COMP_2rACT_PE(COMP_2^3600/(1♦(CONCE^m^ATIONS[COMP_2V(K^^1{COMP_2l)))) 

^ SmEFFLU?^ 31))'REGIONALICOMP_3rFLUX_ADJ(COMP_3r3600*(CARRIER_ 
DJ[COMP_3rACT_PElCOMP_3l*3800/(1♦(CONCE^rrRATIONS[COMP_3I/(Km[COMP_3» 

^ m^'^lS^MJlCOMP 41)rREGIONALICOMP_4rFLUX_ADJ[COMP_41-3600*{CARRIER 
SlCOMP_4rACT_PElCOMPj«r360(W(1*(CONCENTRATIONSlCOMP_4V(KmICOMP^^^ 

^ Si'^FFLU^'aD^^ 51)rREGIONAUCOMP.5rFLUX_ADJ[COMP_5r3600*(CARWER_ 
StCOMP_5rACTiElCX)MFLW3800/ 
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C;^ Km[COMPSl » [1 . 

1 , 
1 . 
1 . 
1] 

r , PASS_PE(COMPS] = (0 . 
1.10E-06. 
2.17E-06. 
4.06E-06 . 
3.80E-06 ] 

O RC[COMP 1] = PASS_PE[COMP_irO 

O RC(COMP 21 = IF PASS_PEICOMP_2]>0 THEN 1 ELSE 0 

Q RC(C0MP"3] = IF PASS PE|COMP_31>0 THEN 2 ELSE 0 

O RC(C0MpI41 = IF PASS_PE[COMP_41>0 THEN 4 ELSE 0 

O RC[C0MP"5] = PASS_PE[COMP_5rO 

O RCSUM = RCfCOMP 21+RqCOMP_31+RC(COMP_4] 

O REGI0NAL(C0MP_1] = PASS_PE[COMP_11*RCSUM*0 

n REGIONALICOMP 2] = if RCSUM=2 then 

(EXP( -9.01 1926 + 2.594378 •L0GN(1/PASS_PEIC0MP_2]) -0.065515 

•L0GN(1 /PASS_PE[COMP_2])^2)r(-1 ) else 

if RCSUMs4 then 

(EXP(.0.369414*LOGN(1/PASS PE(COMP 4lH).23756*LOGN(1/PASS_PEtCOMP.4I)*2-0.00e 

9719'LOGN(1/PASS_PE(COMP_4l)*3))*{-1)el8e 

if RCSUMse then 

0.5'{EXP( -9.011926 ♦ 2.594378 •L0QN(1/PASS_PE[C0MP_3]) -0.065515 
•LOGN(1/PASS_PE(COMP_31)*2))*(-1) 

♦0.5*(EXP( -21.009845 ♦ 4.544238 •L0QN(1/PASS_PEIC0MP_4]) -0.140815 
•L0GN(1/PAS8 PEtCGMP_41)*2))*(.1)el8e 

PASS PE(C0MP_21 
O REGIONALtCOMP 31 a |f RC8UM-1 then 

(EXP( -3.238469 .509131 •L0GN(1/PASS_PEIC0MP_21) -O.022109 

•LOGN(1/PASS_PEICOMP_21)*2))*(-1 ) else 

IfRCSUMMthen 

(EXP{-0.093739T.OGN(1/PASS_PEtCOMP_41)*0.182344*LOGN{1/PASS_PE[COMP_4])*2-0.00 
23631*LOGN{1/PASS_PElCOMP_41)*3))*(-1 ) else 
if RCSUM«Sthan 

0.5*(EXP( -3.238469 ♦ 1.509131 •LOGN(1/PASS_PElCOMP_21) -0.022109 
•L0GN(1/PASS PEICOMP_2])»2))^-1) 

+0.5*(EXP( -15.415683 ♦ 3.543563 •L0GN(1/PASS_PE(C0MP_41) -0.100318 
•L0GN(1/PASS PECC0MP_41)^))*<-1)etoe 
PASS PEtC0MP_3] 
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O REGlONAI4COMP.4)«ifRCSUM»1 th«i 

(EXP( 14.455255 -1.264630 •L0GN(1/PASS.PE(C0MP_2)) ♦ 0.082015 
•LOGN(1/PASS_PE[COMP_2]r2)y(-1 ) else 
if RCSUMs2 then 

(EXP( 11.480333 -0.791109 •L0GN(1/PASS_PEIC0MP.3I) ♦ 0.066063 

•LOGN(1/PASS_PE(COMP_3ir2))*(-1 ) else 

ifRCSUMaSthen 

0.5*{EXP( 14.455255 -1 .264630 •L0GN(1/PASS_PE(C0MP_2]) ♦ 0.082015 
•L0GN(1/PASS PEtCOMP_2])*2)n-1) 

•K).5*(EXP( 1 1 .480333 -0.791 109 *L0GN(1/PASS_PE[C0MP_3]) * 0.066063 
•L0GN(1/PASS PEtCOMP_3])*2))*(-1)el«e 

PASS PEtCOMP_<l 
O REGi6nAL[COMP_51 = PASS_PEtCOMP_Sl ♦RCSUM'O 

SOLUBILIY CALCULATION 

O ADJ_S0LUB[C0MP_1) = if USER_pH(C0MPJl>=pH[C0MP_1] then USER_S0LUB[C0MP_11 



((USER SOLUBICOMP 21-USER_SOLUB[COMP_1iy{USER_pH[COMP_2J-USERj)H(COMP.1] 
))*(pH[COMP 11.USER_pHIC0MP_1])+USER.S0LUB[C0MPJl 
O ADJ SOLUB^OMP 2] = If USER_pH[COMP_21=pHlCOMP_2]then USER_S0LUBICGMP_21 
elself USER_pH[COMP_21<pHtCOMP_21 then 

((USER SOLUBICOMP SHJSER SOLUB[COMP_2iy(USER_pH(COMP_3J-USER_pH[COMP_2) 
)r(pHICOMP 21-USEr3h[COMP_21)*USER_SOLUBICOMP_21 else 

((USER S0LUB[C0MP 2}-USER S0LUBlC0MP_1)y(USER_pH[C0MP_2hUSER_pHtC0MP_1) 
)npH[COMP 2]-USERj)HlC0MPJ])*USER_S0LUBlC0MP_1J 
O ADJ S0LUBP0MP_31 = If USER_pHICOMP_3J=pHICOMP_31 then USER_S0LUBIC0MP_3] 
else1iFUSER_pH[c6MP 3]<pH[COMP_3] then 

({USER SOLUB[COMP 4J.USER SOLUBICOMP_3jy(USER_pH[COMP_41.USER_pH[COMP_^ 
)r(pH(COMP 31-USERj)HlC0MP 3])HJSER_SOLUBlCOMP_3J else 

((USER SOLUBICOMP 3hUSER SOLUBJCOMP 2iy(USER_pHICOMP_31-USER_pHICOMP_21 

)r(pH(CbMP 3}.user3hicomp_2ihjser_solubicomp_2i 

O ADJ SOLUBPOMP 4]»lfUSER_pHIC0MP_4]=pH[C0MP_41thenUSER_S0LUBlC0MP_4) 
elself USER_pH[COMP_4]<pHlC0MP_4] theii 

((USER.SOLUBICOMP 5J.USER_SOLUB[COMP_4iy{USERj>HICOMP_5hUSER_pHICOMP_4J 
))*(pH{COMP 4hUSER3HlCOMP_41>*USER_SOLUBlCOMP_4) else 

((USER SOLUBICOMP 41.USER SOLUBtCOMP_3]y(USER_pHlCOMP_41.USERj)HtCOMP_31 
))*(pHICOMP 4HJSEr3hICOMP_31)*USER_SOLUBICOMP_3J 
O ADJ SOLUBICOMP 51»ifUSERjiHtC0MP_31"pHlC0MP_3JthenUSER_S0LUBlC0MP_3] 
elself USER_pHICOMP_3)<pHIC0MP_31 then 

((USER.SOLUBICOMP 4HJSER 80LU8IC0MP_3]y(USERj)HtC0MP_41-USER_pHIC0MP_31 
)r(pHlCOMP 3)-USEr3hICOMP_31)*USER_SOLUBICOMP_3) else 

({USER SOLUBICOMP SJ-USER S0LUBIC0MP_2]y(USER_pHIC0MP_31-USER_pHtC0MP_21 
))*(pHIcbMP_31-USER3HICOMP_21HJSER_SOLUBICOMP_2J 
Q USER _pHIc6mPS] = I1.5. 



else 



5. 
6.5. 
7 . 
7.51 
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Q USER_S0LUB(C0MPS1 « (31 . 
3.65," 
3.65. 
3.65. 
3.65) 
(n» TRANSIT TIME 

O ADJ TRANSIT TIME[COMP 1) = .5TRANSIT_ADJtC0MP_11'USER_TTJNPUT 
n ADjItRANSIt"tiME(COMP 2) = .25-TRANSIT_ADJtCOMP_2rUSER.TT_INPUT 
n ADJ TRANSIt'tiMEICOMP 3ls1.5TRANSIT_ADJ[COMP_3rUSER_TT_INPUT 
n ADJ"tRANSIT"tIME[COMP 4)=1.5'TRANS!T_A0J[C0MP_41«USER_TT_INPUT 
O A0j"tRANSIT"tIME(COMP,5J = 24TRANSIT.AOJ(COMP_5rUSER_TTJNPUT 
n CUMU TTTCOMP 1] = ADJ_TRANSrr_TIME(C0MP_1J 

O CUMuItTIC0MpI21 » A0J_TRANSrr_TIMEtC0MP_1J*ADJ_TRANSIT_TIMEIC0MP_2l 

^ AoIrT^SIT^TiMElCOMPJl^ADJ.T^ 
3) 

^ AoIrTwSsn'^TiiiS 

31+ADJ_TRANSIT_TIME(COMP_4] 

^ AW TwSsn-^TllS'cOMP 1I*ADJ_TRANSIT_TIMEtC0MP_2]*ADJ_TRANSIT_TIME(C0MP. 

31+Abj TRANSIT TIME[COMP_41*ADJ_TRANSIT_TIMElCOMP_5) 

O TRANSFERS(COMP 1] » LOGN(10yAOJ_TRANSIT_TIMElCOMP_1) 

n TRANSFERS[COMP 2) = LO6N(10yADJ_TRANSrr_T1ME[COMP_2l 

O TRANSFERStCOMP 31 » LOGN(10yADJ_TRANSIT_TIME[COMP_3) 

O TRANSFERStCOMP 4) » LOGN(10)/ADJ_TRANSIT_TIMEtCOMP_4J 

O TRANSFERSIC0MP_5J - LOGNdOyADJ.TRANSTT.TlMElCOMP.Sl 

O USER.TTJNPUT"? 
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All publications and patent applications mentioned in this specification are herein 
incorporated by reference to the same extent as if each individual publication or patent 
application was specifically and individually indicated to be incorporated by reference. 

The invention now being fiilly described, it will be apparent to one of ordinary skill in 
the art that many changes and modifications can be made thereto without departing fi-om the 
spirit or scope of the appended claims. 
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